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I.  INTRODUCTION 

The  detection,  localization  and  tracking  or  submarines  by  underwater  acoustic  signals 
are  greatly  disturbed  by  signal  fluctuations.  These  originate  primarily  in  source/receiver 
motion,  but  also  are  caused  by  acoustic  rays  reflected  from  wind-roughened  surface,  and  by 
scattering  from  inhomogeneities  randomly  dispersed  in  the  ocean.  In  past  years  such  fluctua- 
tions have  been  spatially  and  temporally  averaged  to  mean  values.  However,  random  variations 
of  signals  about  mean  values  have  recently  taken  on  significant  meaning  because  they  strongly 
affect  target  detection  methods  which  use  statistical  decision,  such  as  is  found  in  receiver 
operating  characteristic  (ROC)  curves.  In  addition,  a higher  order  statistical  analyses  of  acous- 
tic fluctuations  provides  a means  to  evaluate  the  reliability  of  experimental  data  and  can  lead  to 
the  construction  of  confidence  limits  to  the  mean  values  currently  used  in  ASW  operations. 


Research  progress  in  acoustic  fluctuations  has  been  reviewed  and  evaluated  in  several 
recently  held  conferences  and  workshops.  Since  1974  there  have  been  four  distinctive  but  re- 
lated workshops  of  this  nature,  one  review  panel,  and  an  Acoustical  Society  of  America  (ASA) 
session  held  which  bear  somewhat  on  acoustic  signal  fluctuations.  On  27-28  March  1974  a 
fluctuation  workshop  was  conducted  at  Naval  Research  Laboratory  (NRL)  hosted  by  Acoustic 
Environmental  Support  Detachment  (AESD)  and  a summary  (Ref.  1)  is  available.  An  Interna- 
tional Workshop  (Ref.  2)  on  Low  Frequency  Propagation  and  Noise  was  sponsored  by  CNO 
(OP-95)  and  supported  by  Chief  of  Naval  Research  (CNR)  and  held  at  Woods  Hole  Oceano- 
graphic Institute  (WHOI)  in  October  1974.  The  David  Taylor  Naval  Ship  Research  and 
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. Development  Center  (DTNSRDC)  was  the  host  for  an  Operations  Research  Models  of 
Fluctuations  Affecting  Sonar  Detection  Workshop  during  19-21  March  1975.  Proceedings  were 
published  (Ref.  3)  in  two  classified  volumes.  On  1-2  October  1975  a workshop  (Ref.  4)  under 
auspices  of  Underwater  Sound  Advisory  Group  (USAG)  was  held  on  Acoustic  Signal  Coher- 
ence at  Naval  Underwater  Systems  Center  (NUSC).  Next,  under  direction  of  chief  of  Naval 
Operations  (CNO:OP95)  a technical  review  panel  for  a coordinated  Ocean  Acoustics  Program 
Plan  met  8-10  November  1977  at  Naval  Ocean  Research  and  Development  Acitivity  (NORDA) 
and  produced  a summary  (Ref.  5).  Lastly,  at  the  December  1977  meeting  of  the  Acoustical 
Society  of  America  (ASA)  at  Miami  a session  on  Underwater  Acoustics  was  devoted  to  Fluc- 
tuations and  Signal  processing.  A review  of  several  of  the  significant  papers  presented  at  that 
session  appears  in  Ref.  6. 

Each  of  the  preceding  efforts  did  not  completely  address  the  signal  fluctuation  prob- 
lems. For  instance  the  one  hosted  by  AESD  was  specifically  directed  to  satisfy  potential  SASS 
contractors.  The  International  Workshop  at  WHOI  was  unclassified  and  many  system  applica- 
tions were  omitted.  The  Operations  Research  effort  at  DTNSRDC  neglected  the  environmental 
effects.  The  one  on  Signal  Coherence  at  NUSC  was  constrained  in  subject  material.  The  scope 
of  the  coordinated  plan  at  NORDA  (Ref.  5)  examined  the  contributions  of  Navy  research, 
exploratory  and  advanced  development  programs  to  the  four  ocean  acoustic  subjects  of  bottom 
interaction,  signal  fluctuation,  coherence  and  directional  noise.  Each  of  these  subjects  was 
divided  into  segments  dealing  with  modeling,  measurement  systems  and  data  acquisition  pro- 
grams. Of  necessity  these  elements  received  a very  broad  treatment  and  generally  were  of  a 
non-technical  nature  and  lacked  specificity. 

Lastly,  the  ASA  session  was  instructive  but  again  because  of  its  unclassified  nature  a 
proper  discussion  and  exposition  of  the  efforts  being  made  by  the  exploratory  and  advanced 
development  workers  could  not  be  made. 


c Current  Navy  sponsored  programs  in  data  acquisition  and  analytical  modeling  of 
fluctuations  of  received  acoustic  signals  show  a great  diversity  of  approach,  data  collection,  sig- 
nal processing  and  final  display.  The  ultimate  use  of  the  results  of  research  to  improve  under- 
sea surveillance  often  falls  short  of  expectations  because  of  a mismatch  between  the  user’s 
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understanding  of  his  requirements  and  the  investigator’s  interpretation  of  those  requirements. 

There  is  and  has  been  a need  to  clarify  both  the  user’s  and  investigator's  concept  of  their 
requirements.  To  this  end  the  Chief  of  Naval  Development  with  the  assistance  of  the  Naval 
Electronic  Systems  Command  sponsored  a classified  workshop  on  Acoustic  Fluctuations. 

The  objectives  of  the  Acoustic  Fluctuations  Workshop  were  to: 

a.  Delineate  the  problem  areas  associated  with  modeling  of  fluctuations  of  acoustic 
signals  and  noise  as  caused  by  the  medium,  its  boundaries  and  motion  of  the  source  or  receiver 
platforms.  Assess  the  impact  of  these  fluctuations  on  system  performance  models  as  translated 
into  model  input  parameters. 

b.  ’ Examine  the  current  programs  addressing  underwater  acoustic  fluctuations. 

Present  program  goals,  technical  approaches  to  meet  the  goals,  and  work  currently  underway. 

c.  Recommend  adjustments  to  program  goals  to  meet  the  input  needs  of  users  of 
acoustic  environmental  fluctuation  models. 

On  22-23  Feb.  1978  on  Acoustic  Fluctuation  Workshop  was  held  at  NRL  and  the 
Agenda  was  essentially  as  per  Figs.  1 and  2. 

Attendance  at  the  Workshop  was  limited  to  the  principal  investigators  and  managers 
from  the  involved  sectors  of  the  Navy;  industry  and  academic  community.  Over  50  participated 
in  the  presentations  that  were  made  and  discussions  that  were  held.  An  editorial  summary  was 
made  for  those  papers  which  were  supplied  by  the  authors.  A synopsis  was  prepared  for  all 
unwritten  papers  and  combined  with  a comprehensive  technical  review  on  Ocean  Science  Fluc- 
tuations into  a separate  report  (Ref.  6). 
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Subsequently,  a working  group*  was  charged  with  the  task  of  formulating  objectives, 
approaches  and  specific  tasks  which  could  serve  to  assist  in  structuring  R&D  programs  on 
acoustic  fluctuations.  This  document  is  the  result  of  those  efforts  and  has  undergone  several 
revisions  one  as  a result  of  consultation  and  convening  of  a special  panel**  where  Navy  labora- 
tories and  industry  were  represented. 


‘ l)r  Sam  Danish.  C.R  Rollins.  J.  Cvbulski  of  NRL 
" Dr  P Cable  (NUSCI.  Dr  J Slcwart  (NOSC).  Mr  C Spnflorcl  (SAD 
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INTRODUCTION 

Fluctuations  of  acoustic  signals  propagating  in  the  ocean  is  a subject  of  considerable  com- 
plexity. This  is  caused  primarily  by  random  inhomogeneities  along  the  propagation  path  coupled 
to  the  comparable  size  wavelengths  of  these  signals.  Its  two  key  elements  of  environment  and 
acoustics  have  been  studied  separately  and  together,  and  a copious  literature  exists.  In  this 
regard  the  physics  of  that  part  of  the  ocean  environment  which  is  thought  to  be  the  cause  of 
fluctuations  is  still  receiving  attention.  It  could  be  used  as  the  basis  of  the  review  to  be  made 
below.  An  alternate  and  complementary  basis  is  to  consider  the  propagating  medium  as  a 
“channel”  along  which  a communication  signal  is  being  transferred  from  a source  to  a distant 
receiver.  Such  a communications  approach  offers  a considerable  advantage  in  that  it  permits 
a large  overview  of  the  theoretical  and  experimental  aspects  of  fluctuations  to  be  formulated  in 
a simple  abstract  way  which  makes  visible  the  otherwise  very  complicated  subject  matter  of  the 
study.  We  shall  adopt  the  communications  approach  in  all  that  follows. 

We  begin  then  with  a distant  source  of  sound  and  assume  there  is  a distribution  of  radiators 
(or  transmitting  antenna)  which  generate  a spatial  and  temporal  distribution  of  sound  Q(rs,tg), 
a portion  Qm  of  which  is  injected  into  the  channel,  thus  constituting  the  channel  input.  The 
channel  then  propagates  the  sound  along  a particular  path,  transforming  it  in  the  process,  accord- 
ing to  an  operation  Tm , so  that  it  arrives  at  a receiver  in  the  state  Tm  Qm . At  the  receiver  the 
received  signal  is  then  combined  with  other  signals  in  the  receiver  antenna  which  transforms 
it  according  to  an  operator  TR,  and  delivers  the  total  acoustic  pressure  p = TRTmQm.  During 
transfer,  as  well  as  during  reception,  the  channel  is  corrupted  by  noise  so  that  both  Tm  and 
Tr  contain  random  parts.  The  signal  reported  by  the  receiving  antenna  will  therefore  be  a 
random  variable,  both  because  the  medium  is  changing  in  a random  way,  and  because  the 
medium  and  the  receiving  antenna  are  noisy. 

In  general,  there  will  be  more  than  one  propagation  path  connecting  source  to  receiver. 

We  call  each  distinct  path  a channel.  Let  there  be  a total  of  M channels  between  the  center 
of  the  source  rs  and  the  ith  element  of  an  N-element  receiver,  located  at  rRw.  Thus  the 
total  received  signal  at  the  i*h  element  hydrophone  is, 

P(?r},  tE,  ts)  = £ Tm(ts IE?„7g>,  tE)  Qj>>  (ts\?s  (0,*)  + N(r£>,  tN) 

m 

Here  the  source  radiation  injected  into  the  m^)  channel  is  labeled  , and  is  a function  of 
the  “source  time”  ts,  and  source  coordinates  rs  (6,  <p),  where  0,  are  spherical  angles.  The 
transmission  operators  Tm  which  is  conceived  to  be  a function  of  the  source  time  ts,  and  an 
environment  matrix  E which  contains  all  the  aspects  of  the  environment  which  govern  the  pro- 
pagation of  sound,  both  deterministic  and  random.  In  particular  E is  a function  of  the  “envi- 
ronment time”  tE,  distinguished  from  ts  by  its  (generally)  different  s^ale,  and  a function  of 
source  and  receiver  location.  The  number  of  components  in  matrix  E depends  on  which  trans- 
mission problem  is  to  be  solved,  and  may  include  such  effects  as  randomly  rough  ocean  surfaces. 
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internal  waves,  eddies,  currents,  volume  scattering,  etc.  The  symbol  N is  the  local  noise  at  the 
ith  receiver  hydrophone.  It  is  a function  of  receiver  location  r^  and  of  a time  tN , which  again 
generally  differs  from  tg  and  W by  its  different  scale.  Actually  N represents  an  acoustic  signal 
generated  by  a noise  source  ) at  location  r^  (0,  ip)  and  propagated  along  one  or  more  chan- 
nels to  the  receiver.  Thus,  the  total  received  pressure  field  at  the  ith  hydrophone  may  be  rep- 
resented in  a different  way  by  the  sum  of  multipaths,  each  of  which  may  be  corrupted  by  noise, 

pfc&tB.ts  )=£  Tm  (t8|E(vrN,^,tE))  I Q(^>(t8|7B(0,^))  + Q(^)(tN  |7n  (0.SP)  • 
m=l 

This  formulation  allows  the  mth  channel  to  carry  signals,  or  noise,  or  both.  When  the  receiver 
is  an  L element  array  designed  for  steering,  the  pressure  in  the  ith  hydrophone  is  multiplied  by 
a steering  operator  A(r^,  t).  The  result  is  cross-multiplied  with  the  received  signal  plus  noise 
in  all  other  hydrophones  at  different  time  delays  r d to  form  products  which  in  turn  are  time 
averaged  to  form  hydrophone  cross-correlation  in  space  and  time.  Each  such  cross-correlation 
is  normalized  with  respect  to  the  time  averaged  power  of  the  ith  arid  the  jth  hydrophone, 
thus  delivering  a typical  term, 

A(r^,  rd)  r/t™  rd)  A(?»,  rd) 
lRp(rW,o)Rp(r$,o)]1/2  - 

The  ratio  of  sum  of  all  these  terms  when  signal  only  is  considered,  to  the  sum  when  noise 
only  is  considered,  divided  by  the  same  ratio  in  a single  hydrophone  as  averaged  over  the  array 
constitutes  the  array  gain  : 


A(r(i>,  Td)  Rp*  (r(i£j) , Td)  A(r^),  rd) 
[R(‘>  (r$,o)R(*>  (r^>,o)]1/2 


[R(p)^R),  °)  rp°  1/2 


A(r$,  rd)  R(p}  (r(i£j) , rd)  A(r^,  rd) 

[r(",(*<B.o)r(;)(«*.0)]«2 


R(p ) (r^)»o)] 1/2 


R^  (o) 
R<>) 


In  actual  practice  the  steering  operators  are  random  functions  of  time  and  space, 


A ->  A(r  $,rd  tA) 

in  which  tA  is  an  array  time,  distinguished  from  all  other  times  by  its  scale.  We  return  to  the 
statistics  of  steering  matrices  later. 

The  Environmental  Matrix  if 
— ► 

In  general  E is  a matrix  whose  elements  Ej , E2  ...  describe  the  various  random  environ- 
mental processes  which  affect  the  propagation  of  sound.  Each  element  is  a random  function 
j of  space  (rs,  rR ) and  time  (tE ).  A transmission  channel  can  be  characterized  as  time-varying, 

or  space-varying,  or  both.  To  begin,  let  the  mth  channel  be  a time-varying  channel  randomly 
perturbed  by  Em  (tE ) and  assume  it  has  an  ensemble  of  realizations.  Then  the  autovariance 
of  Em(tE)is, 

KEm  (tE.  W)  = < (Em  (tE)  - < Em  (tE)  > )(Em  (tE)  - < Em  (tfe)  > ) > 

in  which  the  symbol  ( > represents  ensemble-average.  In  real  physical  events  this  function 
effectively  vanishes  then  tE  - tE'  = Lt,  which  then  defines  the  correlation  time  for  the  Em 
perturbation  of  the  environment.  If  the  channel  Em  (tE ) is  nonstationary  then  Lt  = Lt  (tE ). 
Only  for  stationary  channels  is  Lt  function  of  time  difference  tE  - tE,  independent  of  tE. 

Next  let  the  mth  channel  be  space-varying  channel  randomly  perturbed  by  Em(z)  and 
assume  it  has  an  ensemble  of  realizations  in  the  z-direction  (=  depth).  The  auto  variance  in 
spatial  z is  the  ensemble  average 

K(Z,Z')  = < (Em (Z)  - < Em(Z)  > ) <Em(Z')  - <Em(Z') > )> 

Again,  for  real  physical  events  this  function  effectively  vanishes  then  Z-  Z'  ^ Lz  which  defines 
the  correlation  length  of  the  perturbation  Em(Z).  If  Em(Z)  is  nonstationary  then  Lz  = Lz(Z), 
namely  the  correlation  length  depends  on  z. 

Now  assume  all  J elements  Ej  ( E2  . . . Eg  of  matrix  If  participate  in  the  propagation  of 
sound.  This  might  mean,  for  example,  that  the  transmission  channel  experiences  effects  of 
surface  scattering  (call  this  Ex ),  volume  scattering  (E2),  internal  waves  (E3),  eddies  (E4 ), 
tides  (E5 ),  etc.  One  can  then  form  a covariance  matrix  of  all  perturbations  Em , m = 1,  2 . . ., 

K (tE,tE',rg,  rs,  rR,rR)  * TKElEl  KElEz KElE 


KEjEl  KejE2  KEjEj 
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}* 


* -1 
K'  'j 


in  which  the  diagonal  elements  are  autovariances,  already  defined,  and  the  off-diagonal  ele- 
ments are  cross-variances,  either  time-varying  or  space-varying,  or  both.  Their  form  is  illus- 
trated by  an  example  of  a time-varying  channel. 


KE.E.  = < We)  - < Ei(tE)  > ) (Ej (tE')  - < Ej(tg')  > ) 


For  each  element  in  this  matrix  one  can  define  correlation  times  and  length,  as  illustrated 
earlier. 

From  this  it  is  seen  that  a complete  description  of  the  randomizing  effects  of  the  envi- 
ronment requires  the  construction  of  the  covariance  matrix  of  all  elements  Ej , E2 of  the 

environment  matrix  Iti. 

— ► 

In  most  practical  cases  of  the  propagation  of  sound  only  one  realization  of  E is  available, 
rather  than  an  ensemble.  If  this  sample  is  long  enough  (in  time  or  space)  it  can  be  applied  by 
means  of  the  ergodic  hypothesis  to  form  correlation  averages,  R.  Thus  for  a time-varying 
channel,  in  place  of  KEiEj  the  time  average  over  a single  realization  is  used. 


REjEj^E^Ef)  ” ^i^E ) Ej(tE  + tgj) 

in  which  the  overbar  indicates  an  average  over  time  tE . Again,  if  the  process  Ei(tE)js  sta- 
tionary, the  auto-correlation  Re:E:  is  independent  of  tE.  When  all  J elements  of  E partici- 
pate in  the  propagation  of  sound  one  can  describe  the  environmental  effects  by  a coherence 
matrix,  RE 


— > t — ► — — ► — ► r . 

RE  (^E>  ^E>  rs»  rs  » rR » rR ) — 


Re1e1 

RElE2 

■ • • REiEJ 

E-EjEj 

REjE2 

• • • REjEj  _ 

in  which  the  diagonal  elements  are  the  autocorrelations,  and  the  off-diagonal  elements  are  the 
cross-correlations,  all  elements  being  general  functions  of  spt.ce  and  time.  Correlation  times 
and  correlation  lengths  can  then  be  deduced  from  each  correlation  function  in  the  manner 
noted  earlier  for  the  covariance  matrix. 

The  effects  of  the  environment  need  not  be  confined  to  real  space  and  time.  It  is  often 
advantageous  to  transform  the  matrices  R^  and  into  wavenumber  space  and  frequency 
space.  The  results  are  spectral  density  matrices  whose  elements  ‘j^EjEj  are  obtainable  by 
application  of  Fourier  transformation  through  the  Wiener-Kinchine  theorem, 


EjEj  =2f|ke.e.| 


Here  the  transform  coordinate  pair  fixes  which  space  is  being  transformed . The  matrix  of 
spectral  and  across-spectral  densities  corresponding  to  various  perturbations  Em  then  becomes. 
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1E1  ^ExE2 


<E  jEl  °tf& jE2 


c)^E1Ej ' 


^EjEj 


To  illustrate,  let  Ej  be  the  component  of  E which  deals  with  the  propagation  effects  of  the 
temporal  fluctuations  in  sound-speed  profile,  that  is  Ej  -*  6 C(tE ).  Let  the  Fourier  pair  be 
tE  **■  <I>S.  Then,  for  non-stationary  5 C(tE ), 

#sc(^E.^e)  = 2 J K5C(tE,tE')e-il«-EtE  -4-EtEl  d(tE,tE') 

Only  for  a stationary  process  does c depend  on  4>E  - <tE. 

A further  discussion  of  frequency  space,  and  wavenumber  space  will  be  taken  up  in  a 
later  section. 

The  time  coordinate  of  the  environmental  matrix  E has  been  designated  tE.  A second 
time  coordinate  may  be  introduced  whenever  the  space  coordinates  r£,  r^  themselves  become 
functions^  time.  This  occurs  when  source,  or  receiver,  or  both  are  in  motion.  One  writes 
rs(tsM  )>  rR  ( ^rm  )» >n  which  tSM , tRM  mean  source  time  and  receiver  time  associated  with 
motion.  Motion  of  source  and  receiver  thus  associates  a temporal  significance  to  a purely 
spatial  description  of  the  medium.  By  Fourier  transformation  of  tSM , tRM , source  receiver 
motion  can  be  studied  in  frequency  (or  doppler)  space. 


The  Channel  Operator  Tm 

— ► 

While  the  environmental  matrix  E describes  the  effects  of  a randomizing  environment 
on  the  propagation  of  sound,  the  channel  transmission  itself  is  described  by  the  operator 
Tm  which  transfers  the  sound  from  source  to  receiver.  Thus  for  the  mth  channel  which 
contains  source  Qm  the  delivered  pressure  is, 

Pm  - ^m  Qm 

To  elucidate  the  nature  of  Tm  let  Qm  be  delta-function  is  space  and  time, 

Qm  =5(Wso)  5(t-t0) 

It  is  known  that  the  acoustic  pressure  is  then  the  Green’s  function  G: 

Pm  (rR>  Q = Gm  (rR,  rgo  | tg,  tgo) 
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If  we  write 

Pm(?R.  t)  - J J Gm  (i,3  I tg,  t.')  6(5  -50)  6(tg'  - tgo)  d3^  dtg' 

it  is  then  seen  that  the  transmission  is  an  integral  operator, 

Tm  - J |Gm  (?R,  i;  | tg,  t.')  [ ]d3?g'dtg' 

Thus  in  generalTm  is  a weighted  integral  over  space  and  time  of  the  sources  of  the  pressure 
field,  whether  these  sources  are  true  sources  supplying  energy  to  the  channel,  or  fictitious 
sources  redirecting  energy  (by  scattering)  during  channel  transmission.  In  the  former  case 
the  space  of  integration  is  the  small  volume  surrounding  the  true  source,  while  in  the  latter 
case  it  is  the  space  of  the  fictitious  sources  of  the  channel,  such  as  channel  area  in  the  case 
of  surface  scattering,  or  channel  volume  in  the  case  of  volume  scattering,  etc. 

Additional  characterization  of  the  transmission  Tm  is  obtained  by  study  of  the  inverse 
operator  Tm  1 , 

I'm  1 Pm  ~ Qm 

This  is  the  partial  differential  equation  of  propagation  of  an  acoustic  pressure  in  the  presence 
of  sources  Qm.  When  Qm  is  fictitious  it  is  written  Qm  <x  apm. 


The  equation  of  propagation  is  then 


I'm  1 Pm  "a  Pm 

in  which  a is  a function  of  coordinates,  or  is  constant.  Actual  forms  of  this  equation  in  com- 
mon use  are  the  parabolic  equation,  the  Helmholtz  equation,  the  eikonal  equation,  etc.  A 
great  literature  exists  on  the  subject  of  the  operator’s  Tm , Tm_1,  and  the  solutions  of  the  equa- 
tions they  obey. 


The  Acoustic  Pressure  pm 

The  contribution  of  acoustic  pressure  from  a single  channel  to  the  i1*1  hydrophone  of  the 
receiver  is. 

Pm  (i(r  . V *£>  " Tm  (t.  I E(5,  , tE)  Qm  (tg  Ug  (0,  sfi)) 

Since  pm  is  a random  variable  it  can  be  separated  into  a mean  part  and  a fluctuation  part, 

Pm  ~ ^ ^m  Qm  ^ + ® ^m  Qm 

In  nonstationary  fields  the  mean  part  is  a function  of  space  and  time  while  the  fluctuation  part 
is  a function  of  coordinates,  rather  than  differences  of  coordinates.  The  properties  of  pm  are 
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summed  up  in  its  probability  density  PPm  (or  distribution)  and  its  statistics  (or  probability 
moments).  The  statistics  most  often  used  fall  into  the  category  of  single-point  in  which 
moments  are  calculated  for  one  point  in  space  or  time,  or  both;  and  two-point  in  which 
moments  are  determined  with  two  points  in  space  or  time,  or  both.  Three-point  and  higher 
order  of  statistics  are  used  when  special  properties  of  the  acoustic  field  are  needed. 

The  probability  description  of  pm  is  often  extended  from  real  space  to  transform  space 
by  Fourier  transformation.  In  particular,  the  source  time  t8  is  transformed  to  source  fre- 
quency f,  whenever  narrow  (frequency)  band  sources  are  present.  Then  the  field  pressure 
becomes, 

Pm(f..  ^R.  4e)  " Tm  (f.  I E (r„"r+R,  tE)  Qm  (f,  I "r,  (0,  *)) 

which  at  short  enough  ranges,  or  for  weak  scattering,  can  be  represented  by  a phasor  with 
explicit  amplitude  and  phase. 

Pm  (f.»  *R  • 4E  ) = | Pm  (f.»  *R . *E>  | exP  j*  *m  <f»’  *R » lE  )} 

= Pom  (f..^R>  exP(Xm) 

Xm  = Xm  » rR»  ^e)  = Xm/i  \ + * Xm/o\  - 1 


in  which  pQ  is  the  deterministic  acoustic  pressure  at  the  receiver  in  the  absence  of  any  fluc- 
tuations, and  exp  X is  the  effect  of  fluctuations,  as  described  by  a complex  random  variable. 

One-Point  Statistics 


The  received  acoustic  pressure,  whether  expressed  in  real  space  or  transform  space,  or 
mixed  space,  allows  one  to  construct  a list  of  observables,  which  are  defined  as  real  physical 
events  at  single  points  in  space-time;  (See  following  sheet,  Fig.  III.1.) 

Observables  at  single  points  such  as  these  cover  only  single  channel  transmission.  Multi- 
channel transmission  from  M channels  to  single  points  are  described  by  the  expression, 

M ^ 

P(r(jl)  , t„  tE)  = £ T™  (t,  | E(?.,-rW  , tE)  Qm  (t.  | r.  (0,  *) 
m » l 

In  the  case  of  narrow  band  sources  this  becomes, 

P^ . f».  tE)  = Po  (f»>  ) exp  (X) 

A new  list  of  observables  can  be  constructed  from  these  quantities  by  removal  of  subscript 
m from  the  previous  list. 

Two-Point  Statistics 

We  return  now  to  the  case  of  one  or  several  single  channel  receptions  and  consider  two 
space-time  acoustic  pressures  received  under  the  following  conditions: 


I 


p 

<p> 

p2 

Ip  I 
Ip  I2 

0,  x2 

<0),<X2> 

<02>,<x|> 

In  Ip  I2  = 1 
<I2>  = <(ln  Ip  I2)2) 

<01  > = 2 <X!X2) 

Xi 

2 

<p2>  = [exp  2(<X2))  + exp  2(<X2))] 

In  l<p>  I2 

<X2>  = i(<IXl2>  ± Re<X2  >) 

<XiX2>=|lm  <X2> 

In  (p2> 

TQ  cos  0 

Note  the  following  definitions  of  the  channel  pressure  pmi 

Pm  = TmQm  = <TmQm>  + 5QM  = Ipjexp  i 0m  = p0  exp  Xm 

P = STmQm  = Po  exP  x 
X = Xj  + iX2 

m = channel  or  multipath  number 

Figure  III.l 

List  of  Observables  at  Real  Points  in  Space-time 
(Symbols  defined  in  text) 
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(1)  Two  single  channels  to  two  space  points  r^,  at  fixed  time. 

(2)  One  single  channel  at  one  space  point  at  two  different  times  tE,  tE'. 

(3)  One  single  channel  at  one  space  pointr^  at  two  different  source  frequencies,  fg,  fg . 

(4)  etc. 

The  symbol  etc.  means  a further  listing  of  whatever  combination  of  two  events  in  real  space- 
time  or  transform  space-time  where  “time”  means  a choice  of  tE,  ts,  tSM,  tRM,  as  defined 
above. 

The  random  pressures  (or  one  random  pressure  and  its  conjugate)  are  used  to  form  two- 
point  statistics  in  space-time.  These  are  elements  of  a covariance  matrix  K,  if  ensembles  of 
the  random  pressures  are  available,  or  coherence  matrix  R,  if  only  one  realization  is  available. 
Consider  the  latter  since  they  are  most  common,  and  to  begin  with,  construct  cross  correla- 
tions of  two  channels.  A list  of  these  follows  in  which  m,  n denote  the  mth  channel  and  nth 
channel,  i,  j denote  two  receiver  space  points  , r^,  and  t,  t'  denotes  two  time  points, 


Correlation 

(1)  R[fV>(m),  f«)  '(n)  I t$(m),  t^'(n)] 

(2)  R[t<‘>(n),t«>'(n)  | t$(m),  t^>(n)] 

(3)  R[f<‘>(n),  f<J>(n)  | *$(m),  *$(n)] 

(4)  R[t<‘>(m),  tg'«)(n)  | <t><£(m),  ^'(n)) 

(5)  R[f<i>(m),  f«\n)  I k^>(m),  k<£>  '(n)] 

(6)  R[t<g^m),  t<j)'  (n)  | k<j>(m),  k^'(n)] 

(7)  R[k(;)(m),k(gj),(n)  | k^(m),  k$'  (n)] 

(8)  R[t^(m),  t£>'(m)  I (m),  *<|>'(m)| 
^)(m),  k<j)  (m)  | k<#  (m),  k^>  (m)] 

(9)  R[ty)(m),  t'W(m)  | k<;>  (*SM(m)),  k(8°' 

(4>SM(m))  | k^  (^>RM(m),  k$  (^rm^))) 

(10)  etc. 


Name  of  Correlation  Function 

(1)  Signal  frequency -environmental  time 

(2)  Signal  range  — environmental  time 

(3)  Signal  frequency-environmental 
frequency  (bifrequency) 

(4)  Range-doppler,  or  temporal 
scattering 

(5)  Frequency-angle  of  arrival 

(6)  Range-angle  of  arrival 

(7)  Spatial  (or  angle)  scattering 

(8)  Range-doppler-angle  of  arrival 
autocorrelation  or  ambiguity 
function. 

(9)  Range-doppler  autocorrelation  due 
to  motion  of  receiver. 


i 
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The  real  and  transform  coordinates  in  these  correlations  are: 


(1) 

(1) 

Range  -*■  frequency 

(2) 

~r 

*8  8 

(2) 

Location  vector  -*■  angle  of  arrival 
vector 

(3) 

(3) 

Location  vector  -►  angle  of  arrival 
vector 

(4) 

i 

tE  4*e 

(4) 

Environment  time  ->  environment 
doppler 

(5) 

*SM  ~ ^SM 

(5) 

Motion  time  -*•  motion  doppler 

(6) 

4RM  **  $RM 

(6) 

Motion  time  -*■  motion  doppler 

These  correlations,  and  many  more  of  the  same  type,  are  functions  of  space  and  time  and  thus 
correspond  to  nonstationary  processes.  Since  stationarity  is  important,  we  discuss  it  in  the 
next  section. 


Stationarity 


Each  element  in  the  matrix  E is  a random  process  in  space-time.  The  transmission  chan- 
nel and  the  arrived  signal  are  thus  also  random.  Because  in  general  the  ocean  is  nonstationary, 
and  because  the  construction  of  nonstationary  statistics  to  difficult  at  best,  it  has  been  cus- 
tomary to  consider  the  transmission  process  as  locally  stationary.  To  define  this  concept  it 
will  be  sufficient  to  take  one  of  the  random  variables,  say  the  environmental  matrix  E,  and 
focus  attention  on  one  of  its  coordinates,  say  the  environmental  time  tE . We  will  then  con- 
sider conditions  forjtationarity  in  this  selected  coordinate.  For  simplicity  we  take  only  the 
mth  component  of  E and  represent  it  as  a Fourier  transform, 


) = i*p  [i<t>EtE]  Z(dtE) 


in  which  Z(d<I>E)  is  a random  amplitude  with  mean  value  of  zero.  Now  the  autocorrelation 
of  E_  is, 


REm(tE.tE')  = fj  e'l*EtE  " 4>EtE  ^ ( Z(d4»E)Z*  (d4»'E)> 


In  order  for  Em  to  be  stationary  in  the  coordinate  tE  it  is  required  that  correlation  R depend 
only  on  difference  tE  - tE  and  not  on  tE,  tE'  separately.  Thus  one  must  set, 

< Z(d4>E)  Z*(d4»E’ ) > - - 4>e')  <fOJ>E  - 4*E')  d<i>E  d$p' 

in  which  is  the  frequency  power  spectrum  of  the  random  Em.  This  means  that  the 

random  spectral  amplitudes  are  uncorrelated  at  different  frequencies  4>E  if  E_  is  to  be  sta- 
tionary in  time  tE . 

In  a similar  way  stationary  of  Em(4>E)  in  frequency  requires  that 
( Z(dts)  Z*(dt/) ) = 5(ts  - ts')  F(dt,  dtg') 

so  that  R(4>E,  <I>E')  = R (4>E  - <j>E').  This  is  the  condition  of  uncorrelated  scattering  in  which 
the  scattering  of  that  portion  of  the  signal  at  time  ts  is  uncorrelated  to  the  scattering  at  time 
ts  . Stationary  of  Em  (ra,  rR ) in  space  is  defined  in  the  same  way.  For  example,  the  sta- 
tionary in  z-coordinate  requires  that  the  random  amplitudes  Z(dk7),  Z'(dk2)  be  uncorrelated, 

< Z(dkz)  Z*(dk'z)  > = 6(kz  - k'z)  F(dkz,  dk'z) 


This  means  that  different  ray  paths  arriving  at  a point  are  uncorrelated  in  the  z-component 
of  the  angle  of  arrival. 

Propagation  channels  which  obey  stationary  in  the  limited  sense  describe  here,  and 
which  have  constant  mean  values  in  the  transmission  of  a field  are  wide-sense  stationary. 

Few  channels  of  ocean  propagation  which  are  long  enough  in  space-time  are  stationary.  How- 
ever, over  a “short  enough  interval”  in  space-time  the  increment  in  the  random  ocean  process 
may  be  stationary.  For  example,  let  Em  be  a random  process  with  nonstationary  mean  in 
the  coordinate  time  tE.  Choose  a “small-enough”  increment  in  time,  call  it  TF,  and  form 
the  increment 


AEm  = Em  (4E  +TE)-Em  (tE) 

If  the  low  frequency  swings  in  (tE)  occur  so  slowly  that  the  increment  A Em  is  unaf- 
fected by  them  then  the  nonstationary  function  Em  (tE)  is  defined  to  be  stationary  in  its 
first  increments.  The  second  order  statistics  of  A Em  are  then  described  by  the  structure 
function,  which  is  the  ensemble  average  of  the  square  of  its  zero-mean  part: 

DEm  (TE)  = ((AEm-  <AEm>)2> 

The  properties  of  D are  self  evident:  (1)  D(0)  = 0 (2)  D(«>)  js  finite  (3)  if  E_(tE)  is  a 
random  process  with  zero-mean  then 

°Em  (TE)  = 2 [REm(0)  - REm  (Te)] 
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From  the  definition  of  the  concept  of  structure  function  it  is  seen  that  it  can  be  applied  to 
each  coordinate  of  space,  or  time,  or  both,  and  to  each  random  process  of  ocean  propagation, 
particularly  the  acoustic  processes  itself.  Thus  structure  functions  of  received  pressure  are  in 
use  wherever  nonstationity  is  tentatively  accounted  for. 


Antenna  Function  and  Total  Received  Pressure 

The  pressure  received  by  the  ith  hydrophone  in  an  L element  hydrophone  array  is  cus- 
tomarily modified  by  an  antenna  steering  function  A which  shades  the  received  amplitude 
and  introduces  a time  delay,  rd . In  actual  practice,  A is  a random  process  of  space  and  time, 
with  the  coordinate  of  randomness  in  time  being  tA . Summation  over  all  hydrophones  then 
gives  the  total  received  pressure, 

L 

P(fR  (0,*),tE,ti,tA)=  Z]  A(^(0,^)|r$,rd,tA)p(r^),tE,tg) 


in  which  rR  (6  ,<p)  is  the  coordinate  of  spatial  distribution  of  the  received  pressure  and 
p(rl'),  tE,  tg)  is  a sum  of  multipaths.  A plot  of  |p|  vs  0,  <fi  gives  the  beam  pattern  of  the  re- 
ceived pressure,  which  because  of  both  p(')  and  A,  is  a random  variable.  The  random  aspects 
of  A can  be  studied  in  the  same  way  as  for  the  environmental  matrix  l£,  namely  by  formation 
of  a covariance  matrix  in  the  time  coordinate  tA , or  in  the  space  coordinate^) , or  in  both. 
Alternatively  one  can  form  a coherence  matrix  R of  correlations  of  A in  the  same  coordinates. 

The  coherence  matrix  in  correlations  of  the  total  pressure  allows  one  to  obtain  the  array 
gain  AG  as  a function  of  steering  rd  and  random  perturbation  rA  in  S/N  power  in  the  manner 
suggested  earlier:  ' 


L A(i,  tA ) R(*)  (i,  j,  rd ) A(j,  tA ) 
i,  j [R<*)  (i,  o)  R(»)  (j,  o)] l'2 


[R8  (i.o)  R(»)  (j,o)]i/2 


L A(i,  tA)  Rn  (i,  j,  rd)  A(j,tA) 

AG(Td,tA)*<)  

i,j  [R(N)(i,o)R(N)(j,o)]i/2 


[RN  (i,  o)  R(j,o)p/2 


R(»)  (o) 

Rn  (o) 

in  an  obvious  shortening  of  notation.  Once  again,  since  the  A’s  are  random  we  noted  that 
the  array  gain  is  a random  function  of  space-time. 
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Coherence  and  Fluctuations 

The  temporal  cross  correlation  of 
is, 


R(x1,x2>r)  = 


If  the  signal  at  X2  simultaneously  undergoes  a shift  in  phase,  6Jdt,  the  cross  correlation  fea- 
tures an  additional  parameter  cud , 
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R(xltx2lr,(jd)  = J p(x1 , tg)  p*(jT2,  tg  + r)  e~>LJ  dt*  dt, 

-T/2 

The  normalized  absolute  value  squared  of  this  quantity  is  the  ambiguity  function. 


two  signals,  p(xa , t)  and  p(x2,  t)  over  time  duration  T 
T/2 

/ P(xi.tg)p(x2,tg +T)dts 


-T/2 


If  the  signal  at  ^ is  the  transmitted  signal,  and  that  at  is  the  same  signal  delayed  t units 
and  phase  shifted  cudt  units,  this  function  measures  the  mean  square  departure  of  the  trans- 
mitted signal  from  its  original  waveform.  A three-dimensional  plot  of  lx(r,  tod ) I2  vs  r and 
cud  is  used  in  sonar  detection  and  localization  because  it  exhibits  a peak  at  the  supposed  loca- 
tion where  there  is  a target.  However,  the  peak  is  usually  broad,  and  is  generally  accompanied 
by  other  peaks,  indicating  the  uncertainty  of  the  presence  of  a target  in  the  midst  of  false  alarms. 


The  search  for  target  by  integration  over  the  time  domain  is  advantageously  replaced 
in  many  instances  by  a search  in  the  frequency  domain.  The  method  of  constructing  the 
appropriate  mathematical  surface  in  which  peaks  are  to  be  found  is  briefly  as  follows.  Let  a 
record  of  pressure  p(t)  be  represented  by  N values  p(mAt),  n = 0, 1,  2 . . . N.  The  discrete 
Fourier  transform  of  these  N points  into  N points  in  frequency  yields  the  sequence, 

j-  = |g(o)  + G(  Aco)  + G(2A  oj)  + ....  + G((N  - 1)  A o>)  J 

in  which 


N-  1 

G(mAw)  = 21]  g(mAt)  expr^-mnAwAt  1 

m = 0 L N J 

Each  G is  thus  a sum  of  N exponentials  or  phasors  for  each  value  of  frequency  m Aco.  In  this 
way  signals  at  space  points  Xj , X2  can  be  represented  in  the  frequency  plane.  Similarly  the 
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cross  correlation  of  signals  in  the  time  domain  has  its  analog  in  the  crosspower  spectrum  in  the 
frequency  domain.  In  particular,  the  ambiguity  function  at  a single  frequency  (mAco)  becomes, 


l£  ** 


S^x2  (mAoj.r,  a>d)  = 


GX1  (m  Acj)  G*^2  (mAu,  t) 


N i N . 

g|G^(mAw)|2  2Z|  G(mAcu,  t)|2 


A plot  of  vs  r and  tod  is  a topological  surface  of  the  cross-power  spectrum  between  sig- 

nals at  Xj  and  X2  calculated  on  the  basis  of  a single  frequency.  The  highest  peaks  of  all  peaks 
of  this  surface  identifies  possible  targets,  but  with  ambiguity  caused  by  the  br"  ^ness  of  the 
peaks  and  by  the  presence  of  false  alarms  (smaller  peaks). 

— ► — ► “► 

The  locations  X2 , X2  can  refer  to  several  points  of  measurement  of  signals:  X2  can  be  the 
source  location  and  X2  the  point  of  reception;  or,  Xj , X2  can  be  the  location  of  two  hydro- 
phones in  an  array,  or  two  arrays  in  a multi-array  receiving  system.  In  all  cases  plots  of  ,x2 
vs  r and  tod  are  labeled  coherence  surfaces  of  signals  at  XltX2,  t is  called  the  range,  or  the 
time  difference  of  arrival,  and  cud  is  called  the  differential  Doppler. 

The  factor  exp  (-jlcud)  in  Sxlx2  is  significant.  Since  £ (or  discrete  time)  and  cod  are  a 
Fourier  pair,  it  is  seen  that  increasing  the  limit  of  £ from  N to  some  higher  integer  simultan- 
eously decreases  the  width  of  cod,  that  is,  it  increases  the  resolution  in  Doppler.  One  says 
that  increase  in  the  integration  time(meaning  the  sum  over  £)increases  the  resolution  in  u>d . 
Similarly  in  the  form, 


G(mAco,  r)  = Y1  g(ni(A  t + r))  exp  ( — — mnAcu(At  + r)) 
t—t  N 


if  one  increases  the  limiting  value  of  m which  is  equivalent  to  increasing  the  bandwidth  of  the 
processing,  the  resolution  in  differential  time  of  arrival  t is  also  increased. 

The  mathematical  content  of  coherence  surfaces  Ix^i  ^2 I2  and  ^2  so  far  have  not 
accounted  for  ocean  fluctuations.  In  actual  measurements  in  the  ocean  these  surfaces  have 
been  found  to  fluctuate  as  functions  of  environmental  time  tE  and  multichannel  interferences 
ZTm  caused  by  relative  motion  of  source  and  target,  whether  opening  or  closing.  The  study 
of  coherence  surfaces  may  thus  be  distinct  from  the  study  of  fluctuations.  The  latter  is  essen- 
tially a study  in  statistics,  the  former  is  essentially  deterministic. 

Modeling  Fluctuations  in  Signal  Excess  in  Sonar  Detection  Models 

Each  term  SEj  in  the  equation  of  signal  excess  represents  a random  variable  in  space-time. 
For  a fixed  point  in  space  the  stochastic  process  SEj(t)  can  be  modeled  in  several  ways  in 
accordance  with  the  requirements  of  the  underlying  physical  process.  In  most  applications 
the  mean  SEj(t)  although  the  function  of  time,  is  considered  deterministic.  Attention  is  then 
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focused  on  the  fluctuating  part  X(t).  Given  the  great  diversity  of  causes  of  ocean  randomness 
it  would  seem  that  the  construction  of  a single  model  of  X(t)  is  not  feasible.  However,  if  the 
underlying  physical  process  leads  to  a normal  probability  distribution  of  X(t)  and  an  auto- 
correlation which  is  exponential  in  time,  and  if  a plot  of  X(t)  vs.  t looks  like  a random  walk 
with  discrete  steps  of  relatively  small  size,  one  could  usefully  adopt  a random  walk  model  for 
X(t). 

Fig.  III.2  shows  two  realizations  of  a random  function  Y(t)  pictured  as  a random  walk: 
the  abscissa  is  Y(t),  quantized  into  distinct  excursion  steps,  or  states,  Y2,  Yj , Y0,  Yj , 

Y2 while  the  ordinate  is  quantized  in  time,  0,  rE,  2te> Realization  1 

consists  of  the  walk  Yx,  Y0,  Y.j , Y0,  Y.lt  Y.2,  Y.j  while  Realization  2 describes  another 
walk  Y_! , Y.2,  Y.! , Y0,  Y.j , Y0,  Yj  . The  probability  of  the  occurrence  of  Realization  1 
is  modeled  as  a Markov  chain, 

p J 1 1 = al  Pl.o  Po,-i  P-i,0  Po.-i  P-1,-2  P-2,-1 

in  which  a1  is  the  initial  probability  that  Y(t)  is  in  state  Yj ; px  0 is  the  conditional  (or  tran- 
sition ) probability  of  going  from  state  Y1  to  state  Y0;  p0  j is  the  transition  or  proba- 
bility of  going  from  state  Y0  to  state  Ylt  etc.  The  probability  of  occurrence  of  Realization  2 
(P  {2}),  is  similarly  constructed.  In  numerical  work  one  must  assign  probability  distributions  to 
aj,  andpjj,  depending  on  the  physics.  Since  these  choices  are  critical  a variety  of  them  have 
been  proposed  to  fit  measurements.  Suppose  the  chain  of  a + 1 states  exhibits  a tendency  to 
cluster  about  a central  state  Yc,  and  that  the  movement  from  a high  number  state  toward  the 
central  state  is  always  more  likely  to  happen  than  away  from  it.  The  matrix  of  transition 
probabilities  is  then  modeled  as, 
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A Markov  chain  of  Y(t)  that  features  these  transition  probabilities  is  an  Ehrenfest  random  walk. 
The  entries  mean  the  following:  if  Y(t)  is  initially  in  state  Y0  it  has  zero  probability  of  remain- 
ing there,  p00  = 0.  In  contrast  it  moves  to  state  Y1  with  certainty,  p01  = 1.  If  Y(t)  is  in  state 
Yx  it  has  zero  probability  of  remaining  there,  px  1 = 0.  It  either  moves  to  state  Y0  with  proba- 
bility a-1,  or  to  state  Y2  with  probability  1-a-l.  It  has  zero  probability  of  moving  to  any  other 
state  than  the  one  immediately  to  the  left  or  to  the  right.  Since  a is  usually  large  the  movement 
favors  Y2.  If  Y(t)  is  in  state  Y2  it  again  has  zero  probability  of  remaining  there.  Instead  it 
moves  to  state  Yx  with  probability  2a- 1,  and  to  Y3  with  probability  1-28-1.  Once  again  Y3  is 
more  likely.  Now  suppose  Y(t)  is  in  state  Yj,  j > a/2.  Then  the  movements  Yj.j  has  probability 
j/a,  which  is  greater  than  1/2,  and  to  Yj+1  with  probability  1-j/a  which  is  less  than  1/2.  The 
movement  back  to  a lower  state  is  therefore  most  likely.  Thus  Y(t)  in  this  model  is  a diffusion 
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with  a central  force.  The  random  walk  pictured  by  Y(t)  clusters  around  a mean  value  a/4, 
with  a standard  deviation  of-y/  a/4.  In  order  to  model  the  fluctuation  X(t)  it  must  have  the 
mean  value  removed.  When  this  is  done  one  can  form  the  ratios, 

SE-SE  _ Y-a/4 

°SE  OSE  V a/4 

in  which  OgE  is  the  standard  deviation  of  the  signal  excess  (an  experimentally  determined 
quantity).  Now  the  probability  distribution  of  Y(mtE),  m = 0, 1,  2 ...  is  known  to  be  bino- 
mial with  mean  a/4,  and  variance  a/4  (see  Feller).  By  choosing  the  number  a+1  of  states  of  Y 
Y one  can  construct  a random  number  generator  from  which  realizations  of  sequences 
Y0,  Yj  . . . Ya  can  be  obtained,  modeled  as  Ehrenfest  random  walks,  with  mean  a/4  and 
variance  a/4.  From  these  sequences  one  can  obtain  realizations  X0,  Xj , . . . Xa  by  multi- 
plication with  Ogg.  Thus  any  sequence  {X;}  can  be  generated  as  a Monte  Carlo  trial  begin- 
ning with  the  state  Yj,  which  itself  is  selected  according  to  some  probability  distribution 
(usually  normal). 


Time  Interval  of  Stepping,  te 

A random  walk  in  time  requires  specification  of  a stepping  interval  rE.  Now  if  X(t)  is 
taken  as  a general  Gauss-Markov  process  its  normalized  correlation-function  is  exponential, 

At 

Px.(t 2 - ) - e rR  A t = 1 12  - tj  | 


in  which  rR  is  the  “relaxation  time”,  or  that  time  interval  t2  - ta  when  px  = e-1 . In  contrast 
the  autocorrelation  function  for  the  Ehrenfest  random  walk  is 


in  which  S is  the  number  of  steps  of  Y(t)  to  be  used  in  forming  pE.  If  Y(t)  is  to  model  X(t) 
one  requires. 


= e 


For  a single  step  one  sets  S = 1.  This  then  defines  the  interval  rE  = At. 

a - 2 . 

= e " teItr 


te  =-tr  In 


*W.  Feller,  “An  Introduction  to  Probability  Theory  and  Its  Applications,”,  Vol.  I. 
Wiley,  New  York  (1943). 
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Thus  the  stepping  interval  in  time  (=te)  is  determinable  as  soon  as  the  relaxation  time  rR  of  the 
underlying  Gauss-Markov  physical  process  is  known,  and  a number  a + 1 of  states  is  specified. 
Random  realizations  of  Y(mrE)  are  then  supplied  by  the  random  number  generator  in  the  form 
of  steps  of  the  Ehrenfest  random  walk,  from  which  various  X(mrE)  are  determined. 

Environment , Target  and  Sensor  Fluctuations 

In  more  refined  models  of  detections  of  sonar  the  total  fluctuation  X(t)  is  taken  to  be  a sum 
of  fluctuations  of  terms  in  the  equation  for  signal  excess  SE  = SE  + X(t).  Typically, 

X(t)=XE(t)+Xt(t)+Xs  (t) 

in  which  environment  (or  transmission  loss),  target,  and  sensor  (hydrophone  or  array)  are  repre- 
sented. Each  Xs  can  be  modeled  as  a statistically  independent  random  walk  with  measured  stan- 
dard deviation  ax/j,  measured  relaxation  7-r^j,  and  selected  stepping  interval  r En\  as  determined 
from  and  the  selected  number  of  states  (a+l)^).  The  variance  of  the  total  X(t)  is  then, 

°SE2  = °E2  + °R2  + °S2 

The  variations  of  X;  with  time  t make  it  sensitive  to  frequency.  Thus  the  assignment  of  0x(j) 
and  will  vary  with  frequency. 

SYSTEM  SIMULATION  (SS) 

Underwater  surveillance  systems  are  visualized  in  Fig.  III. 2.  The  key  features  are:  (1) 
number  and  location  of  signal  sources  Sn,  m = 1,  2...;  (2)  number  and  location  of  noise 
sources  Nm,  m = 1,2...;  (3)  spectrum  of  the  signal  source,  or  noise  source;  (4)  number  of 
sensors,  Ah('),  i = 1,2...;  A v<*\  i = 1,2...;  (5)  number  and  character  of  transmission  channels. 

This  visualization  in  its  entirety  is  simulated  (not  duplicated  in  a physical  sense)  in  its 
dynamic  operations  (detection,  localization,  tracking)  by  the  computer  program  SS.  The 
simulation  is  on  a very  large  scale:  it  handles  a large  number  of  separate  sources  of  signal  and 
noise,  several  frequency  lines  per  received  target,  and  a large  number  of  sensors.  Simulation 
begins  with  a sonar  equation  for  each  source-channel-receiver  to  determine  signal-to -noise 
ratio.  Since  each  term  in  the  equation  is  a random  variable  there  is  a need  to  specify  its  proba- 
bility distribution  and  statistics.  The  mean  value  of  each  term  is  obtained  by  interpolation 
from  a table  of  mean  values  previously  stored  in  the  memory  of  the  computer  from  accumu- 
lated field  data.  Its  value  is  stationary  with  time.  The  fluctuating  component  of  each  term  is 
modeled  as  a Gauss-Markov  process  (or  chain),  meaning  that  in  successive  realizations  of  time 
record  of  length  O to  T,  the  probability  between  immediately  adjoining  points  is  Markov. 

Since  the  sonar  equation  is  in  logarithmic  units  the  probability  distribution  of  the  Gauss- 
Markov  chain  is  taken  as  lognormal  with  a variance  (or  standard  deviation)  to  be  supplied  to 
SS  as  drawn  from  field  data.  The  second  order  characterization  of  the  Gauss-Markov  simula- 
tion is  the  temporal  autocorrelation  of  a single  realization.  Since  this  correlation  is  assumed 
exponential  with  a characteristic  relaxation  time  rR  field  data  must  again  be  used  to  furnish 
this  to  SS  as  an  input.  Other  acoustic  inputs  to  SS  are  listed  in  the  proceedings  paper  by 
Fretwell. 
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When  all  the  acoustic  input  to  PSEUDO  is  completed  the  computer  program  selects  from  a 
random  number  generator  based  on  Gauss-Markov  statistics  and  from  the  input  field  data  a set  of 
random  numbers  to  represent  realizations  of  fluctuation  processes  in  target  strength,  transmis- 
sion loss,  ambient  noise,  beam  noise  and  signal  gain.  Mean  values  of  these  same  quantities  are 
obtained  from  an  interior  stored  table  by  simple  look-up  guided  by  the  acoustic  input.  PSEUDO 
then  completes  one  run,  delivering  a single  realization  of  received  acoustic  pressure  versus  time. 
The  program  is  then  returned  to  initial  conditions,  a new  set  of  random  number  is  drawn,  and  a 
second  realization  is  completed.  The  procedure  is  repeated  many  times  to  build  up  a sequence 
of  realizations  of  sufficient  number  to  allow  ensemble  averages  to  be  made  with  confidence,  and 
hence  to  allow  a statistical  description  of  the  SNR  of  the  received  signal  to  be  made.  From  this 
system  performance  can  be  estimated:  (1)  detection  is  declared  when  the  SNR  rises  above  a 
threshold  (statistically  selected),  (2)  classification  is  determined  by  signal  to  noise  ratio 
criteria,  (3)  localization  and  tracking  are  accomplished  by  use  of  the  SNR’s  in  separate 
algorithms  programmed  into  PSEUDO.  PSEUDO  is  used  primarily  as  a comparator  of 
rival  or  competing  underwater  surveillance  systems  declaring  some  to  be  better  than 
others  within  the  assumptions  built  into  its  algorithms.  Since  it  does  not  duplicate  physi- 
cal processes  in  the  ocean  its  validity  must  be  sought  for  in  the  context  of  Monte  Carlo 
trial  on  simulated  random  processes  rather  than  in  realization  by  realization  agreement 
with  physical  measurement. 


Solving  for  Channel  Transmission  Operator  Tm , Tm-1 

The  ocean  is  stratified  in  the  depth  coordinate  z,  and  the  properties  of  stratification  vary 
with  the  range  coordinate  X.  At  a point  ?in  the  medium  and  at  a particular  frequency  co  the 
pressure  field  can  be  represented  by  an  angular  (spatial)  spectrum  of  plane  waves  in  the  trans- 
verse wavenumber  kt  and  transverse  coordinate  p = jy  + kz 

p(? <*>)=![  x)  e‘kt  * P d^t 

— oo 

— k * * 

kt  =jky  +kkz 


The  plane  wave  amplitudes  F are  generally  not  expressible  in  closed  form.  However,  for  slow 
medium  variations  F may  be  approximated  by  the  WKB  representation, 

F(kt,co;x)  = A(kt,  cu;x)  exp-^ij  kx  (kt,  co;  £)  d£  j 
Thus  the  angular  spectrum  representation  of  the  pressure  field  becomes, 

P(r,t)  = e~iwt  dw  j ff  A(kt,cj;x)exp|  i kt'~p+J'  kx  (kt,  «;  £)  d{  jjdkt 

Here  the  sjjunje  of  the  field  is  not  explicitly  represented.  If  a point  source  in  space  is 
present,  Q = 5(p  - p')  6(x  - x'),  the  amplitudes  F are  the  longitudinal  Green’s  functions  G, 
which  in  the  WKB  approximation,  become 

F = G(kt,cu;x,x')  = A(kt,  to;x,x')exp|i  / / [ kx  (kt,  u;  $)  - kx  (kt,  co;  $')  jj"d  f d £' 
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The  four-dimensional  space  time  Green’s  function  is  then  constructible  in_the  form, 
, , CCr  ->  , dk[  dcu 

G(r,  r;  t,  t')  =JJJ  0(1^,  «;x,  x')  eikt(o  ’ * > e’iw  <*"*  > ~^3 


When  the  source  is  distributed  in  space/time,  Q(rg,  tg),  the  resultant  pressure  field  at  a specific 
frequency  w is  r 

p(rR , u)  = / G(rR ,rs  |u)  Q(r  | w)  dr, 

J vol 
or 

V(rn,cj)=f^{[  A(kt,w;xR,xg)expji  jf  Rk,  (kt,  w;  £)  dt|e,kt  * <PR - P»> 

X Q(r^  | «)  dkt  dr^ 

where 

^ = i x,  + J y«  + kzg  = i xg  + p. 


The  relation  of  these  forms  to  the  operator  Tm  is  seen  in  the  definition. 

Pm  ~ TmQm 

Thus,  in  the  WKB  approximation,  the  channel  operator  Tm  becomes, 

Tm  -*  [ [[  A(kt,o>;  x,  xg)  exp  i{  f kx  (kt,w;  £)  d£ 

vol  '-X. 


Xeikt'(PR  -P«) 


| dkt  drg 


A principal  feature  of  this  operator  is  integration  over  the  range  (X)  component  of  the  propa- 
gation vector  to  give  the  WKB  approximation  to  phase. 

This  form  of  Tm  is  deterministic.  However,  in  ocean  propagation  the  pressure  at7R  is  a 
random  variable.  To  represent  this  condition  the  amplitudes  A are  made  random  variables. 
Thus  Tm  for  the  general  transmission  channel  is  a random  variable. 

The  Basic  Tracking  Problem  in  System  Performance  Modeling 

A submarine  radiating  specific  sounds  follows  a prescribed  underwater  track  in  a 
selected  area  of  the  ocean  (Fig.  III.3).  Several  listening  receivers  (A  through  F)  are  re- 
quired to  track  the  path  of  the  submarine.  The  received  signals  and  noise  at  these  arrays 
form  time  series  (Fig.  III.4).  These  exhibit  random  behavior  about  some  mean.  At 
several  times  in  each  record  there  are  sudden  drop-outs,  where  occurrences  are  random. 

At  other  times  the  signals  fade  out  completely.  Let  a threshold  TH  be  set  for  each  time 
record,  and  assume  that  when  SNR  is  above  TH  the  submarine  is  detected.  Because  of  the 
random  nature  of  these  records  several  probability  statements  on  detection  are  required  to 
be  formulated.  Two  examples  are: 

27 


Fig.  III.  4 

Received  SNR’s  of  Basic  Tracking 
Problem 


SNRn 
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(1)  What  is  the  probability  distribution  that  a specific  array  will  track  the  target  at  least 
X minutes?  (2)  What  is  the  joint  probability  that  two  (three,  four,  etc.)  arrays  will  track 
the  same  target  at  least  X minutes.  The  answers  to  these  questions  and  others  in  existing 
systems  are  of  course  obtained  by  trial.  However,  proposed  new  systems  and  proposed 
improvements  on  old  systems  must  be  critically  evaluated  between  competing  candidates, 
making  system  performance  prediction  mandatory.  The  construction  of  prediction  algo- 
rithms becomes  thus  a basic  problem. 

A first  approach  to  system  performance  modeling  is  to  duplicate  the  underlying 
physics  directly  by  mathematical  formulas  and  numerical  analysis.  In  this  procedure  the 
acoustic  pressure  arriving  at  arrays  A through  F by  multipath  channels  through  a randomiz- 
ing ocean  is  obtained  by  physical  analysis.  This  yields  a spatial  distribution  of  mean  value, 
which  is  locally  a function  of  time,  and  spatial  and  temporal  auto-and  cross-correlations 
(or  mutual  coherence)  of  the  fluctuating  component  of  the  field.  Correlation  lengths 
correlation  times,  power  spectra,  etc.,  are  predicted  from  a knowledge  of  the  environment. 
This  approach  is  shown  in  Fig.  III.  5 
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Figure  III.  5 


Using  this  model  one  computes  acoustic  intensity  and  thus  one  obtains  transmission  loss. 
Other  components  of  the  sonar  equation  are  similarly  treated  e.g.,  analysis  of  underlying 
physics  is  made  for  ambient  noise,  and  a model  constructed.  Also,  by  means  of  statistical 
antenna  theory  one  constructs  a mathematical  model  of  the  array  gain,  beamwidth,  etc. 
which  includes  spatial  and  temporal  fluctuations. 

All  these  analytic  models  of  sonar  equation  components  contribute  to  the  SNR  at  each 
array.  Since  each  SNR  is  a random  process  it  will  have  a predicted  probability  distribution, 
from  which  all  the  statistical  moments  can  be  drawn.  Thus  the  statistics  of  SNR  at  each 
array  can  be  predicted  for  each  position  of  the  submarine  along  the  track.  Furthermore,  the 
models  also  contain  algorithms  for  predicting  the  spatial  and  temporal  civ  a correlation  of 
total  SNR  between  arrays.  This  gives  estimates  of  SNR  interarray  mutual  coherence,  SNR 
correlation  time  and  correlation  length,  which  in  turn  can  be  used  to  calculate  the  real 
physical  tracking  statistics  of  real  systems,  and  hence  can  be  used  for  system  performance 
prediction. 


A second  approach  to  system  performance  modeling  is  to  devise  a simulation  algorithm 
which  models  the  tracking  problem  as  a random  process,  chosen  for  convenience  as  plausible, 
but  not  specifically  addressed  to  the  underlying  physics.  For  each  track  position  the  acoustic 
intensity  (or  transmission  loss)  from  submarine  to  arrays  A through  F is  simulated  as  the  sum 
of  a deterministic  mean  value  and  a random  (fluctuating)  value.  The  mean  value  at  each  array 
is  obtained  from  entries  in  a previously  prepared  table  which  are  based  on  past  experience  and 
past  modeling,  in  similar  ocean  areas  at  similar  ranges  between  source  and  receiver.  The  fluctu- 
ating part  of  received  acoustic  intensity  at  each  array  is  modeled  as  a plausible  temporal  ran- 
dom process.  Currently,  the  popular  choice  is  a Gauss-Markov  process  equivalent  to  a random 
walk,  implemented  in  the  form  of  a random  number  generator.  Each  number  delivered  by  this 
generator  is  a step  of  random  amplitude  to  be  added  to  the  mean  value  of  transmission  loss 
already  obtained  for  the  specific  track  position  and  specific  listening  array,  the  total  being  a 
realization  of  random  transmission  loss.  At  the  same  track  position,  other  components  of  the 
sonar  equation  are  similarly  realized,  and  hence  a single  realization  of  signal  excess  is  obtained 
for  each  array.  Upon  moving  to  the  next  position,  and  eventually  to  all  positions,  a final 
sequence  of  simulated  time  series  for  signal  excess  at  all  the  arrays  is  obtained,  Fig.  III.6.  The 
result  may  be  called  Realization  Set  of  SE  1.  After  completion  of  the  track,  the  simulator 
returns  again  to  position  1 and  repeats  the  entire  tracking  of  the  course  to  yield  Realization 
Set  of  SE  2.  Repetition  then  commences  a third  time.  After  a large  number  of  realizations 
are  obtained  the  simulator  algorithm  processes  the  ensemble  of  them  to  furnish  all  the  tracking 
statistics  needed.  However,  it  is  understood  that  the  statistics  so  obtained  form  a plausible 
facsimile  of  the  real  tracking  process.  They  cannot  be  used  to  understand,  or  predict,  acoustic 
events  .in  the  ocean,  or  physical  tracking  of  real  submarines. 
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Figure  III.  6 
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PERFORMANCE  CRITERIA  AND  STOCHASTIC  SIGNAL  EXCESS 

The  performance  of  underseas  surveillance  systems  in  their  task  of  detecting,  classi- 
fying and  tracking  hostile  submarines  is  modeled  on  several  criteria.  Two  of  these  cur- 
rently in  use  are  cumulative  probability  of  detection  and  mean  time  to  hold  contact.  A 
simple  model  of  each  is  to  write  it  as  the  sum  of  a deterministic  part  and  a random  part, 
both  being  functions  of  time.  For  example,  a parameter  of  these  criteria,  signal  excess 
S(t),  is  written, 

(1)  S(t)  = M(t)  + £(t)  (2)  M = SL  - TL  - (NL-DI  + DT) 

in  which  M(t)  is  deterministic,  based  on  prior  knowledge  of  source  level  SL,  transmission 
loss  TL,  noise  level  NL,  directivity  index  DI,  and  detection  threshold  DT;  and  £(t)  is  ran- 
dom, representing  the  effects  of  random  changes  in  ocean  environment,  target  character- 
istics, sonar  equipment,  operator  efficiency,  etc.  The  calculation  of  M(t)  can  be  done  in 
the  conventional  manner.  The  stochastic  term  is  more  elusive  in  that  it  requires  a prob- 
abilistic statement  tailored  to  the  mode  of  operation  (whether  active  or  passive).  For 
example,  during  passive  detection  one  requires  an  estimate  of  the  probability  that  signal 
excess  is  greater  than  zero  during  some  time  interval: 

Pr{S(t)  > 0,  0 < t < T}  . 

Likewise,  during  tracking  by  a number  m of  sensors  each  of  which  can  hold  the  target 
only  part  of  the  time,  the  operator  needs  an  estimate  of  the  probability  that  some  j 
sensors  of  the  total  m,  during  some  interval  T,  will  exhibit  a signal  excess, 

Pr{Sj(t)  >0,  1 < j < m,  0 < t < T}  . 

Naturally,  a prediction  model  of  these  probabilities  would  be  of  great  importance.  An 
oft-quoted  one  (now  mainly  historical)  is  called  the  X,  o model.  In  it  the  stochastic  signal 
excess  S(t)  is  taken  to  be  normally  distributed  with  zero  mean  and  standard  deviation  o. 
The  probability  that  a signal  excess  will  occur  in  an  interval  x is  written  4*(x).  Since  not 
all  intervals  will  generate  favorable  detection  levels,  their  occurrence  is  assumed  to  be  a 
Poisson,  or  shot,  process.  The  model  for  probability  of  detection  is  then  written  as 

Pd(t)  = 1 - [1  - $(a)]e'^(*)‘ 
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in  which  X is  the  rate  of  favorable  arrivals.  It  predicts  that  when  signal  excess  is  negative  the 
probability  of  detection  is  not  zero,  even  in  short  observation  times.  Indeed  the  model  pre- 
dicts that  for  signal  excess  of  -a  dB  the  probability  at  the  end  of  6.2  relaxation  times 
T(T  = 1/X)  is  nearly  0.7,  a large  value. 

In  the  years  1969-1972  an  effort  was  made  by  NSRDC  in  conjunction  with  Com- 
mander Submarine  Development  Two  to  use  data  of  at-sea-exercises  to  estimate  X,  a.  The 
results  are  somewhat  unsatisfactory,  the  fundamental  difficulty  being  the  lack  of  con- 
sistency from  one  exercise  to  the  next.*  Thus,  while  the  problem  of  modeling  cumulative 
detection  probability  for  stochastic  processes  has  been  current  for  a decade,  it  still  con- 
tinues. Part  of  the  continuing  problem  is  that  in  all  models  the  stochastic  variable  of  signal 
excess  cannot  be  calculated  from  a simple  enough  algorithm  in  such  a way  to  be  both  accurate 
and  practical  at  the  same  time.  It  is  a focal  issue  of  this  workshop  to  assess  current  ability  to 
provide  such  an  algorithm,  tailored  to  the  needs  of  the  system  performance  modeler. 

COMPONENTS  OF  THE  PROBLEM  OF  DESCRIBING  £(t) 

The  stochastic  part  of  signal  excess  is  an  effect  of  ocean  environment,  and  source- 
receiver  characteristics.  The  components  are  pictured  as  shown  in  Figs.  1(a)  through  1(c). 

Figure  1(a)  shows  a background,  or  equilibrium  ocean,  having  surface  and  bottom 
with  a distribution  of  sound  speed  profiles.  Figure  1(b)  shows  a perturbed  ocean  in  which 
the  profiles  undergo  random  variations  in  space  and  time.  Figure  1(c)  shows  volume  dis- 
tributed sources  and  receivers  in  horizontal  and  vertical  motion.  The  general  problem  of 
describing  the  random  part  of  signal  part  of  signal  excess  £(t)  consists  in  assembling  all 
three  pieces.  Figs.  1(a),  1(b),  and  1(c),  constructing  a model  of  the  propagation  of  an 
acoustic  signal  in  the  assembly,  and  preoicting  the  statistical  characteristics  of  the  received 
signal.  This  signal  is  then  inserted  in  the  sonar  equation,  a threshold  set,  a signal  excess 
calculated,  and  from  it  the  cumulative  probability  of  detection  is  predicted. 

To  keep  within  reasonable  bounds  we  restrict  attention  hereafter  to  the  acoustic 
problem  alone,  which  concentrates  on  finding  the  statistical  characteristics  of  the  received 
signal  as  a function  of  the  components  shown  in  Figs.  1(a),  1(b),  and  1(c)  and  of  the  fre- 
quency content  of  the  acoustic  signal.  We  now  address  the  first  issue. 
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THE  BACKGROUND,  OR  EQUILIBRIUM,  OCEAN 

To  begin  the  problem  of  modeling  (that  is,  understanding)  fluctuations  in  the  ocean,  we 
must  inquire,  which  ocean?  The  acoustician’s  ocean  is  contained  in  dynamic  relations  be- 
tween seven  variables;  3 components  of  convection  velocity,  u,  v,  w;  mass  density  p,  total 
pressure  p,  salinity  S,  and  temperature  T.  Seven  equations  are  need  to  show  the  interrela- 
tions of  these  variables.  The  first  is  an  equation  of  state  relating  p,  p,  S,  T.  The  second  is 
an  equation  of  continuity  which  accounts  for  temporal  changes  in  mass  density.  The  third 
is  an  equation  of  diffusion  of  salinity  S.  The  fourth,  fifth  and  sixth  are  three  dynamic 
equations  of  motion  in  u,  v,  w.  Finally,  the  seventh  is  the  equation  of  heat  flux. 


♦During  the  same  time  more  elaborate  models  were  constructed  by  Navy  laboratories  and  Navy  contractors, 
the  chief  being  APSURV,  PSEUDO,  APAIR,  etc.  These  will  be  discussed  by  other  speakers.  They  contain 
to  this  day  several  difficulties  which  make  their  predictions  uncertain,  and  which  this  workshop  will  seek 
to  resolve. 
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Naturally  we  would  prefer  constructing  models  on  the  most  realistic  description. 
However,  the  ocean  description  in  this  case  leads  to  very  cumbersome,  perhaps  intractable, 
mathematics.  But  there  are  minimum  statements.  An  irreducible  minimum  choice  is  that 
the  ocean  is  channeled,  inhomogeneous  and  anisotropic.  Secondary  but  important 
statements  Eire  that  the  ocean  surface  is  in  random  motion,  and  the  ocean  bottom  is  ran- 
domly rough.  Other  choices  may  be  that  the  ocean  is  isohaline  and  adiabatic  (meaning 
no  diffusion  of  salinity,  and  no  heat  transfer).  However,  no  matter  what  choices  are 
made  the  crucial  issue  in  describing  the  background  ocean  is  not  only  fidelity  to  reality, 
but  one  other  thing,  namely  the  internal  consistency  among  all  assumptions.  Clearly  the 
choice  of  a background  ocean  which  is  an  assemblage  of  desired  properties  not  necessarily 
consistent,  will  lend  an  air  of  arbitrariness  to  all  subsequent  calculations  of  signal  excess. 

Of  equal  importance  are  the  boundary  conditions  associated  with  the  seven  equations. 
These  occur  at  (1)  discontinuities  of  surface  and  bottom,  (2)  discontinuities  in  the  volume 
of  velocity,  salinity,  temperature,  pressure  of  mass  density.  Such  conditions  determine  in 
great  measure  the  specific  form  solutions  of  the  equations,  as  they  apply  to  the  ocean, 
and  in  particular  the  sound  speed  profiles. 


OCEAN  PERTURBATIONS  DUE  TO  DRIFTING  INHOMO- 
GENITIES  AND  THEIR  SPACE  AND  TIME  SCALES 

A second  problem  is  the  proper  description  of  ocean  perturbations.  Since  there  are 
several  causes,  the  question  is,  which  perturbations?  Let  us  consider  here  perturbations 
inside  a volume  V caused  by  inhomogenities,  Fig.  2(a).  Choose  one  variable  of  the  seven 
of  the  background  ocean,  or  choose  some  scalar  quantity,  say  sound  speed,  which  is  a 
combination  of  them.  Let  h(x,  t)  represent  the  choice  of  parameter.  A measurement  of 
h actually  made  at  sea  may  show  it  to  be  a random  variable,  Fig.  2(b).  We  ask:  from 
this  record  can  we  determine  the  characteristic  sizes  or  space  scales  of  h?  One  answer  might 
be,  try  Fourier  decomposition,  viz,  at  time  tx , 

htf.tx)  = H(it,  tj  )e‘K*x 

J (2n)3 

H(K,tx)  = f d3xh(x,t1)e'iK-?  . 

Each  sinusoidal  component  has  a spectral  weight  H(K,  tj ) d3K.  This  could  be  a character- 
istic scale.  More  conveniently,  the  ocean  volume  V can  be  pictured  as  a sample  member  of  a 
periodic  train.  Then  H(K,  tj)  is  a line  spectrum,  which  could  also  be  considered  a scale. 
Other  mathematical  measures  of  characteristic  size  are  in  vogue.  But  it  is  more  physical 
to  use  second-order,  or  energy -related,  statistics.  Thus  to  obtain  spatial  scales  one  forms 
the  autocorrelation, 

T(Ax,  tx)  = ^ j'  h(x,  tx)h(x  + Ax,  tj)  d3x  (units:  field2) 

from  which  the  correlation  function  (or  coefficient)  is  defined  to  be 
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^ T(Ax) 

p(Ax,t1)  T(Q)  • 

Several  characteristic  sizes  are  based  on  the  functional  dependence  of  p : for  example,  the 
value  of  separation  Ax  where  p = e-1  is  taken  to  be  the  correlation  length  of  random 
process  h;  or  the  integral  J = /p( Ax,  tj)  d3Ax  (units:  m3)  is  taken  to  be  the  integral 
scale  of  h.  Other  measures  of  characteristic  sizes  are  obtained  by  Fourier  inversion  of  the 
autocorrelation,  which  is  the  “power  spectrum”  of  the  random  spatial  process  h: 

W(k,  tx)  = j" T(Ax,  tjje-'1^'^^  d3Ax  . (units:  field2xm3) 

The  power  spectrum  in  space-time  is  similarly  obtained, 

W(k,  u>)  = J' T(Ax,  r)  exp  (-iK  • Ax  + icor)  d3  Ax  dr  . 


The  spectral  energy  over  a range  of  wave  numbers  from  Kx  to  K2  is 

K2 
Ki 


->  ->  r -»■  ->■ 

S(K1K2,t1)  = / W(K,  tj)  d3K  . 

J tr 


(units:  field2) 


The  magnitude  of  this  number  tells  how  much  energy  is  associated  with  the  Fourier- 
derived  wavelengths  in  the  range  Kj12ir  and  K212tt.  Significance  of  power  spectra  will 
be  taken  up  after  a discussion  of  the  temporal  power  spectrum. 

We  consider  time  scales  by  fixing  point  X!  and  allowing  time  to  vary.  The  entire 
discussion  is  repeated  using  period  J as  analog  of  volume  V,  r as  analog  of  Ax,  and  co  as 
analog  of  2.  One  thus  obtains  the  "temporal  power  spectrum  at  Xx , W(3tj , oo)  (units: 
field2xS),  the  correlation  T(x1,t),  and  the  correlation  length  Lt. 
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Space  and  time  separation  can  be  taken  together,  and  a space-time  power  spectrum 
W(it,  o>)  of  the  perturbation  process  h constructed.  Figure  3 shows  a one-dimensional 
spectrum  of  this  type.  The  significance  of  this  plot  can  be  understood  if  we  take  the 
perturbation  parameter  to  be  in  the  sound  speed  and  the  cause  to  be  the  vertical  displace- 
ments f(z)  of  internal  waves.  An  integration  of  the  spectrum  over  frequency  range 
to  co2  and  spatial  wavenumber  range  Kj , K2  of  internal  waves  then  generates  the  random 
variable  (f2(z)>.  Several  models  of  the  power  spectra  of  displacement  and  of  fluctuations 
in  sound  speed  profiles  are  in  current  use,  mostly  constructed  by  Garrett  and  Munk.  They 
are  random  functions  of  depth  z and  range  R.  From  them  one  can  construct  the  variation 
of  mean  values  of  displacement  with  depth  and  range  which  is  the  basic  component  in 
calculating  the  rms  value  of  the  random  change  in  sound  speed  profile.  These  models 
predict  coherence  lengths  Lv  of  5c/c  in  the  vertical,  and  LH  in  the  horizontal.  Figure  4 
is  a pictorial  statement  of  such  type  perturbations  of  ocean  sound  speed  profile.  It  snows 
schematically  a coherence  length  in  depth  which  is  much  smaller  than  the  coherence 
length  in  range.  The  complexity  of  this  sketch  invites  efforts  at  simple  modeling.  Here 
we  see  the  background  modeled  by  a bilinear  profile,  and  the  perturbations  of  these  pro- 
files as  functions  of  range.  One  can  add  additional  perturbations  as  functions  of  time. 
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Fig.  3— Space-time  power  spectrum 
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Another  important  ocean  parameter  perturbation  is  random  surface  waves.  Here  we  see 
approximate  amplitude  spectra  of  these  waves  for  which  a useful  model  is  that  of  Pierson- 
Moskowitz. 


FLUCTUATIONS  IN  SOUND  SPEED  PROFILES  CAUSED 
BY  SOURCE-RECEIVER  MOTION 

On  selected  time  scales  an  ocean  of  inhomogeneities  may  appear  frozen.  A sound 
wave  generated  by  a fixed  source  arrives  at  a fixed  receiver  with  a determinate  phase. 
Suppose,  however,  the  source-receiver  combination  are  in  relative  motion.  Then  since  the 
perturbation  in  sound  speed  is  random  in  space,  although  frozen  in  time,  the  effect  on  the 
receiver  sensor  is  the  same  as  if  it  were  fixed  and  the  frozen  ocean  convected  past  it  at 
the  speed  of  the  relative  motion.  The  sensor  is  effectively  scanning  the  random  scalar 
field  can  be  constructed.  An  additional  effect  is  the  change  in  the  random  path  over 
which  acoustic  signals  are  propagated.  This  is  discussed  later. 
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ACOUSTIC  PROPAGATION 

When  a source-receiver  pair  is  inserted  in  the  random  ocean  the  propagation  of  the 
acoustic  signal  becomes  randomized.  Several  decades  of  effort  have  been  spent  in  investi- 
gating propagation,  and  many  models  are  currently  in  use.  Let  us  choose  a ray-trace  model 
for  this  review.  Evidently,  the  effect  of  random  inhomogenities  depends  on  the  selected  ray, 
but  a first  question  is,  which  ray ? In  addition,  in  a realistic  ocean  model  a second  question  is 
how  many  rays  must  be  considered  to  contribute  to  the  field  at  a point?  Consider  a point 
source  of  sound  fixed  in  the  volume  V,  and  radiating  a CW  signal.  The  resultant  acoustic 
field  is  governed  by  the  elliptical  equation  of  Helmholtz  in  the  scalar  acoustic  pressure, 

[V2  + k2(x,  t)]  p(x,  t)  = V(x,  t)p(x,  t) , k(x)  = co/c(x,  t) 

in  which  V(x,  t)  is  the  perturbation  in  the  propagation  wave  number  k(?),  which  in  turn 
is  the  perturbation  in  the  sound  speed.  Neglecting  for  a moment  temporal  effects  caused 
by  the  changes  of  the  environment  we  find  the  appropriate  Green’s  function,  or  response 
at  vector  location  y due  to  a unit  delta  function  source  at  x, 


[V2  + k2(x)]  G(x|y)  = 53(x-y). 


This  is  a deterministic  statement  of  rays  propagating  in  space.  It  shows  that  one  or  several 
of  them  reach  y,  depending  on  the  nature  of  k(x).  Suppose  k(x)  = k = constant.  Then 
only  one  ray  from  the  point  source  reaches  from  x to  y along  the  path  of  a straight  line, 
Fig.  5(a).  Since  the  medium  is  here  assumed  isotropic  and  homogeneous  the  Green’s  func- 
tion is  G = eikl*-y|/4ir|j?-'^|.  If,  however,  k depends  on  location  vector  then  several  rays 
originating  a source  vector  it  may  reach  location  vector  y,  Fig.  5(b).  In  this  case  the  field 
at  y is  the  sum  of  n rays, 


n 

Gtfl?)  = £ Gitftf)  . 

We  recognize  this  to  be  a statement  of  a sound  channel.  The  deterministic  rays  1,  2,  3, 
etc.,  can  be  calculated  to  sufficient  approximation  by  geometric  optics.  If  the  source  is 
at  the  origin,  this  means  each  ray  has  the  form, 


Gill'S)  = K$r3)  exp 


Consider  next  that  the  perturbation  V(x,  t)  of  the  medium  is  random  so  that  each  ray  exe- 
cutes a sort  of  randbm  walk  filling  a tube.  Fig.  5(c).  The  acoustic  field  pressure  at  y is  now  a 
sum  of  contributions  of  those  rays  along  zig-zag  tracks  within  tubes  which  connect  ? to  y . 
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p(y)  = £ Pi(y)  • 

i 


An  enlarged  view  of  a ray-tube  is  shown  in  Fig.  5(d).  Along  cross-section  a - a each  ray  is 
identified  by  two  symbols,  q,  y,  that  is,  pi(q;y).  When  q = y,  pi(y;y)  means  the  unperturbed 
ray.  Outside  the  tube  pi(q,y)  vanishes. 


I 


A Green’s  function  associated  with  each  ray  inside  the  tube  is  written  Gj(q;y|o),  in 
which  the  source  of  the  rays  is  at  the  origin.  The  field  at  q couples  with  the  perturbation 
of  sound  speed  V(q),  and  the  scattered  field  is  propagated  from  q to  $ by  G(y|^).  Thus, 
each  ray  tube  contributes  the  geometric  optics  ray  multiplied  by  a complex  phase  factor 
which  sums  up  the  contributions  of  all  zig-zag  rays  within  the  tube  to  the  field  at  f, 

Pi$)  = Gi(jf)3)e* 


X = 


-ir  f d^Gfjn^MVtfiGitf,?^]. 

'l(yio)  J j tube  ray 


This  important  result  is  a generalized  form  of  the  well-known  Rytov  approximation.  The 
field  Gj(y|o)  is  the  deterministic  part  of  the  ith  contributing  ray.  The  function  ex  is  the 
random  part.  One  sees  that  the  random  part  requires  an  integration  of  the  effects  of  the 
random  phase  changes  along  the  ray-tube  paths. 

Returning  to  our  question,  which  ray?  We  must  reformulate  it  and  ask,  which  ray- 
tube?  Secondly,  how  many  ray  tubes? 


OBSERVABLES  IN  ACOUSTIC  PROPAGATION 

Let  us  accept  the  generalized  Rytov  form  for  the  ray-tube.  For  each  ray-tube  we 
ask  the  question,  what  observables  in  the  field  measurement  are  useful  in  categorizing  the 
nature  of  the  field  at  the  receiver?  Write,  for  convenience. 

Pi  = Poe*  = IP  I e,<#>  , X = Xj  + iX2 

in  which  p;  is  the  acoustic  pressure  contributed  by  the  jth  tube,  p0  is  the  deterministic 
ray-tube  contribution,  <p  ;s  a random  phase,  X} , X2  are  quadrature  components  of  the 
integral  factor  X.  A list  of  observables  is: 


(l2)  = ((lo8e  IPl2)2) 

(P?)  = J lp0l2pxp(2(x2)) 

(l)  = loge  1 Po  1 

+ exp^2(x*2))] 

I = ^ (iPl2) 

X 
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II 

X 

1+ 
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XxX2  = -|lm{x2). 

(<t>2)  ■ (x|)  , (01)  = 2<x1x2) 
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Of  particular  interest  in  modeling  are  (X),  and  <ex>.  Some  modelers  assume  the  8c  fluc- 
tuations to  be  Gaussian  distributed  in  space  and  time.  This  they  take  <X>  = 0 and  <ex>  to 
be  el(l/2)<x2  >].  Attention  is  then  focused  on  finding  a model  for  (X(y)2>.  In  addition, 
focus  is  made  on  p2  as  defined  above,  hence  (|X|2>.  Other  modelers  create  an  entire 
theory  on  the  observable, 


(I2)  - (if  = 2((|X|2)  + Re<X2)]. 


This  entity  takes  a prominent  part  in  the  determining  regimes  of  fluctuations,  now  to  be 
discussed. 


II 


DIFFRACTION  PARAMETER  A,  STRENGTH  PARAMETER  $2 
FOR  STRAIGHT-LINE  RAYS 

We  noted  earlier  that  the  deviation  6 c/c  of  sound  speed  is  a random  function  of 
space  and  time,  and  that  by  use  of  autocorrelations  or  power  spectra  we  can  define  three 
correlation  lengths:  vertical  Ly,  horizontal  Lh,  transverse  Ly ; and  one  correlation  time  Lj-. 
These  specify  the  spatial  and  temporal  characteristic  sizes  of  6 c/c. 

A CW  sound  wave  of  wavelength  X passing  through  a field  of  5c/c  is  diffracted  by 
the  inhomogenities.  A question  arises:  how  important  is  diffraction?  One  can  try  to 
decide  by  calculating  a diffraction  parameter  which  is  the  square  of  the  ratio  of  the  radius 
of  the  first  Fresnel  zone  Rp  of  a point  source  at  range  R to  a selected  coherence  length, 
L,  Fig.  6.  In  a homogeneous  ocean,  where  propagation  is  along  straight  lines,  the  ratio  A 
is  a function  of  range.  The  average  A is. 


lAV 


i rR  /Rp(x)\  dx 
R J \ L / 2tt  ' 


6L2k 
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Similarly,  an  acoustic  wave  passing  through  a field  of  6c/c  experiences  a random  change 
of  phase.  Again  the  question  arises,  how  strong  is  the  phase  change?  A convenient  meas- 
ure is  the  ensemble  average  4>2  of  the  square  of  the  integrated  phase  along  a ray,  is 
called  the  strength  parameter.  In  a homogeneous  medium. 


The  product  <I>2A  is  a measure  of  the  combined  effects  of  diffraction  and  random  phase  change. 

In  a realistic  description  of  the  ocean  acoustic  rays  travel  in  curved  paths.  The  cor- 
relation vector  length  L has  a component  Lp  along  the  ray  path.  We  wish  to  calculate  4>2 
and  A.  A typical  case  is  shown  in  Fig.  7.  Ray  1 is  deterministic  from  the  origin  to  jf.  Ray  2 
reaches  a point  3,  distance  a from  Ray  1,  and  instead  of  continuing  to  point  C is  deflected  at 
angle  0 to  pass  on  to  y along  path  4.  The  difference  AS  in  path  length  between  paths  2,  4 
and  path  1 is  a geometric  construct.  In  3-space 

[AS]  - [A] 


in  which  [A]  is  a matrix.  In  one  dimension  A(x)  = dd/da. 

Using  these  definitions  one  finds  that  in  the  case  of  a deterministic  sound  channel, 
an  integration  along  the  ray  path  is  needed  to  find  4>2,  and  A: 


4,2 ' **  jf  11,1  ((! 

C dx  I 


5c  (z  ray)V 


4>2A  = kn 


5c (z  ray )\ 


L (0,zray) 


L (0,zray) 


kol  A|Ly 


i 
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This  says  that  at  each  increment  in  range,  the  segment  at  depth  z-ray  experiences  a con- 
tribution to  given  by  the  integrand  shown.  The  important  point  to  notice  is  that  4>2A 
is  just  the  quantity  (I2>  - <I>2.  It’s  significance  will  appear  later.  Generally,  however,  the 
integration  must  be  done  by  numerical  means.  The  parameters  4>2A  are  crucial  in  estab- 
lishing regimes  of  signal  statistics. 


i 

! 
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SIGNAL  STATISTICS 

In  the  presence  of  a determinate  channeled  ocean  the  received  signal  is 

p(5?,w,t)  = p0(2,  w)e-iwt . 

When  fluctuations  are  included  the  signal  is 

p(2,  w,t)  = p(5t,  co,  t)e"iwt . 

It  is  convenient  to  define  a reduced  wave  function, 

p(x,w,t) 

0 (x,  CO,  t)  = — • 

Po(x,  w) 

As  before,  several  real  functions  of  0 are  observables.  To  list  them  one  writes^ 

i p = |p|e‘0  fin  «//  = fin  |p|  + i 0 

I = I0I2  = Ipl2  = e8nI  • 

The  list  is 

(1)  I (3)  fin  jpl  (5)  Re  0 

(2)  fin  I (4)  0 (6)  Im  0 . 


Let  us  picture  0 as  a m-multi-dimensional  plot  of  space,  time  and  frequency.  Fig.  8(a). 
In  general  there  are  3 space  coordinates;  1 frequency  coordinate;  1 time  coordinate;  and 
perhaps  1 range  coordinate.  At  point  1 , 0 = 0 (zj , tj , « j ),  and  at  point  2,  0(z2,  t2,  w2) 

A similar  chart  for  the  three  variables,  angle  of  arrival,  doppler  shift,  and  pulse  time, 
(which  are  Fourier  conjugates  of  Fig.  8(a))  is  shown  in  Fig.  8(b).  The  important  idea  is 
that  4/  is  a random  variable,  so  that  a complete  description  of  the  m-dimensional  field 
requires  all  n-point  correlations,  n “ 1,  2,  ...»  °°.  It  is  customary  to  use  one-point  and 
two-point  correlations  of  signal  statistics.  A list  follows: 


One-Point  Statistics 

These  are  statistics  constructed  on  measurements  at  a single  point,  say  X! , along  the 
time  axis  in  the  m-multi-dimensional  plot.  For  example,  measurements  of  4>  = 0(xj > t,  coj) 
are  used  to  construct  the  random  intensity  I,  fin  I,  fin  Ipl,  0,  etc.,  as  time  progresses.  An 
important  application  of  one-point  functions  in  detection  is  that  of  Urick,  which  is  now 
discussed. 
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Urick’s  Phenomenological  Model  Based  on  Intensity  Statistics  al,  One  Spatial  Point 

Fluctuations  of  received  signal  are  caused  by  multipath  propagation  of  an  initially 
CW  signal  in  an  inhomogeneous  moving  medium.  No  matter  what  the  cause  the  direct 
path  component  decreases  with  range,  and  the  scattered  or  diffracted  components  in- 
crease with  range. 

A useful  theoretical  model  is  the  probability  distribution  of  the  envelope  of  a sine 
wave  plus  narrow  band  Gaussian  random  noise,  or  the  probability  distribution  of  a con- 
stant vector,  summed  with  a random  vector  whose  x;  y components  are  Gaussian  distrib- 
uted time-functions  (Fig.  8(c)).  Here 

X(t)  = P cos 'P  + x(t) 

Y (t)  = P sin  'P  + y(t) . 

The  solution  of  this  problem  is  well-known  to  be  the  Rician  distribution.  Fig.  8(d).  Here 
P is  the  magnitude  of  the  constant  victor,  and  T is  the  fractional  relative  power  of  the 
random  component, 

_ = RAYLEIGH  POWER 
1 TOTAL  POWER 

in  which  Rayleigh  power  corresponds  to  zero  constant  vector,  and  V is  the  vector  sum, 


Cumulative  distributions,  obtained  by  integration  over  a specific  Rician  distribution  curve 
from  -°o  to  V,  expressed  in  units  of  V2/V2,  where  V is  the  mean  value  of  V over  the  dis- 
tribution, is  shown  in  Fig.  8(e).  A similar  presentation  with  logarithmic  horizontal  scale 
is  shown  in  Fig.  8(f).  A great  number  of  measurements  made  at  sea  have  been  assembled 
and  compared  with  these  theoretical  predictions.  One  can  gauge  from  this  the  popularity 
of  intensity  statistics  at  one  spatial  point. 


Two-Point  Functions 

Measurements  are  made  at  2 points  in  m-space.  They  are  separated  in  space  (Az,  Ay), 
in  time  At  = t2  - tj,  or  in  frequency  Aoj.  One  measures,  for  example, 

\p(Zl,  tlf  cj)  , i//(zx  + Azj,  t + At,  + Ac o)  . 

Cross-correlation,  or  covariance,  of  \p  is  written  as  an  average  over  space,  time  and  fre- 
quency 

(lMzl»  tj,  Wj)^(zi  +Az1,t1  + Atj , cuj  +ACJ!)). 

As  noted  earlier,  for  each  separation  variable  Az,  At,  Aoo,  etc.,  a plot  of  i J/$)  vs 

separation  can  be  used  to  define  a coherence  scale.  Similarly,  by  Fourier  transformation 


if  (z. 
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Fig.  8(e)— Rician  curves  with  T as  parameter 
on  semi-log  coordinates 


Fig.  8(f)— Rician  curves  with  T as  parameter  on  prob- 
ability coordinates  and  a logarithmic  horizontal  scale 
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Thus  the  conjugate  of  separation  in  depth  is  vertical  angle  of  arrival,  etc.  Conjugates 
are  very  useful  in  applied  work. 

Three-Point  Functions 


A cross-power  spectrum  constructed  on  2-points  may  be  a function  of  time.  For 
example,  take  the  time  separation  correlation  <i//i(tj ) ^(tj  + At)>  and  its  cross-power 
spectrum  of  doppler  shift  IP(7).  Measured  over  a given  duration,  T.  Break  T into  several 
smaller  durations,  or  integration  times,  and  construct  W(7)  for  each  of  these.  Plots  of 
Wx(7),  W2(7),  etc.,  show  different  spread  of  W among  them,  and  show  the  central  doppler 
shift  wandering,  Fig.  9.  A description  of  spread  and  wander  can  be  obtained  by  use  of 
three-point  functions. 


Signal  statistics  contained  in  1-,  2-  and  3-point  functions  form  the  basis  for  calcu- 
lating the  coherence  scales  or  their  conjugates  as  they  vary  along  ray  paths.  Thus  path 
lengths  of  disturbed  rays  become  functions  of  A,  <!>. 


SPREAD 


WANDER 


Figure  9 
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Earlier  it  was  stated  that  when  the  medium  was  fluctuating  the  sound  propagates  in 
ray  tubes  rather  than  rays.  Each  tube  contains  a number  of  separate  rays  which  start  at 
the  source  and  end  at  the  receiver.  They  are  micro  multipaths.  Ray  tubes  themselves 
that  meet  at  a receiver  are  multipaths  proper.  The  question  arises:  what  is  the  spread  of 
the  micromultipaths  compared  to  correlation  lengths?  Figure  10  shows  comparisons. 

Take  the  case  of  vertical  separation,  Az,  and  calculate  its  2-point  correlations,  its 
cross-power  spectrum  and  thus,  its  vertical  correlation  length,  Lv.  Also  calculate  the 
strength  parameter  4>  and  the  diffraction  parameter  A.  It  is  found  that  when  A4>  > 1 
the  micromultipaths  are  spread  over  a region  larger  than  Lv  and  are  thus  uncorrelated. 
The  regime  is  called  fully  saturated.  When  A4>2  > 1 the  micromultipaths  are  spread  over 
a region  smaller  than  Lv  and  are  correlated.  The  regime  is  called  partially  saturated. 
When  4>  < 1 there  are  no  extra  micropaths.  The  regime  is  called  unsaturated.  In 
this  regime  calculation  of  ray  statistics  can  be  accomplished  by  small  modifications  of 
geometric  optics. 


A (diffraction) 


Fig.  10  — Regions  in  A - <1>  space 


This  diagram  can  be  useful  in  finding  the  regime  in  which  a particular  sound-transmission 
experiment  is  operating.  Choose  a specific  unperturbed  ray,  then  calculate  4>,  A by  inte- 
grating amplitude  and  phase  accumulation  along  the  ray  path  caused  by  random  fluctua- 
tions in  the  medium.  Then  convert  the  information  into  a frequency  range  plot  of  the 
type  shown  in  Fig.  11. 

The  usefulness  of  A - <1>  space  is  illustrated  for  three  cases:  (1)  time  separation,  Fig. 

12,  in  which  the  tip  of  the  signal  phase  shows  a temporal  history  of  different  type  depend- 
ing on  A,  4>  for  a chosen  ray;  (2)  space  separation.  Fig.  13,  in  which  the  ocean  is  frozen,  and 
2-point  statistics  are  determined  for  each  separation,  and  the  signal  phase  executes  a dance 
as  the  spatial  separation  is  increased;  (3)  pulse  time  wherein  a delta  pulse  emitted  at  the 
source  arrives  at  the  receiver  with  a different  spread  in  time,  depending  on  A,  4>,  Fig.  14. 

Actually  the  JASON  theory  presented  here  has  strong  analogs  in  the  theory  of  FM 
modulation.  Standard  references,  say  Taub  and  Schilling  “Principles  of  Communications 
Systems”  describe  phasor  diagrams  of  wide  band  and  narrow  band  FM  signals  (Fig.  15(b). 
Here  two  sidebands  consituting  a modulation  of  a carrier  form  a resultant  R.  The  prop- 
erties of  this  modulation,  as  discussed  in  this  text,  parallel  concepts  shown  in  Figs.  11, 

12.  and  13. 


A ► 

Fig.  13— Space-separation  behavior  of  the  wave-function  phasor.  The  curves  indicate  a typical 
behavior  of  the  phasor,  where  the  parameter  along  the  curve  is  spatial  separation. 


ANTENNAS  AND  FLUCTUATIONS 

Inhomogeneous  medium  impose  limits  on  the  achievable  directive  gain.  Consider  the 
case  of  a wave  front  which  is  initially  plane,  then  t-ansverses  an  inhomogeneous  medium 
resulting  in  random  distortions  of  phase  and  amplitude.  An  antenna  length  L receives  the 
random  phase  signal.  The  mean  directive  gain  is  written. 

0-5°  jj 

-1 

in  which  G0  is  the  determinate  gain  in  the  absence  of  fluctuations,  x = 2z/L,  and  0 is  a 
random  phase,  mean  zero,  variance  a and  correlation  distance  p.  The  question  is  posed, 
is  there  a practical  limit  to  the  mean  directive  gain  as  the  length  L increases?  The  answer 
depends  on  a,  p,  L.  There  is  a wide  literature  on  the  subject,  contributed  by  radar  and 
underwater  surveillance  research.  NRL  has  worked  on  this  problem  over  several  years. 

A specific  answer  of  radar,  tailored  to  atmospheric  inhomogeneities  is  shown  in  Fig.  16. 
Other  characteristics  of  statistical  antenna  are  mean  beam  width,  mean  side-lobe  levels, 
mean  directivity,  etc.;  variance  and  correlations  in  directive  gain;  beam  width  direction  of 
the  principal  maximum,  etc.  All  these  similarly  depend  on  a,  p,  L. 
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CONCLUSION 
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Performance  criteria  of  underwater  surveillance  systems  in  current  use  are  cumulative 
probability  of  detection  and  mean  time  to  hold  contact.  Both  are  calculated  from  a 
knowledge  of  signal  excess.  In  a fluctuating  medium  where  source-receiver  may  be  in 
motion,  signal  excess  is  a random  variable  caused  by  perturbations  in  ray  paths  due  to 
surface,  bottom,  random  volume  inhomogeneities,  noise,  antenna  characteristics,  operator 
performance,  etc.  The  present  review  has  dealt  with  fluctuations  in  transmission  loss  and 
source-receiver  motions.  The  key  issues  were  couched  in  terms  of  modeling  the  fluctua- 
tion effects.  They  are: 

• which  background  ocean? 

• which  perturbation  of  the  background? 

• which  space  and  time  scales  of  fluctuations,  and  their  scales  compared  to  acoustic 
wavelength? 

• which  micromultipaths  in  a ray  tube? 

• which  combination  of  ray-tubes  (=  multi  paths)? 

• which  observable? 

• which  signal  processing  of  which  observable? 

• which  signal  statistical  regimes  of  the  A - 4>  chart? 

Other  issues  are: 

• coherence  degradation  as  a function  of  platform  motioiynultipath  structure, 
surface  and  bottom  roughness,  acoustic  wavelength,  observation  time. 

• statistical  independence  of  terms  in  the  sonar  equation. 


ADDITIONAL  ISSUES 


The  key  issues  in  research  on  acoustic  fluctuations  caused  volume  scattering  just  dis- 
cussed are  relevant  to  time  scales  of  hours  to  days.  On  a shorter  time  scale,  say  less  than 
one  hour,  the  dominant  cause  of  acoustic  fluctuations  is  source/receiver  motion.  Over  long 
enough  ranges,  such  motion  makes  the  received  signals  show  time-varying  multipath  inter- 
ference. Coherence  length  and  coherence  time  become  parameters  for  characterizing  the 
length  scale  and  time  scale  over  which  these  signals  can  be  effectively  averaged.  Degradation 
of  coherence  between  separated  receiving  sensors  as  range  and  time  increase  is  a frequently 
observed  phenomenon  for  high  frequencies,  but  more  rarely  observed  for  very  low  frequen- 
cies in  deep  ocean  paths.  Thus,  a key  issue  in  acoustic  fluctuation  research  is  coherence 
degradation  as  a function  of  platform  velocity,  multipath  structure,  surface  and  bottom 
roughness,  acoustic  wavelength,  and  duration  of  observation. 
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The  results  of  research  on  acoustic  fluctuations  in  the  ocean  have  accumulated  in  the 
last  decade  in  various  scientific  journals,  government  laboratory  reports  and  personal  notes. 
Only  a portion  of  these  has  reached  potential  users  in  systems  applications.  The  unused 
portion  is  large  and  clearly  needs  to  be  applied  to  practical  problems  in  submarine  detection. 
However,  even  if  all  of  it  were  used  there  would  still  be  unresolved  issues  requiring  further 
research.  In  this  connection,  a key  problem  remaining  is  the  statistical  interdependence  of 
terms  in  the  sonar  equation  for  signal  excess.  For  example,  the  mutual  coupling  of  ambient 
noise  and  directivity  is  recognized  to  be  significant.  Further  analysis  and  measurement  are 
needed  to  fix  magnitudes  of  this  coupling,  and  other  coupling  of  terms  in  the  power 
equation,  and  to  evaluate  their  importance  in  detection  algorithms. 
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PSEUDO-SURVEILLANCE  SYSTEM  PERFORMANCE  SIMULATOR 


L.  Fretwell 
(See  Volume  2) 


FLUCTUATIONS  FOR  THE  APAIR , APSURF,  APSUB  MODELS 


R.  Flum 
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REVIEW  OF  BEAM-NOISE  FLUCTUATION  MODELS 

R.  C.  Cavanagh 
Science  Applications,  Inc. 

McLean,  Virginia  22101 


SUMMARY 

For  system  performance  estimation  and  the  interpretation  of  measured  data,  there  is  a 
need  to  predict  and  simulate  the  statistical  properties  of  the  response  of  a horizontally- 
directive  acoustic  array  to  ambient  noise.  In  current  use  in  the  Navy  there  are  a number  of 
such  “beam-noise”  models  which  predict  the  properties  of  noise  caused  by  the  dominant, 
prevailing  source  at  low  frequencies:  surface  shipping.  A study  was  performed,  with  ob- 
jective: to  review  these  models  in  some  detail  and  then  to  recommend,  for  specific  applica- 
tions, approaches  which  utilize  the  best  features  of  each  model. 

Nine  Navy -sponsored  models  were  thus  investigated,  and  all  were  found  to  be  based  on 
the  fundamental  hypothesis  that  the  beam  noise  is  the  convolution  of  the  array  beam  pat- 
tern with  the  sum  of  intensities  from  the  individual  ship  sources.  The  models  are  then 
distinguishable  by  their  treatments  of  the  ship  sources,  the  transmission  loss  (TL),  and  the 
array’s  response.  Two  categories  were  natural:  Analytical  models  which  calculate  the  sta- 
tistical properties  of  noise  directly  from  those  of  the  components  (source  level,  TL,  etc.); 
and  Brute-Force  models  which  use  simulation  or  Monte  Carlo  techniques.  The  ability  to 
calculate  statistics  over  all  possible  values  of  a parameter  (e.g.,  ship  locations  for  a given 
density,  all  possible  source  levels),  or  over  values  which  might  occur  in  a short  time  period 
(say  24  hours)  led  to  a second  type  of  classification:  “grand”  and  “short-term”  ensembles. 

The  nine  models  have  been  reviewed  and  are  described  in  a final  report*  in  terms  of 
approach,  treatment  of  ships  and  TL,  receiver  submodel,  implementation,  advantages,  and 
shortcomings.  Tables  jure  given  to  summarize  these  findings. 

For  various  applications  of  beam-noise  models,  the  statistical  quantities  and  types  of 
ensembling  required  have  been  proposed  (e.g.,  first-order  statistics  with  full  ensembles,  or 
multi-array  statistics),  and  the  most  promising  models  and  techniques  are  identified  for 
each.  General  recommendations  for  an  advanced-development  approach  are  given. 

The  work  was  supported  by  the  Long  Range  Acoustic  Propagation  Project  (NORDA 
Code  600). 


•Review  of  Models  of  Beam-Noise  Statistics,”  Final  Report  to  LRAPP,  SAI-78-696-WA;  Ocean  Science 
Division,  Science  Applications,  Inc.,  McLean,  VA;  November  1977. 
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EXPLANATORY  NOTES  AND  THE  FIGURES 

1..  Nine  beam-noise  fluctuation  models  were  reviewed.  Fig.  1. 

2.  In  general,  the  models  concentrate  on  certain  fluctuation  mechanisms,  as  shown.  Fig. 
2.  No  model  explicitly  treats  all  of  these,  but  the  list  gives  the  ones  most  often  treated  in 
rough  order.  Although  wind-related  noise  can  be  very  important,  all  of  the  models  concen- 
trate on  the  surface-ship  component.  Take  note  of  the  time  scales. 

3.  Before  discussing  how  the  models  predict,  consider  what  they  predict,  Fig.  3.  For 
each  application,  as  shown,  specific  beam-noise  statistics  are  required.  In  defining  such  sta- 
tistics, care  must  be  taken  to  specify  a time  scale.  There  can  be  a substantial  difference  be- 
tween temporal  statistics  assembled  over  a day  and  over  a week  or  season. 

4.  To  illustrate  the  last  point,  histograms  of  modeled-beam  noise  are  shown,  Fig.  4. 
Each  corresponds  to  a 12-hour  sample,  but  the  three  are  taken  over  one  week’s  period  of  time. 
The  histogram  corresponding  to  one  week’s  ensemble  is  the  average  of  the  individual  histo- 
grams, and  here  would  be  much  broader  than  any  of  the  short-term  components.  The  sug- 
gestion is  that  it  would  be  a mistake  to  use  a long-term  (or  “grand”)  ensemble  distribution  of 
noise  in  a model  to  predict  such  a Measure  of  Effectiveness  (MOE)  as  the  24-hour  cumulative 
detection  probabilities. 


5.  Consider  next  how  the  models  calculate  noise  fluctuations:  Formula  (1-1),  Fig.  5,  is 
common  to  all  of  the  models  reviewed.  If  each  variable  were  known  and  deterministic,  then 
a standard  “point”  model  such  as  RANDI  would  result.  The  form  of  the  equation  should  not 
suggest  that  the  factors  are  uncoupled— for  in  most  treatments  at  least  the  array  given  depends 
on  transmission  angles.  No  model  actually  uses  the  time  dependence  of  SL  nor  performs  the 
sum  in  phase. 

Models  are  thus  distinguished  by  their  treatment  of  the  ship  source  field  and  by 
the  method  of  calculating  the  statistics  of  Noise  (N). 


6.  To  help  differentiate  model  types, two  approaches  are  defined,  Fig.  6.  The  term 
“Brute-Force”  is  not  meant  to  suggest  that  the  approach  is  unsophisticated  or  mathematically 

trivial,  nor  is  “Analytic”  meant  to  convey  more  than  the  notion  that  the  calculation  is  per- 
formed without  direct  realizations  of  the  random  variables. 

The  first  four  models  of  Fig.  1 are  Analytic,  the  others  Brute  Force. 

7.  Consider  next  some  statistical  quantities.  Fig.  7,  and  how  they  might  best  be  pre- 
dicted by  the  models  reviewed. 

• First  order  statistics,  with  long-term  ensembling,  can  be  most  efficiently  pre- 
dicted by  an  Analytic  model. 


• Statistics  for  limited-time  ensembles  are  difficult  to  generate  with  Analytic 
models,  but  are  natural  for  the  Brute-Force  models. 


• Of  the  Analytic  models,  only  the  BTL  model  is  capable  of  predicting  the  auto- 
covariance function  for  the  general  case.  Certain  of  the  Brute-Force  models 
routinely  predict  it. 


64 


\ * 


! 


• For  higher  order  temporal  statistics  (such  as  level-crossing  properties)  and/or 
limited-time  ensembling,  the  Analytic  approaches  can  seldom  be  applied.  Hence, 
a Brute-Force  model  is  indicated. 

• Multi-beam  and  multi-site  statistics  (e.g.,  cross-correlations)  are  predicted  by 
several  Brute-Force  models,  but  only  by  one  Analytic  approach  (BTL). 

In  summary,  no  single  model  will  satisfy  all  requirements;  a carefully  combination  of 
the  best  of  each  type  seems  the  best  approach. 

8.  Selected  problem  areas  for  a beam-noise  prediction  capability  are  noted.  Fig.  8. 


A - USI  Model  (Underwater  Systems,  Inc.) 

B — BTL  Model  (Bell  Laboratories) 

C — BBN  Model  (Bolt,  Beranek  and  Newman,  Inc.) 
D — WAGNER  Model  (Wagner  Associates) 

E - NABTAM  Model  (ORI,  USI,  NORDA,  et  al) 

F - OSBN  Model  (Science  Applications,  Inc.) 

G - BEAMPL  (NORDA) 

H - SIAM  I (NRL) 

I -SIAM  II  (NRL) 

Figure  1 
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Mechanism 

• Ships  on  and  off  beams 

• Average  TL  variations 

• Detailed  TL  fluctuations 

• Array  response  variations 

• Ship  source  level  variations 

• Wind  noise  variations 


Typical  Time  for  a 
Significant  Change 

Minutes  for  Nearby  Ships,  Hours 
to  Days  for  Distant  Ships 

Hours 

Minutes  to  Hours 
Minutes  to  Hours 
Hours  to  Days 
Hours  to  Days 


Figure  2 
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SAMPLE  PROBABILITY 


TYPICAL  APPLICATIONS 


• Instantaneous  PD 

• CPD  over  hours  or  a day 

• Distribution  of  Holding  Times  over  a Day 


BEAM-NOISE  STATISTICS 

• Distribution  of  Beam-Noise  Level  (over  a Short  Time  Period  and  Small  Ocean  Region) 

• Statistical  Dependence  of  Noise  at  Two  Points  in  Time 

• Distribution  of  Beam-Free  Times 

• Distribution  of  Time  Intervals  for  Which  Noise  is  Below  a Threshold 

Figure  3 


Figure  4 
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N(t)  =2^  SLj(t)  ■ Tj  (t)  • AGj(t), 


J(t)  is  the  number  of  ships  at  time  t, 

SKj(t)  is  the  source  intensity  of  the  j-th  ship  at  time  t, 

Tj(t)  is  the  transmission  ratio  for  ship  j at  time  t (i.e.,  TL  = -10  log  (Tj)), 

AGj(t)  is  the  array  response  for  the  arrivals  from  ship  j at  time  t. 

Figure  5 


BEAM-NOISE  FLUCTUATION  MODELING 

• Analytic  Approach:  - treats  variables  as  well-defined  random  processes 

- calculates  statistical  properties  "directly" 

- typically  employs  convolutions 

• can  yield  efficient  model  for  predicting  selected  noise  statistics 

• is  difficult  to  set  up  and  to  do  higher-order  statistics 

• Brute  Force  Approach:  - evaluation  of  statistics  is  numerical,  using  realizations  of  ship  tracks,  etc. 

• can  be  used  to  estimate  just  about  any  statistic 

• can  use  all  available  information 

• is  often  very  expensive  and  requires  special  analysis  routines 

Figure  6 


STATISTICAL  QUANTITIES 
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• One-Dimensional  Density  Functions  and  Moments,  Grand  Ensemble 

• One-Dimensional  Density  Functions  and  Moments,  Limited-Time  Ensemble 

• Basic  Temporal  Statistics,  Grand  Ensemble 

• Higher-Order  Temporal  Statistics 

• Multi-Beam  and  Multi-Array  Statistics 

Figure  7 


SOME  PROBLEMS 


• Inputs  for  Detection  Simulations:  Need  approach  to  predict  signal-plus-noise 

in  one  "Bin"  and  noise  in  another 

• Statistics  for  Time-Dependent  Array  Response:  None  of  the  models  treat  a 

moving  array. 

• Ship  Information:  In  spite  of  some  recent  progress,  still  lack  information  on 

source  levels  and  locations 

• Model  Evaluation:  Convincing  evaluations  of  models  for  mean  noise  are  scarce. 

• Wind-Dependent  Noise:  Array  system  response  to  wind-generated  noise  is  not 

understood. 


Figure  8 
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THE  EFFECT  OF  ACOUSTIC  FLUCTUATIONS  IN  THE  OCEAN  UPON  COHERENCE 


(See  Volume  2) 


THE  EFFECTS  OF  FLUCTUATING  SIGNALS  AND  NOISE 
ON  DETECTION  PERFORMANCE 

J.  C.  Heine  and  J.  R.  Nitsche 
Bolt  Beranek  and  Newman  Inc. 

Arlington,  Virginia  22209 


ABSTRACT 

Receiver  operator  characteristics  derived  for  post  detection  integration 
receivers  are  significantly  perturbed  if  the  signal  power  fluctuates  with  a charac- 
teristic time  roughly  two  or  more  times  the  receiver  integration  time.  In  this 
paper,  the  effects  on  detection  performance  of  slow  fluctuations  in  both  signal 
power  and  in  noise  power,  such  as  occur  in  the  10-200  Hz  region  due  to  mer- 
chant shipping,  are  explicitly  considered  and  are  compared  to  results  for  signal 
fluctuations  alone.  Furthermore,  degradation  in  detection  performance  in  the 
region  of  PD  > 0.5  is  shown  to  occur  and  the  sensitivity  of  the  degradation  to 
the  statistics  of  the  noise  power  fluctuations  is  demonstrated.  The  effects  of 
the  noise  statistics  on  experimental  determination  of  sonar  system  performance 
is  discussed. 


INTRODUCTION 

For  discussion  purposes,  the  system  will  consist  of  an  ideal  single  channel  quadratic 
detector  which  is  preceded  by  a narrow  band  filter  and  followed  by  an  ideal  integrator. 
The  input  to  this  temporal  processor  is  the  output  from  a single  beam  of  an  underwater 
acoustic  array. 

One  measure  of  the  detection  performance  for  such  a system  is  obtained  from  the 
system’s  receiver  operator  characteristic  (ROC)  curve.  Classically,  the  system’s  ROC 
curves  are  obtained  in  terms  of  ideal  inputs;  namely,  stable  sinusoidal  signals  in  white 
Gaussian  noise. 

It  is  well  known  that  the  real  ocean  environment  is  not  characterized  as  white  Gaus- 
sian noise.  Moreover,  because  of  multipath  propagation  primarily,  the  signals  are  not 
received  as  stable  sinusoids. 

In  this  paper,  we  examine  the  impact  of  fluctuating  signals  and  noise  (in  short,  fluc- 
tuating output  signal  to  noise  ratios)  on  system  performance  prediction. 

In  particular,  we  are  herein  concerned  with  the  effects  of  signal  and  noise  fluctua- 
tions which  have  characteristic  times  on  the  order  of  two  or  three  times  the  receiver 
integration  time. 

In  the  body  of  the  paper,  we  first  identify  the  sources  of  fluctuating  output  signal  to 
noise  ratios  (SNR’s)  and  estimate  their  characteristic  times. 

The  transition  curve,  which  is  here  defined  to  be  a plot  of  detection  probability  (PD) 
vs  SNR  for  a given  false  alarm  probability  (PFA)  and  integration  time  (T),  is  one  means 
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for  assessing  system  performance.  In  order  to  establish  a basis  for  comparison,  we  next 
review  the  development  of  the  ROC  or  transition  curve  for  the  constant  signal  case. 

Our  method  of  approach,  in  the  manner  of  A.  D.  Whalen  [1]  * is  to  then  consider  the 
SNR  as  a random  variable  so  that  PD  is  a function  of  a random  variable.  Given  the  prob- 
ability density  function  of  SNR,  namely  p (SNR),  we  can  determine  a set  of  modified 
transition  curves  which  will  be  average  Pj>  (i.e.,  E[PD])  vs  E[SNR].  The  principal  task 
then  is  to  determine  the  forms  of  p (SNR). 


Finally,  we  interpret  the  resulting  set  of  transition  curves  in  the  context  of  system 
performance  prediction. 


RECEIVER  OUTPUT  FLUCTUATIONS 


We  shall  consider  a single  channel  NB  filter-square-law  detector-integrator  as  our  tem- 
poral processor  or  receiver  as  sketched  in  Fig.  1.  To  determine  whether  a signal  is  present, 
the  output  of  the  integrator  is  compared  to  a threshold  VT.  The  level  of  the  threshold  is 
established  on  the  basis  of  a desired  false  alarm  probability  (Pp^)  constraint. 


The  receiver  input  r(t)  will  be  assumed  to  be  obtained  from  a single  beam  of  an 
underwater  acoustic  array.  For  now,  assume  that  r(t)  is  ocean  ambient  noise  and  that 
we  are  interested  in  the  general  spectral  range  of  10  to  200  Hz. 

The  receiver  output  Z(t)  plotted  as  a function  of  time  would  typically  appear  as 
sketched  in  Fig.  2.  As  indicated  in  Fig.  2,  several  noise  powe'r  fluctuation  time  scales  are 
often  observable.  These  characteristic  times  are  identified  as: 


a short-term  fluctuation  associated  with  the  response  of  the  receiver  (controlled 
by  the  input  filter  bandwidth) 


a medium-term  fluctuation  associated  with  the  motion  of  a specific  ship  popu- 
lation in  the  ocean  space  sampled  by  the  spatial  processor 


tl  = a long-term  fluctuation  associated  with  the  rate  of  passage  of  ships  in  the  area. 


These  characteristic  times  are  listed  in  Table  1 


THRESHOLD  (VT) 
COMPARATOR 


Fig.  1— Single  channel  receiver 


•References  are  identified  at  the  end  of  the  paper 


r HRF  . 

— v 1 y — 

- LPF  f(t)  i/*Tf (t)dt  . 

r(t)  • 
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TIME  [hours] 

Fig.  2— Characteristic  fluctuation  times  of  receiver  output 


Table  1 —Characteristic  Fluctuation  Times 


Characteristic  Fluctuation 
Times 

Source  of  Fluctuation 

Typical  Times 
(sec) 

ts  Short  [S] 

Filter  response 

1-10 

tm  Medium  [M] 

Moving  ships  (e.g.,  TL 
and  R) 

102-103 

rL  Long  [L] 

Different  ships  at  different 
look  times 

104-105 

Tjjl  Ultra  long  [UL] 

Seasonal  fluctuations  in  TL 

Seasonal 

The  latter  three  noise  processes  are  random  processes.  The  ones  of  particular  interest 
are  tm  and  rL  and  especially  rM  because  it  is  on  the  order  of  the  receiver  integration 
time.  It  has  been  shown  (2]  that 


where 


rM  ~ 5/R 


5 = correlation  length  of  the  transmission  loss  condition 
R = range  rate  of  passing  ship  relative  to  the  spatial  sensor. 

Heine  [3]  has  estimated  the  probability  density  function  for  Z(t)  (see  Fig.  1)  for  such  an 
M-term  noise  process  and  we  shall  use  it  shortly. 

With  regard  to  the  signal,  if  a sine  wave  of  constant  amplitude  is  transmitted  to  the 
underwater  acoustic  sensor  especially  over  the  appreciable  ranges  of  interest,  the  received 
signal  amplitude  will  not  usually  be  constant.  The  random  amplitude  variations  at  the 
receiver  input  are  primarily  due  to  multipath  acoustic  propagation.  If  the  number  of 
paths  is  sufficiently  large,  the  received  signals  will  look  very  much  like  a sample  function 
of  a narrow-band  Gaussian  noise  process  [4] . The  envelope  of  such  a process  has  a 
Rayleigh  probability  density  function  (pdf)  and  the  signal  is  termed  a Rayleigh  fading 
signal. 


RECEIVER  OUTPUT  CHARACTERISTIC  (ROC)  CURVE 
FOR  STATIONARY  SIGNALS 

The  ROC  curve  (or  alternatively,  the  transition  curve)  is  one  of  the  more  useful  ways 
of  estimating  receiver  detection  performance.  Because  it  is  used  to  define  the  recognition 
differential  (NRD)  and  because  it  will  serve  as  the  basis  for  comparison,  we  will  briefly 
review  the  development  of  the  ROC  for  the  case  of  detecting  a stable  sinusoidal  signal  in 
the  presence  of  additive  stationary  white  Gaussian  noise.  For  this  case,  the  input  to  the 
receiver  is  given  as 

r(t)  = A cos  (oo0t  + 0)  + n(t) 

A = constant  amplitude 
coq  = constant  carrier  frequency 

P(A)(*>  = 0 < 6 < 2rr. 

The  noise  is  white  Gaussian  (with  a power  density  spectrum  Sn(f ) = N0/2;  < f < °°) 

and  it  can  be  represented  in  quadrature  form  [5]  as 

n(t)  = x(t)cosco0t  + y(t)sinoj0t. 

A detection  decision  is  made  every  T-seconds  at  the  output  of  the  integrator.  The 
integrator  output  at  the  sample  time  is 


j..- 
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zo<T)  = b Jo  [a2  * 2Ax(t)1  dt  + b J tx2(t)  + y2(t)1  dt  • 

With  regard  to  the  decision  process,  there  are  two  hypotheses: 

Ho  noise  alone  is  present 
H}  signal  plus  noise  are  present. 

When  Z0t  exceeds  the  threshold  VT,  we  decide  that  signal  is  present;  when  VT  is  not 
exceeded,  we  decide  that  signal  is  absent. 

Since  the  time-bandwidth  product  of  the  receiver  is  large,  we  assume  that  Z0(T)  is 
Gaussian  by  the  Central  Limit  Theorem.  We  can  then  express  two  probability  density 
functions  of  the  output,  depending  upon  the  hypotheses. 


PH»,Z)  ‘ 757,  -[-§] 
P">(Z)-  vsb; exp  [- ’ 


where 


• We  made  a change  of  variables  Z(t)  = Z0(t)  - Rx(0)  for  convenience  where 
Rx(0)  * total  noise  power 

• pHo(Z)  and  pHl  (Z)  are  termed  likelihood  functions 

• Am  is  the  difference  in  the  means  of  these  two  pdf’s. 

The  receiver  output  pdf’s  as  defined  above  are  depicted  in  Fig.  3.  The  probabilities  of 
detection  and  false  alarm  are  obtained  by  suitably  integrating  the  respective  likelihood 
functions. 


Z - zc(t)  - Rx(0) 


no  signal  present 


signal  present 


Z • V.  • 

T O 


Fig.  3 — Receiver  output  probability  density  functions 
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r°°  r°°  i r 721 

'A‘40PH«<z)dz‘40^o“pr^jdz 

JD  = / PHl(Z)dZ=  / — ^ exp[-(-Z~^m)2]  dZ. 

Ka,  ^Kai  ox  a2 


It  may  be  shown  that 


Am  = S = A2/2  = average  signal  power 
2 N0B 

og  = = var  [Z(T)/H0] 


U1  ^ 

By  making  the  change  of  variables 


2 = ^0  + N|B 


and  tl-4 
1 2o?  2 2 o§ 


we  may  obtain  the  forms 


i r Ko0  - s 

P°'2erf' 


We  can  further  express  PD  in  terms  of  the  output  SNR  where  we  define 


S _ S 


average  signal  power 


N0  N noise  power  per  unit  bandwidth 


1 r>/K  B - VT  R 


PD  = ~2  erfc 


|,4R  + 2B 


A plot  of  PD  vs  R is  termed  a transition  curve.  Such  a curve  for  the  receiver  of  interest 
is  shown  in  Fig.  4 for  the  nonfluctuating  signal  case.  The  parameter  K is  chosen  on  the 
basis  of  a specified  acceptable  PFA . As  an  aside,  it  is  pointed  out  that  K is  not  a very 
strong  function  of  PFA  as  noted  in  the  following  table. 


■ • i 
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Figure  4 is  labeled  as  average  PD  vs  average  R in  anticipation  of  the  fluctuating  R 
cases  to  be  considered  next.  For  the  constant  signal  (thus  constant  R)  case,  of  course  no 
averaging  was  done. 

Finally,  before  we  leave  this  case,  it  is  well  to  point  out  that  the  PD[R]  that  was 
derived  is  not  restricted  to  the  small  signal  or  low  R case.  The  forms  are  slightly  different 
as  noted  below 


pd[R1 


= 2 erfc 


K y/B  - y/T 


small  R case . 


ROC  CURVES  FOR  FLUCTUATING  SIGNAL  TO  NOISE  RATIOS 

In  what  follows  we  shall  consider  the  receiver  output  signal  to  noise  ratio  to  be  a 
random  variable.  Consequently  PD[R]  will  be  a function  of  a random  variable  and  the 
expected  or  average  PD  will  be 


E[Pd(R)]  = f PD(R)p(R)dR, 


where 


E[  ] = the  statistical  expectation  operator 

•I  [K  JB  - R' 

Pd(R)  = 2 erir  ^4R  + 2B 


p(R)  = the  probability  density  function  of  R. 

Thus,  the  transition  curves  for  the  fluctuating  R case  will  be  the  standard  transition  curve 
for  nonfluctuating  R averaged  over  all  R or  smoothed  by  p(R).  The  problem  is  then  to 
determine  the  appropriate  analytic  forms  tor  p(R). 

We  are  herein  concerned  with  medium  term  (see  Table  1)  fluctuations  in  R.  That  is, 
fluctuations  in  R that  are  such  that  R can  be  considered  to  be  constant  over  any  T second 
interval,  but  varying  from  interval  to  interval.  The  fluctuations  in  R can  be  due  to  fluc- 
tuating signal  only  or  to  a combination  of  fluctuating  signals  and  fluctuating  noise.  We 
shall  compare  the  effects  of  both. 


Fluctuating  Signal  Only 

The  case  of  a slowly  fluctuating  signal  has  been  considered  by  a number  of  investi- 
gators [1,4,6] . The  amplitude  A of  a Rayleigh  fading  signal  has  a probability  density 
function  (pdf)  given  by 


t 


77 


where 


r 

i* 


P(A) 


a r a2 

= ¥ exp  [-  2S 


0 < A < «>, 


[a21 

S = E T = average  signal  power. 
The  average  PD  is  thus 


1 

rv  a21 

K°0  - T 

A 

A2' 

2 erfc 

_ V?  °i 

• -g  exp 

2S. 

dA  , 


and  in  terms  of  the  average  signal  to  noise  ratio  (R  = A2/2N),  we  have 

'K  VB  - v'T  R' 


/•  oo 

E[Pd(R)]  = J |erfc 

— OO 

where  R = E[R] . 


L V4R  + 2B  J R 


1 e'R/R  dR  , 


The  average  PD  vs  average  R is  plotted  in  Fig.  5 and  the  standard  transition  curve  is 
shown  also  for  comparison.  We  can  interpret  the  new  transition  curve  as  a smoothing  of 
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the  standard  curve  by  p(R).  In  this  case  p(R)  is  an  exponential  pdf  and  is  shown  plotted 
in  Fig.  6 for  various  values  of  R.  Whereas  the  standard  transition  curve  shows  a rather 
abrupt  “transition”  from  low  to  high  probabilities  of  detection  (i.e.,  over  a AR  = 3 dB), 
the  smoothed  curve  shows  the  same  transition  (i.e.,  from  PD  = 0.2  to  0.8)  occurring  over 
a considerably  greater  range  of  SNR,  namely  AR  = 9.5  dB. 


Fluctuating  Signal  and  Noise 

We  now  extend  the  previous  development  by  considering  a random  signal  to  noise 
ratio  due  to  fluctuations  in  both  signal  and  noise.  Again  we  restrict  our  attention  to  the 
i characteristic  fluctuation  time  rM  ~ 102  - 103  seconds. 

The  signal  power  fluctuations  will  again  be  assumed  to  be  exponential.  The  noise 
power  fluctuations  have  been  hypothesized  as  being  gamma  distributed  [3]  with 


X_ 

P(N)  = FIT]  <XnN>  exp  [-XnN]  , 


r = E2[N]/var  [N]  = half  the  number  of  degrees  of  freedom 
K = r/E[N]. 

The  signal  to  noise  ratio  R = A2/2N  can  be  obtained  by  a direct  transformation  since  we 
have  p(A)  and  p(N)  and  they  are  assumed  to  be  independent.  Also  since  R is  the  ratio 
of  two  Chi-squared  pdf’s,  R will  have  an  F-pdf. 


p(R)  = 


[l  + 41' 

rRj 


A family  of  these  pdf’s  is  shown  in  Fig.  7.  They  are  plotted  parametrically  for  various 
mean  signal  to  noise  ratios  (R)  and  for  2 and  6 degrees  of  freedom.  The  average  Pp  for 
the  case  of  fluctuating  signals  and  noise  is  thus 

K v/B  - VT  R 


4R  + 2B 
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Fig.  7— Probability  density  functions  of  R for  fluctuating  signals  and  noise 
for  various  R and  degrees  of  freedom 


AVERAGE  PROBABILITY  OF  OETECTIOR  E[P0(R)] 
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A set  of  transition  curves  (for  2,  6,  and  10  degrees  of  freedom)  for  the  case  of  fluctuating 
signals  and  noise  is  shown  in  Fig.  8.  The  standard  curve  is  shown  for  comparison.  As  was 
true  for  the  fluctuating  signal  alone  case,  we  can  interpret  these  new  transition  curves  as 
smoothed  or  averaged  versions  of  the  constant  version. 

Observe  that  as  the  number  of  degrees  of  freedom  n increases,  the  fluctuating  signal 
and  noise  curves  approach  the  fluctuating  signal  alone  curve.  This  asymptotic  behavior 
for  large  n can  be  understood  by  recalling  that  p(N)  was  assumed  to  be  gamma  and  that 
its  variance  decreases  as  N increases.  In  the  limit  N large,  the  density  p(R)  becomes  ex- 
ponential; as  is  p(S)  for  signal  fluctuations  alone.  Figure  9 compares  the  exponential  pdf 
(appropriate  to  fluctuating  signal  alone)  with  the  F-density  with  two  degrees  of  freedom 
(appropriate  to  the  extreme  case  of  signal  and  noise  fluctuations). 


IMPLICATIONS  ON  PERFORMANCE  PREDICTION 

We  have  observed  that  the  transition  curves  for  the  fluctuating  signal  and  noise  case  are 
smoothed  versions  of  the  constant  amplitude  case.  The  amount  of  smoothing  was  seen  to  be 
a function  of  the  numbers  of  degrees  of  freedom  of  p(R).  To  demonstrate  the  possible  im- 
pact that  such  smoothing  could  have  on  system  performance  prediction,  consider  Fig.  10. 

A pair  of  dashed  and  solid  curves  were  calculated  for  two  different  spatial  sensors; 
one  spatial  sensor  provided  a theoretical  5 dB  improvement  in  average  signal  to  noise 
ratio,  R.  Assume  this  improvement  was  due  solely  to  the  fact  that  the  acoustic  aperture 
of  one  array  was  three  times  the  aperture  of  the  other.  One  curve  in  a pair  assumes 
constant  R;  the  other  assumes  a fluctuating  R due  to  fluctuating  signal  and  noise. 


Fig.  8— Transition  curves  for  the  fluctuating  signal  and  noise  case 
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Fig.  10— Comparison  of  transition  curves  for  increased  array  length 


Consider  the  predicted  improvement  in  detection  performance  as  would  be  calculated 
under  the  constant-R  assumption.  This  is  shown  in  the  upper  half  of  Fig.  11.  Note  that 
at  R 10  dB,  the  PD  for  the  long  array  is  the  PD  for  the  short  array  plus  80%;  a sig- 

nificant improvement  in  detection  performance.  Furthermore,  an  improvement  of  about 
40  percentage  points  is  obtained  over  the  range  -12.5  dB  < R < -7.5  dB. 

Now  consider  the  predicted  improvement  for  the  fluctuating  R case.  The  results  are 
shown  in  the  lower  half  of  Fig.  11.  For  this  case,  an  improvement  is  seen  for  all  R,  but 
the  improvement  is  an  almost  constant  20  percentage  points.  It  is  not  immediately  clear 
that  such  an  improvement  is  significant,  especially  when  viewed  in  terms  of  the  overall 
costs  of  increasing  the  array  length  by  a factor  of  three. 


SUMMARY 

We  have  discussed  the  effects  of  fluctuating  output  signal  to  noise  ratios  (R)  on  sys- 
tem performance  and  performance  prediction.  The  effects  due  to  fluctuating  signal  alone 
were  compared  to  those  due  to  fluctuating  signal  and  noise.  The  fluctuations  considered 
were  assumed  to  have  characteristic  times  on  the  order  of  two  or  more  receiver  integra- 
tion times  and  we  confined  our  attention  to  the  10-200  Hz  spectral  range.  The  p(R)  for 
fluctuating  signal  alone  was  exponential  and  the  p(R)  for  fluctuating  signal  and  noise  was 
F-type. 
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Fig.  11— Comparison  of  improvements  in  detection  probability  with  increase  in  array  length 


The  following  summary  statements  are  made  in  the  context  of  the  stated  assumptions. 

• Under  the  assumptions  of  constant  signal  in  white  Gaussian  noise,  the  transition 
curve  is  very  sharp;  i.e.,  0.2  < PD  < 0.8  for  a AR  = 3 dB.  Hence  this  would  indicate  a 
strong  emphasis  on  accurate  measurements  of  mean  values.  Actual  transition  curves  are 
considerably  more  smooth,  implying  a need  for  measuring  the  statistics  of  R. 

• The  fluctuating  R,  regardless  of  the  source  of  fluctuations,  results  in  an  improved 
detection  performance  at  low  R (i.e.,  long  ranges)  and  decreased  performance  at  shorter 
ranges.  Stated  in  other  terms,  the  R required  at  a given  range  for  the  fluctuating  R is  on 
the  order  of  5 dB  more  than  that  required  for  nonfluctuating  R at  a PD  = 0.8. 

• The  fact  that  the  fluctuating  signal  and  noise  case  approaches  the  fluctuating 
signal  alone  case  as  the  number  of  degrees  of  freedom  increases  implies  that  less  emphasis 
need  be  placed  on  obtaining  noise  statistics  or  conversely  more  emphasis  should  be  placed 
on  signal  statistics.  This  comment  is  more  appropriate  for  short  range  detections,  less 
appropriate  for  long  range  detections  (especially  with  high  resolution  spatial  sensors). 


• Predicted  performance  improvements  based  on  nonfluctuating  R can  lead  to  gross 

ovpr  pfitimatps 
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NUMERICAL  MODELS  OF  ACOUSTIC  PROPAGATION 
THROUGH  INTERNAL  WAVES 


H.  A.  DeFerrari  and  R.  Leung 
University  of  Miami 

Rosenthal  School  of  Marine  and  Atmospheric  Science 
Miami,  Florida  33149 


A theory  of  propagation  through  a fluctuating  ocean  (FDMZ  1977)*  predicts  statistics 
of  acoustic  variability  in  regions  of  A - 4»  space,  where  A and  are  the  diffraction  and 
strength  parameters.  Statistics  for  a single  path  or  for  large  numbers  of  micro-  or  deter- 
ministic multipaths  are  derived.  Applications  to  ocean  propagation  usually  employ  the 
Garrett-Munk  internal  wave  model  (GM72)t  to  predict  the  sound-speed  fluctuations  which 
are,  in  turn,  used  as  input  to  the  propagation  model.  Once  the  internal  wave  statistics 
are  specified,  A - 4>  parameters  are  computed  by  integration  along  the  path  of  propaga- 
tion which  requires  a detailed  ray  analysis. 

We  have  developed  general  numerical  methods  for  the  computation  of  A and  <t>.  For 
arbitrary  depth  dependence  of  sound-speed  C and  buoyancy  frequency  N,  a numerical 
model  computes  internal  wave  modes  subject  to  the  boundary  condition  of  vanishing  verti- 
cal velocity  at  the  surface  and  bottom.  The  approximate  forms  of  the  wave  functions 
used  in  GM72  give  good  approximation  in  the  interior,  away  from  the  boundaries,  and 
for  the  lower  and  more  energetic  frequencies  of  the  internal  wave  spectrum.  Our  numer- 
ically computed  wave  functions,  dispersion  relationships  and  spectra  are  valid  near  the 
surface  and  for  the  higher  frequencies  of  the  internal  wave  spectrum,  although  they  are 
considerably  more  complicated  to  compute.  Thus  we  anticipate  better  approximation  for 
ray  paths  which  have  turning  points  near  the  surface. 

Ray  paths  are  determined  with  conventional  ray  tracing  equations  by  fitting  the 
sound-speed  profile  with  a large  number  of  linear  segments.  The  ray  model  is  used  to 
approximate  the  path  of  propagation  only,  and  therefore,  errors  are  not  introduced  by 
the  discontinuous  derivatives  at  the  interception  of  the  linear  segments.  With  the  sound- 
speed  fluctuations  determined  from  the  internal  wave  model  and  the  ray  path  known,  4» 
is  determined  by  numerically  evaluating  the  integral. 

4>2  = q§  J dx  j'dx'p[z(x),  z(x')]  . 

Computations  of  A are  tedious,  requiring  the  determination  of  Rf  the  Fresnel  zone 
radius  along  the  propagation  path.  For  the  calculation  of  Rf  we  have  used  the  ray  tube 
method  (FDMZ  1977).  To  employ  the  method,  it  was  necessary  to  solve  the  differential 
equation. 

*Flatte,  Dashen,  Munk,  Zachariasen  Model  1977 
tGarrett  Munk  Model  1972. 
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9xx£(x)  + Uo[zray(x)]J(x)  = 0 

where,  £(x)  is  a small  displacement  from  the  unperturbed  ray  and  U0  = (c2  - c§)/2c2. 


r 


Measured  values  of  Uq(z)  are  fit  with  a cubic  spline,  which  insures  unique  and  con- 
tinuous second  derivatives  at  every  depth.  Finite  difference  methods  are  used  to  solve 
the  equation  for  any  specified  ray  path.  A is  then  computed  by  numerical  evaluation  of 
the  integral. 

A = o'2 q$  J dx  J~dx’p[z(x),  z(x')]  (g0AL2) 
where  Lv  is  the  vertical  coherence  length. 

The  models  described  have  recently  been  completed  and  we  have  begun  to  system- 
atically examine  the  effects  of  source-receiver  geometry,  sound-speed  profile,  and  Brunt- 
Vaisala  profile.  Initial  results  indicate  that  A parameter  is  very  sensitive  to  ray  geometry. 
The  Fresnel  zone  radius  changes  greatly  along  the  propagation  path,  having  large  contribu- 
tions to  the  A integral  near  convergences  and  caustics,  and  little  near  the  turning  point. 
Source-receiver  depth  can  result  in  order  of  magnitude  changes  even  at  the  same  source- 
receiver  range.  Rays  having  turning  points  near  the  surface  incur  significantly  less  <1>  varia- 
tion than  predicted  with  (GM72).  Examples  are  presented  to  illustrate  these  points. 

Our  objectives  are  to  extend  results  for  the  single-ray  case  to  the  special  cases  of 
two  and  four  paths,  which  are  typical  in  long-range  deep-ocean  propagation  and  result 
from  combinations  of  up-going  and  down-going  rays.  Such  paths  very  often  have  nearly 
identical  travel  times  and  cannot  be  separated  by  temporal  processing.  Direct  observation 
of  statistics  of  the  combined  paths  is  easily  accomplished  in  ocean  experiments,  whereas, 
isolation  of  a single  path  is  often  impossible  for  frequencies  below  500  Hz  at  ranges  in 
excess  of  100  km. 


SINGLE  PATH-PHASE  AND  AMPLITUDE  FLUCTUATIONS 


Terry  E.  Ewart 
Applied  Physics  Laboratory 
University  of  Washington 
Seattle,  Washington  98X95 


ABSTRACT 

A discussion  of  this  signal  processing  techniques  and  the  results  from  two 
acoustics  transmission  experiments  will  be  presented  (the  range  of  acoustic  fre- 
quencies is  2 to  16  KHz).  The  focus  of  these  experiments  is  to  measure  the 
stochastic  behavior  of  the  phase  and  amplitude  of  a single  Fermat  path,  to  meas- 
ure the  index  of  refraction  fluctuations  in  time  and  space  over  the  transmission 
path  and  compare  the  results  with  theoretical  predictions  of  acoustic  scattering 
based  on  the  index  of  refraction  fluctuation  measurements.  We  have  concen- 
trated on  the  range  of  time  and  space  scales  appropriate  to  the  tides,  internal 
waves  and  fire  structure.  We  show  that  the  signal  processing  technique  used  to 
separate  single  paths  from  a multipath  environment  work  to  the  limit  of  one 
wave  length  path  separation  (which  is  the  limit  of  distinct  Format  path  separa- 
tion). Acoustic  frequency  dependent  effects  are  observed  in  the  amplitude  fluc- 
tuations. These  effects  and  others  are  discussed  in  the  context  of  the  environ- 
mental observations.  The  results  from  the  second  experiment,  MATE  (Mid- 
Ocean  Acoustic  Transmission  Experiment),  are  preliminary,  but  serve  to  indicate 
the  focus  of  the  research  effort. 


PREDICTION  OF  DETECTION  PERFORMANCE 


Magnus  Moll 

Bolt  Beranek  and  Newman  Inc. 
1 701  North  Fort  Myer  Drive 
Arlington,  Virginia  22209 


ABSTRACT 


Statistical  methods  of  systems  analysis  are  employed  to  predict  the  per- 
formance of  a passive  sonar  receiver  operating  in  environments  with  fluctuating 
signals  and  noise.  The  receiver  analog  employed  is  a multichannel  processor 
that  produces  for  each  channel  a decision  threshold  that  is  determined  by  the 
outputs  of  other  channels.  The  receiver  inputs  are  assumed  to  be  a compound 
random  process  exhibiting  characteristics  observed  in  experimental  studies  of 
ambient  noise.  The  parameters  of  these  processes  including  their  relaxation 
times  can  be  selected  to  suit  a particular  set  of  environmental  conditions.  The 
principal  objective  is  to  achieve  realistic  predictions  of  the  probability  of  de- 
tection by  single  observations.  Results  obtained  to  show  that  receiver  perform- 
ance depends  on  both  the  time  scale  of  the  fluctuation  and  its  extent.  The 
effort  is  supported  by  the  Naval  Analysis  Programs,  Code  431,  Office  of  Naval 
Research. 


SUMMARY 

Statistical  methods  of  systems  analysis  are  employed  to  predict  the  performance  of  a 
passive  sonar  receiver  operating  in  environments  with  fluctuating  acoustic  inputs.  The  use 
of  this  approach  requires  an  analog  of  the  sonar  receiver  and  statistical  characterization  of 
acoustic  signals  and  noise. 

A multi-channel  analog  was  employed  to  represent  either  a multi-beam  receiver  and 
recording  display  or  a frequency  analysis  device.  Each  channel  includes  a rectifier  and  an 
integrator.  For  detection  decisions,  the  output  of  each  integrator  is  compared  to  a thresh- 
old value  that  is  a linear  combination  of  the  outputs  of  the  integrators  of  other  channels. 
These  thresholds  therefore  adjust  to  variations  in  noise  background  level. 

The  random  process  model  employed  has  characteristics  that  conform  to  the  findings 
of  several  experimental  studies  of  underwater  ambient  noise: 

1.  Estimates  of  first-order  probability  measures  from  short  data  segments  supported, 
more  often  than  not,  the  hypothesis  that  the  amplitude  of  noise  is  normally  distributed. 

2.  However,  the  number  of  departures  from  first-order  normality  increased  with 
sample  length. 

3.  In  one  set  of  experiments,  about  28%  of  the  estimates  of  kurtosis  fell  outside  the 
95%  confidence  limits  for  a Gaussian  variate,  and  almost  all  of  these  were  on  the  high  side. 
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These  traits  would  be  exhibited  by  a class  of  compound  random  processes  represented 


N(t)  = JW)  G(t)  (1) 

where 

P(t)  is  a non-negative  stationary  random  process 

G(t)  is  a zero-mean  unit-variance  stationary  Gaussian  process,  statistically  independ- 
ent of  P(t). 

It  is  also  assumed  that  the  relaxation  time  of  P(t)  is  much  greater  than  that  of  G(t). 

If  the  last  assumption  is  satisfied,  then  the  bandwidth  of  the  process  is  essentially  deter- 
mined by  G(t).  It  is  clear  from  (1)  that  given  P(t),  the  process  N(t)  is  conditionally  a 
nonstationary  Gaussian  process.  It  is  easily  shown  that  all  odd-order  moments  of  N (t) 
are  zero. 


Since  each  channel  of  the  receiver  analog  squares  its  input,  the  properties  of 

N2(t)  = P(t)G2(t) 


(2) 


affect  the  performance  of  the  post-rectification  smoother.  Analysis  shows  that  the  relaxa- 
tion time  of  N2(t)  is  effectively  determined  by  P(t),  which  can  be  considered  as  a fluctu- 
ating power  envelope.  Starting  with  (2),  it  can  be  shown  that  the  kurtosis  of  N (t)  is 


where 


**4N 

°N 


H4  N is  the  fourth  central  moment  of  N(t) 
is  the  variance  of  N(t) 
o£  is  the  variance  of  P(t) 
mP  is  the  average  value  of  P(t). 

The  kurtosis  for  a normal  variate  is  3. 


(3) 


For  the  receiver  channel  of  interest,  the  test  statistic  that  provides  the  basis  for  detec- 
tion decisions  is  the  difference  between  the  output  of  its  integrator  and  the  threshold  func- 
tion. For  channel  0,  the  test  statistic  can  be  expressed  as 


Z0  = T-l 


tS0(u)  + N0(u)V 


E c0iNf(u), 


du 


(4) 


where 


! ’ 
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S0(t)  is  the  signal  input  to  channel  0 
Nj(t)  is  the  noise  input  to  channel  i,  i = 0,  1,  2,  ....  n 
c0i  is  the  weighting  coefficient  for  channel  i in  forming  the  threshold  for  channel  0. 

If  the  output  of  the  integrator  of  channel  0 exceeds  its  threshold  value,  then  Z0  > 0.  Thus 
the  probability  of  detection  is 


PDo  = P[Z0  > 0|S0(t)  > 0]  = / dzfZo[z0|S0(t)>OJ  (5) 

Jo 

where 

fZo(z0l)  is  the  conditional  probability  density  function  for  the  test  statistic  Z0. 

The  coefficients  c0j  were  selected  to  maximize  detection  performance  in  isotropic  noise 
and  to  achieve  a specified  probability  of  detection  when  the  noise  is  Gaussian. 

The  principal  burden  of  the  analysis  is  to  derive  the  probability  density  function  for 
the  test  statistic.  This  is  a difficult  problem  in  that  the  inputs  to  the  averagers  are  non- 
Gaussian.  The  approach  chosen  was  to  approximate  the  required  density  by 

m 

h(x)  = g(x)  ®ix'  (6) 

i = 0 

where 

g(x)  is  a probability  density  function  selected  as  the  basis  of  the  expansion. 

The  coefficients  a;  are  selected  so  that  the  first  m moments  of  h(x)  match  those  of  the  test 
statistic  Z0. 

The  use  of  the  method  therefore  requires  the  derivation  of  low-order  moments  of  the 
test  statistic  Z0.  Obtaining  the  first  moment  is  trivial,  the  second  moment  interesting, 
and  the  third  moment  extremely  burdensome.  The  analysis  show  that  the  calculation  of 
a given  order  of  moment  requires  knowledge  of  the  corresponding  order  of  moment  func- 
tions of  the  power  envelope  processes.  Expressions  for  moments  one  through  three  have 
been  derived. 

A special  case  with  the  following  features  will  be  discussed: 

• All  channels  have  identical  parameters. 

• All  channels  have  the  same  power  envelope  function  PN(t)  for  noise. 

• A large  number  of  averager  outputs  are  combined  to  form  the  detection  threshold 
function. 

• The  signal  is  Gaussian. 
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m 


For  this  case,  the  mean  value  of  the  test  statistic  is  approximately 


mZ  - Ps  ' 


v/WT 


where 


ps  is  average  signal  power 

pN  = E[PN(t)]  = average  noise  power 

7 is  a parameter  whose  value  is  determined  by  the  false  alarm  probability  for  the 
case  of  Gaussian  noise 

W is  input  bandwidth 

T is  post-rectification  averaging  time. 

When  signal  is  absent,  mz  is  obviously  negative. 

For  the  same  case,  the  variance  of  the  test  statistic  is  approximately 

a\  = (WT)-1[(ps  + pN)2  + ffp(l  + 72Ip)]  (8) 

where  a p is  the  variance  of  the  noise  envelope  process 


rT 

IP  = T_1  / dt(l-  |t|)kp(t) 
J .fT 


where 


kp(t)  is  the  autocovariance  coefficient  function  for  the  power  envelope  process. 

The  second  term  of  (8)  shows  the  effects  of  the  power  envelope  fluctuation.  The  contri- 
bution is  proportional  to  the  variance  of  the  envelope  process.  Expression  (9)  shows  that 
the  variance  of  the  test  statistic  is  affected  by  the  relationship  between  the  post-rectification 
averaging  time  and  the  relaxation  time  of  the  power  envelope  process.  If  the  latter  is  much 
larger  than  the  former,  then  I_  ->•  1,  the  maximum  value.  This  assumption  is  implicit  in  the 
usual  analysis  of  the  effects  of  fluctuation. 

An  expression  for  the  third  moment  has  also  been  derived,  and  all  three  moments  will 
be  employed  to  determine  the  coefficients  of  (6)  for  m = 3.  That  result  in  turn  will  be 
employed  in  (5)  to  examine  the  dependence  of  detection  performance  on  the  statistical 
parameters  of  the  power  fluctuation  envelope. 

For  some  cases  of  interest,  a normal  probability  density  function  can  be  used  as  the 
basis  function  g(x)  in  (6).  For  a third-order  approximation,  the  probability  of  detection 
can  then  be  expressed  as 


I 


where 


P(  ) is  the  normal  distribution  function  (0,  1) 

Z(2)(  ) is  the  second  derivative  of  the  normal  density  function  (0,  1) 

M3Z  is  the  third  central  moment  of  the  test  statistic  Z0 . 

The  form  (10)  is  convenient  in  that  tables  of  both  of  the  required  functions  are  readily 
available. 

The  result  (10)  has  been  applied  to  a case  in  which  (op  t mP)2  = 1/2.  Substitution 
of  that  value  in  (3)  yields  a kurtosis  of  4.5  for  the  noise  process.  It  was  further  assumed 
that  the  autocovariance  function  for  the  power  envelope  is  exponential.  The  results  of 
the  calculation  are  plotted  as  a set  of  detection  transition  curves,  which  show  that  the 
spread  of  the  curves  increase  with  the  ratio  D/T,  where  D is  the  relaxation  time  of  the 
envelope  process,  and  T is  the  post-rectification  averaging  time  (See  Fig.  1). 

Mean  and  variance  of  the  test  statistic  are  given  by  (7)  and  (8)  respectively.  The 
former  requires  the  mean  value  of  the  power  envelope,  and  the  latter  additionally  requires 
its  autocovariance  function.  Such  data  have  been  acquired  experimentally.  Calculation 
of  the  third  moment  of  the  test  statistic  requires  a third  moment  function  for  the  power 
envelope,  a function  that  is  not  commonly  estimated.  We  are  investigating  the  possibility 
of  approximating  this  function  on  the  basis  of  the  third  moment  and  autocovariance  of 
the  power  envelope.  If  that  effort  is  successful,  then  more  realistic  performance  predic- 
tions can  be  based  on  a modest  amount  of  experimental  data. 


probability  of  detection 


D:  envelope  relaxation  time 
T:  averaging  time 
noise  kurtosis:  4.5 


Vo  4T  = 0.05 
no  envelope  fluctuation 


signal  to  noise  level  minus  recognition  differential  for  normal  case,  dB 


Transition  curves 


ACOUSTIC  FLUCTUATION  MODELING  FOR  SYSTEM 
PERFORMANCE  ESTIMATION 
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SUMMARY 

A critical  part  of  the  analysis  and  estimation  of  sonar  system  performance  is  the  char- 
acterization of  the  acoustic  signal  and  background  noise.  Current  techniques  in  perform- 
ance and  engagement  modeling  employ  random  processes  to  simulate  the  long  and  short- 
term variations  of  these  parameters  in  time.  The  usual  justification  is  that  the  complex 
mechanisms  govering  the  generation,  transmission,  reception,  and  processor  response  of  sig- 
nals and  noise  are  sufficiently  “random”  to  be  modeled  by  stochastic  processes  determined 
by  a few  input  parameters.  Since  such  an  approach  is  much  more  efficient  than  measuring 
or  modeling  the  details  of  the  acoustic  environment,  it  is  widely  utilized.  However,  it  has 
not  undergone  a comprehensive  evaluation. 

SAI  has  performed  an  investigation  into  several  aspects  of  the  validity  of  the  random- 
process  approach.  Among  the  topics  considered  are: 

(a)  the  basic  statistical  properties  of  acoustic  signal  and  noise  fluctuations 

(b)  the  ability  of  the  random-process  models  to  properly  simulate  the  important 
properties  of  the  acoustic  variables 

(c)  the  ability  of  the  stochastic  models  to  accurately  predict  measures  of  system  ef- 
fectiveness 

(d)  the  choice  of  a stochastic  model  and  estimation  of  its  inputs,  for  a particular 
scenario. 

The  ideal  way  to  address  these  would  be  to  compare  stochastic  simulations  with  measured 
acoustic  and  performance  data  for  many  ocean  environments,  sonar  systems,  and  engage- 
ment geometries.  However,  the  scope  of  the  study  and  the  availability  of  appropriate 
measurement  data  required  the  investigation  to  focus  on  a simplified  test  case  for  a repre- 
sentative acoustic  environment,  sonar  system,  and  target. 

A single  ocean  environment  (in  the  Northeast  Pacific  during  Summer)  and  target 
scenario  have  been  selected  for  the  test  case.  To  ensure  that  the  effects  of  the  environ- 
mental acoustic  parameters  could  be  gauged,  the  study  assumed  a stable  narrowband 
source  (at  25  Hz),  an  idealized  towed  array  receiver,  and  a processing  system  which  auto- 
matically detects  when  the  beam-output  signal-to-noise  ratio  (SNR)  exceeds  a threshold. 
Evaluation  of  the  random-process  approach  under  these  conditions  then  requires  the  de- 
tailed space/time  properties  of  the  transmission  loss  (TL)  and  ambient  noise.  A novel 
aspect  of  the  present  investigation  is  the  use  of  state-of-the-art  acoustic  propagation  and 
ambient-noise  models  to  fabricate  such  a controlled,  representative  acoustic  environment. 


For  the  test  scenario  and  target  conditions,  there  is  good  reason  to  attribute  fluctua- 
tions in  TL  to  relative  receiver  movement  through  the  source’s  multipath  interference  field 
(assuming  little  influence  from  ocean  dynamics).  The  acoustic  TL  “data”  were  generated 
from  models  which  fully  account  for  this  mechanism.  Similarly,  a special  ambient  noise 
model  was  developed  to  simulate  the  dominant  noise  fluctuation  mechanism:  movement  of 
many  contributing  surface  ships  through  the  interference  field. 

Once  signal  and  noise  time  series  were  constructed  with  the  acoustic  models,  several  of 
the  more  popular  random-process  models  were  used  to  simulate  the  series.  Input  param- 
eters (the  trend  in  TL,  the  variance  and  decorrelation  scales  of  TL  and  noise)  were  deter- 
mined directly  from  the  acoustic  data.  Comparisons  were  then  made  of  the  statistical  prop- 
erties of  the  acoustically  and  stochastically  modeled  signal  and  noise,  as  well  as  of  the 
detection  properties  of  the  signal-to-noise  ratio. 

The  principal  conclusion  of  the  study,  subject  to  the  limitations  in  the  test  environ- 
ment and  data,  is  stated  generally  as: 

• Given  proper  input  parameters,  the  stochastic-process  approach  can  produce  ade- 
quate simulations  of  signal  and  noise.  The  accuracy  is  sufficient  to  provide  useful 
estimates  of  detection  probabilities  and  detection  histories. 

This  result  is  tempered  by  the  observations  that: 

• The  accuracy  and  utility  of  the  approach  are  greatly  dependent  on  the  availability 
of  accurate  input  data  (e.g.,  the  variance  of  beam-noise  fluctuations).  Such  data 
may  be  difficult  to  obtain,  and  without  them  the  predictions  will  usually  be  poor. 

• The  accuracy  also  depends  on  the  particular  random  process  used.  This  study  con- 
centrated on  the  popular  ones,  and  the  agreement  with  acoustic  data  was  just  ade- 
quate. However,  significant  improvements  can  be  realized  when  the  basic  statistical 
properties  (e.g.,  marginal  density)  are  tailored  to  the  variable  and  case  at  hand. 

A report*  presents  the  details  of  the  comparisons  on  which  these  judgements  are  based,  and 
discusses  the  estimation  of  input  data,  the  choice  of  a particular  random-process  model,  and 
the  prediction  of  measures  of  effectiveness. 

The  work  was  supported  by  the  Naval  Analysis  Program  of  the  Office  of  Naval 
Research  (ONR  431). 

EXPLANATORY  NOTES  FOR  THE  FIGURES,  1 THROUGH  24 

1.  Concentrate  on  the  fluctuation  submodel  of  a performance  model.  Usually  a 
structured,  random-process  model  depending  on  2 or  3 input  parameters,  it  provides  a time 
series  of  fluctuations  which  in  turn  is  added  to  the  signal  or  noise  or  signal-excess  time 
series. 


*Cananagh,  R.C.,  “Acoustic  Fluctuation  Modeling  and  System  Performance  Estimation,”  Final  Report, 
Contract  N00014-76-C-0753,  Science  Applications,  Inc.;  McLean,  VA,  January  1978. 


97 


i 


* 


2.  A proper  test  of  the  random-process  approach  involves  comparison  of  acoustic 
data  with  stochastic  simulations.  Agreement  is  gauged  in  terms  of  detection-related  prop- 
erties (e.g.,  level-crossings)  of  the  time  series. 

3.  Ideally,  the  test  would  involve  comparisons  of  measured  data  with  simulated  data 
for  many  environments,  sensors  and  conditions.  The  approach  taken  in  the  study  was  to 
focus  on  one  environment  and  system.  Also,  state-of-the-art  acoustic-model  predictions 
were  used  in  place  of  measured  data. 

4.  & 5.  Four  questions  were  addressed  in  the  study,  and  four  general  conclusions  drawn. 

6.  A data  base  of  acoustically  modeled  data  was  generated.  Note  that  the  Parabolic 
Equation  model  and  a coherent  version  of  the  FACT  model  were  used  to  construct  complex 
pressure  fields  for  waterborne  and  bottom-bounce  paths  respectively  ; the  two  fields  were 
summed  in  phase.  A detailed  simulation  of  shipping  noise  was  constructed  with  the  DSBN 
noise  model.  In  each  case,  the  ocean  medium  is  static,  so  that  all  fluctuations  are  due  to 
source/receiver  motion  through  the  detailed  multipath  interference  field. 


7.  Three  target  tracks  were  considered,  two 
beam  and  one  radially  closing. 


‘tangential”  to  the  array’s  broadside 


8.  The  “acoustic’ 
‘stochastic”  data. 


data  base  is  shown.  It  serves  as  the  control  data  for  the  test  of  the 


9.  Three  random-process  models  were  considered,  and  their  properties  are 
summarized. 

10.  The  random-process  models  provide  an  additive,  dB  fluctuator  X(t)  for  the  vari- 
able Z(t),  be  it  signal  or  noise  or  signal-to-noise  ratio.  Note  that  signal  fluctuations  are  de- 
rived from  transmission  loss  (TL)  and  range  rate. 

11.  Although  the  picture  does  not  have  proper  resolution,  the  acoustically -modeled 
TL  fluctuations  in  range  are  compared  for  various  scales  of  trend  removal.  The  “trend” 
might  be  thought  of  as  the  deterministic  part,  presumably  a quantity  which  can  be  pre- 
dicted with  confidence. 


12.&13. 

shown. 


Some  statistical  properties  of  the  acoustically-modeled  TL  fluctuations  are 


14.&  15.  Compare  properties  inherent  to  the  random-process  models  with  those  of  13.  and 
14.  Notable  discrepancies  include  the  symmetry  properties  of  the  distributions  of  TL  and 
the  form  of  the  autocovariance.  Best  agreement  was  for  the  Jump  process  added  to  the 
2-nm  intensity -averaged  TL. 

16.  Acoustically  modeled,  omnidirectional  ambient  noise  had  properties  consistent 
with  those  inherent  to  the  random  processes. 

17.  Beam  noise,  on  the  other  hand,  can  be  a complicated  function  of  steering  angle 
•uid  beam  width,  as  illustrated. 


18.  Some  properties  of  the  acoustically  modeled  beam-noise  fluctuations  are  shown. 
Comparison  with  stochastic  data  showed  reasonable  good  agreement  when  the  Jump  model 
was  used  with  accurate  inputs. 

19.  Finally,  consider  SNR  and  a detector  which  detects  when  SNR  exceeds  a thresh- 
old (the  case  when  there  is  much  incoherent  averaging  of  filter  outputs).  For  the  special 
case  at  hand,  fluctuations  are  driven  by: 

(a)  Signal  fluctuations  caused  by  changes  in  target  range  (R(t)) 

(b)  Signal  fluctuation  caused  by  changes  in  target  bearing  (0(t))  which  affect 
array  signal  gain 

(c)  Noise  fluctuations  caused  by  ship  movements  in  R and  0. 

20.  Acoustically  simulated  SNR  for  the  tangential  target  (Track  1)  as  seen  on  two 
neighboring  beams  is  shown.  The  variations  are  dominated  by  noise  and  the  signal  gain. 

21.  For  the  same  noise  sources  as  20,  the  acoustically  modeled  SNR  for  the  radial 
target  (track  3)  is  shown  for  both  short  and  long  range  cases.  In  the  former,  the  rapid 
bottom-bounce  variations  are  apparent. 

22.,  23.  & 24.  Detection-related  statistics  for  the  acoustic  and  stochastic  data  are 
compared.  The  best  stochastic  models  and  perfectly  accurate  inputs  were  used.  The 
judgement  of  “adequate”  agreement  is  based  on  such  comparisons  as  these. 
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Figure  1 
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TEST  OF  THE  RANDOM-PROCESS  APPROACH 


GIVEN  REAL  SIGNAL  AND  NOISE  FLUCTUATION  DATA:  ACOUSTIC  DATA 


EXTRACT  INPUT  PARAMETERS  FOR  RANDOM-PROCESS  MODELS 


fi  USE  POPULAR  RANDOM-PROCESS  MODELS  TO  SIMULATE  SIGNAL  AND  NOISE 


time  series:  STOCHASTIC  DATA 


• question:  can  the  random-process  approach  adequately  simulate 


SIGNAL,  NOISE,  SNR? 


Figure  2 


APPROACH 


Ideal:  - 


MANY  ENVIRONMENTS,  GEOMETRIES,  SYSTEMS,  TARGET  PROPERTIES, 
FREQUENCIES,  ETC,  ETC. 


Practical:  - one  ocean  environment  (n.  pacific) 

- one  system  type  (towed  array) 

- IDEALIZED  DETECTOR  (SNR) 

- SINGLE  FREQUENCY  (25  Hz) 


- FEW  TARGET  TRACKS 


USE  OF  ACOUSTIC  MODELS  TO  PRODUCE  CHARACTERISTIC  SIGNAL  AND 
NOISE  TIME-SERIES  DATA 


Procedure : 


GENERATE  ACOUSTIC  DATA  FOR  TEST  CASES 

EXTRACT  FLUCTUATION  PROPERTIES  FOR  INPUT  TO  STOCHASTIC  MODEL 
SIMULATE  DATA  WITH  RANDOM-PROCESS  MODELS  TO  YIELD  STOCHASTIC  DATA 
COMPARE  BASIC  STATISTICS,  LEVEL-CORSSING  PROPERTIES,  DETECTION  HISTORIES 


Figure  3 
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QUESTIONS 


(1)  What  are  basic  statistical  properties  of  acoustic  signal  and  noise  data? 


(2)  Can  random-process  models  simulate  acoustic  data? 


(3)  Can 


APPROACH  ADEQUATELY  PREDICT  MOE's? 


(4)  HOW  CAN  PREDICTIVE  ACOUSTIC  MODELING  AID  RANDOM-PROCESS  APPROACH? 

Figure  4 

CONCLUSIONS 

(1)  Acoustic  data  have  expected  properties. 

(2)  Given  accurate  inputs,  random-process  models  give  adequate 

SIMULATION  OF  SIGNAL,  NOISE,  AND  SNR. 

(3)  Random-process  approach  can  be  improved. 

(4)  Prediction  of  input  data  is  biggest  problem. 

Figure  5 


• Array: 


• Signal: 


• Noise: 


ACOUSTIC  DATA 


TLCR) 

PE  + CFACT 

Simplified  environment 
1/6  NM  SAMPLING  RATE 

SHADED,  HORIZONTAL  LINE 

8°  AND  4°  BEAMWIDTHS,  STEERED  NEAR  BROADSIDE  AT  0°  T 
30  AND  « DB  SI DELOBE  SUPPRESSION 

PERFECTLY  STRAIGHT  AND  STATIONARY  WITH  IDEAL,  PLANE- 
WAVE  RESPONSE,  AG(*) 

SI  CONSTANT 

Constant  velocity  target  tracks:  (R(t)^(t))  with 
5 KNOT  SPEED 

S(T)  = Sl-TLCR(r))  + AG(*(t)) 

LIKE  SIGNAL  MODEL,  BUT  MANY  SHIP  SOURCES 
SHIPS  DISTRIBUTED  ON  LANES  ROUGHLY  EAST-WEST 
RANDOM  COURSE,  SPEED,  SOURCE  LEVEL,  ARRIVAL  TIMES 
INCOHERENT  SUM  OF  SOURCES 

Figure  6 
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Figure  7 


Transmission  loss  versus  Ranee 


- sampled  every  1/6  nm  from  0 to  500  nm 

- for  receiver  depth  of  1000  feet  and 
source  depths  of  60  and  300  feet. 

- at  frequency  25  Hz 

Signal  Level  versus  Time 


- sampled  every  2 minutes 

- for  target  tracks  1,  2,  and  3 

- for  source  depth  of  300  feet  and  re- 
ceiver depth  at  1000  feet 

- for  array  response  parameters  of  Table 
2-1 

Beam  Noise  versus  Time 


- 50  replications,  each  of  duration  10 
hours  and  one  replication  of  duration 
35  hours,  sampled  every  2 minutes 

- for  ship  parameters  described  in  2.2 

Siemal-to-Noise-Ratio  versus  Time 


- as  constructed  from  Signal  and  Noise 
time  series 

Ran ge-Averaged  TL 


- 2 nm  and  5 nm  intensity  averages  for 
above  cases 

- A+BlogR  fit  (LMS)  over  0-100  and  100-500  nm 

Figure  8 


USE  OF  RANDOM-PROCESS  MODELS 


Variable:  Z(t)  (dB's) 

Random  Process:  X(t) 

Simulation:  Z(t)  = Z(t)  + X(t) 

Z(t)  = Z 


TL(R)  = TL(R)  + X(R) 


TLCR(t) ) = TL(R(t) ) + X(R(t)) 

Figure  10 


longer  waiting  times  for  TL<L 

Figure  12 
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STOCHASTICALLY  MODELED  TL  FLUCTUATIONS 


- LOG-NORMAL  INTENSITY  DISTRIBUTION 


- SMOOTH  SPECTRUM 


- SYMMETRIC  LEVEL-CROSSING  PROPERTIES 


Figure  14 
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Figure  15 
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ACOUSTICALLY  MODELED  OMNI  NOISE 


- Number  of  ships:  about  150-160 

- Typical  o over  10  hours:  less  than  3 dB 

- Autocovariance:  exponent ial,  with  t between  0.5  and  2 hours 

- Skewness:  greater  than  1 

- Distribution:  when  not  too  many  nearby  ships,  chi-square  (high  d.f.) 

or  log-normal 

Figure  16 
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Beam  Noise  Time  Series  (Case  2) 


SIGNAL-TO-NOISE  RATIO  (SNR) 


SNR(t)  = SL  - TL(R(t))  + AG(*(t))  - N(t) 


track  1 


TRACK  3 


TRACK  2 


RECEIVER 


TARGET  TRACKS 


Acoustically  Simulated  SNR  for  Radial  Target  at  5 Knots  and  155  dB  Source  Level 
Array  Steering  * 90°  (0°T),  Target  Bearing  «*  0°T 
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larison  of  Acoustically  and  Stochastically  Modeled  Detection  Properties 
5-knot  Radial  Target  on  Broadside  Beam  of  Width  8°,  at  Ranges  100-50  m 

Figure  23 


arison  of  Acoustically  and  Stochastically  Modeled  Detection  Properties 
5-knot  Radial  Target  on  Broadside  Beam  of  Width  8°,  at  Range  300-250  m 
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BEAM  OUTPUT  FLUCTUATIONS  ON  TWO  TOWED  ARRAYS 


Andrew  G.  Fabula 
(See  Volume  2) 


FLUCTUATIONS  DUE  TO  RANGE  RATE 


1 


Ira  Dyer 

Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

(Editorial  Note:  A full  text  of  this  paper  has  not  been  received.  The  following  editorial 
summary  is  provided  to  give  guidance  to  the  reader.) 

The  purpose  of  this  research  is  to  determine  the  effects  of  platform  motion  on  the 
received  signal.  Assume  source/receiver  are  in  relative  motion  V and  a signal  at  single  fre- 
quency co0  is  transmitted.  Let  the  ocean  be  hypothesized  “frozen.”  It  is  then  required 
to  find  the  frequency  shift  wa  as  a function  of  cross-over  angle  0O,  and  velocity  V.  Evi- 
dently, 0O  depends  on  the  sound  speed  profile  of  the  channel.  Choose  a convenient  pro- 
file and  with  it  the  energy  E(0O)  of  the  ray  at  cross-over  angle  0O.  The  N normalized  auto- 
correlation 7(7-)  of  the  sinusoid  at  delay  r is  then  calculated.  From  this  one  derives  an 
expression  for  the  N normalized  power  spectrum  F(w ) of  the  fluctuation  of  the  received 
signal  due  to  range  rate  as  a function  of  the  energy  of  a ray  traversing  the  sound  speed  pro- 
files selected.  From  F(w ) one  can  obtain  mean,  variance  of  this  fluctuation.  Three  types 
of  channel  analyzed  are  (1)  isospeed,  (2)  bilinear,  (3)  Munk.  The  conclusion  is  that  the 
range  rate  spectrum  is  not  sensitive  to  static  background  model  of  the  sound  speed  profile. 
Next,  a comparison  is  made  between  fluctuation  due  to  range  rate  and  that  due  to  internal 
waves.  From  this  comparison  a critical  range-rate  is  found  which  provides  a frequency 
shift  just  equal  to  that  of  internal  waves. 
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COMPARISON  WITH  INTERNAL  WAVE  PHASE  RATE 


(1)  in  each  path,  for  internal  wave  perturbations: 

"rms  58  57  x 10"8  uo  v^km  . rad/sec 

(2)  over  all  paths,  for  frozen  ocean  range-rate: 

wrmg  = aft  = ak0|  V|  ^2/2  a < 1 

(3)  can  define  a critical  range-rate,  Vc,  from: 


y - 2C0  Vrms 

c a wo 


(4)  for  a = \p2  = 6 X 1(T2  (\[/  = 14°)  , R = 103  kr> 

V.  ~ 0.2  m/sec  ~ 0.4  knot 


THE  IMPORTANCE  OF  SOURCE  MOTION,  RECEIVER  ORIENTATION, 
AND  THE  OCEAN  ENVIRONMENT 


William  Jobst 

Institute  of  Acoustic  Research 
Miami,  Florida 


ABSTRACT 


(U)  When  a -ign”1  transmitted  from  a moving  source  is  observed  at  an  ar- 
ray of  sensors,  the  signal  field  fluctuates  as  functions  of  space,  time,  and  fre- 
quency. The  fluctuations  may  be  characterized  by  the  statistical  momenta  of 
the  received  signal  field-mean,  variance,  and  joint  second  order  moments 
typically  called  coherence.  We  examine  coherence  as  a function  of  space,  time 
and  frequency  using  data  from  recent  experiments  in  the  Atlantic  and  Pacific, 
is  shown  that  the  fluctuations  in  the  received  signal  field  are  related  to  source 
motion,  to  receiver  location  and  orientation,  and  to  the  ocean  environment. 


I 

II 
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IMPACT  OF  SOURCE  MOTIONAL  FLUCTUATIONS 
ON  INTERARRAY  SIGNAL  COHERENCE 

Albert  A.  Gerlach 
(See  Volume  2) 
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ACOUSTIC  FLUCTUATIONS 


The  Woods  Hole  Oceanographic  Institution 
and  the 

Institute  of  Geophysics  and  Planetary  Physics  of  the 
Scripps  institution  of  Oceanography 

Robert  C.  Spindel 

The  Woods  Hole  Oceanographic  Institution 


The  Woods  Hole  Oceanographic  Institution  and  the  Institute  of  Geophysics  and 
Planetary  Physics  of  the  Scripps  Institution  of  Oceanography  are  engaged  in  a cooperative 
research  program  to  study  the  effects  of  oceanic  variations  on  acoustic  propagation.  The 
objective  of  the  joint  program  is  to  identify  and  quantify  the  role  of  ocean  processes  (in- 
ternal wave  induced  temperature  and  current  fluctuations,  intrusive  effects  and  finestruc- 
ture,  currents,  mesoscale  eddies,  etc.)  in  determining  the  variability  of  sound  transmission. 
The  ultimate  goals  are  to  further  our  understanding  of  the  fundamental  limits  on  signal 
coherence  in  space  and  time  and  to  develop  techniques  for  monitoring  the  ocean  acous- 
tically. Scientists  at  both  institutions  have  been  active  in  acoustic  fluctuation  research  for 
a number  of  years.  Woods  Hole  has  concentrated  on  low  frequency  (100-400  Hz),  nar- 
rowband, long  range  (10-1500  km)  transmissions  whereas  IGPP  has  focused  on  high  fre- 
quency (2250  Hz),  wideband,  short  range  (25  km)  propagation.  The  former  studies  have 
emphasized  moored  sources  and  receivers  and  have  resulted  in  fluctuation  data  as  a func- 
tion of  range,  depth  and  frequency.  The  latter,  conducted  from  drifting  ships,  have 
yielded  amplitude  and  travel  time  data  along  selected  single  paths. 


THE  WOODS  HOLE  PROGRAM 

The  Woods  Hole  deep  ocean  experiments  have  used  primarily  moored  sources  and 
receivers  to  minimize  interactions  with  the  ocean  boundaries.  Acoustic  transmissions  have 
been  narrowband  so  that  acoustic  travel  time  variations  are  manifest  as  phase  (as  well  as 
amplitude)  fluctuations  in  the  received  signal.  A precision  Doppler  acoustic  position- 
keeping system  has  been  developed  to  permit  coherent  phase  measurements  over  ranges  of 
several  hundreds  of  kilometers.  Phase  and  amplitude  fluctuations  are  measured  using  re- 
mote digital  receivers  that  record  continuous  time  series  of  up  to  60  days  duration.  The 
receivers  contain  a digital  Doppler  tracking  device  for  0.1  meter  position-keeping.  The 
moored  acoustic  sources  have  a comparable  lifetime  and  transmit  continuous  tones  at  a 
level  of  176  dB  re  1 pPa. 

Relating  observed  fluctuations  to  ocean  processes  is  complicated  by  the  multipath 
field.  The  random  superposition  of  a large  number  of  multipaths  results  in  phase  (or 
travel  time)  fluctuations  that  obscure  those  due  to  ocean  variations.  Even  the  phase  of 
a single  (ray)  path  may  fluctuate  widly  because  the  sound  velocity  finestructure  can  split 
the  ray  into  a number  of  fluctuating  raylets  or  micromultipaths.  Thus  a single  transmitted 
pulse  can  be  received  as  a cluster  of  fluctuating  pulses  corresponding  to  totally  chaotic 
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interference.  The  particular  region  in  range-acoustic  frequency  space  where  this  occurs  is 
known  as  the  saturated  scattering  region  and  roughly  corresponds  to  long  ranges  and  high 
frequencies.  For  short  ranges  and  low  frequencies  scattering  is  negligible  and  phase  (travel 
time)  variations  are  relatable  to  mesoscale  ocean  processes.  A systematics  of  acoustic  fluc- 
tuation regimes  based  on  acoustic-internal  wave  interaction  has  been  developed  (Dashen, 
Flatte,  Munk,  Zachariasen)  to  predict  those  frequency,  time  and  range  conditions  under 
which  micromultipaths  are  expected  to  dominate  the  received  signal,  and  those  regimes 
where  phase  (travel  time)  measurements  can  be  related  directly  to  observable  mesoscale 
phenomena. 

The  Woods  Hole  effort  has  been  devoted  to  perfecting  techniques  to  permit  coherent 
phase  and  travel  time  measurements  to  be  made  from  deep  ocean  moored  systems.  The 
I analysis  of  Woods  Hole  continuous  tone  phase  and  phase  rate  fluctuation  data  has  under- 

scored the  important  role  that  internal  waves  play  in  acoustic  variability  with  periods  of 
several  minutes  to  a day.  Woods  Hole  experiments  have  been  accompanied  by  environ- 
mental ocean  measurements  such  as  temperature  and  current  in  order  to  permit  the  cor- 
relation of  acoustic  data  with  observed  ocean  variations.  In  addition,  WHOI  experiments 
have  been  designed  to  delineate  weak  and  strong  ocean  scattering  regimes  so  that  the 
parameters  of  an  ocean  acoustic  monitoring  system  can  be  defined. 


THE  IGPP  PROGRAM 

IGPP  has  performed  short  range,  high  frequency,  ship-to-ship  experiments  using 
broadband  pulse  compression  codes  to  measure  directly  acoustic  travel  time  fluctuations. 
Pulse  transmissions  have  allowed  observation  of  strong  micromultipath  regions  (near  sur- 
face paths)  and  weak  scattering  regions  (deep  paths).  Reciprocal  pulse  transmissions  have 
enabled  discrimination  between  currents  and  variable  finestructure.  IGPP  is  now  con- 
structing an  encapsulated  version  of  their  digital  processing  apparatus  so  that  the  2250  Hz 
transceivers  can  be  moored.  This  eliminates  travel  time  variations  due  to  ship  drift. 


SOME  NUMBERS 

It  is  the  goal  of  the  collaboration  between  IGPP  and  WHOI  to  demonstrate  the  fea- 
sibility of  monitoring  mesoscale  ocean  processes  via  acoustic  methods.  Direct  travel  time 
measurements  of  transmitted  pulses  (or  regenerated  pulses  from  wideband  transmissions) 
is  the  preferred  mode  of  operation  since  it  avoids  the  problem  of  multipath  interference 
and  resulting  phase  discontinuities  inherent  in  continuous  tone  systems.  A typical  meso- 
scale feature  has  a temperature  variation  of  about  2°C,  a fractional  change  in  sound  speed 
of  8c/c  5 X 10-3,  and  a spatial  extent  of  about  100  km.  This  large  temperature  anom- 
aly is  confined,  for  the  most  part,  to  the  upper  500  m.  Within  this  scale  there  will  be 
two  upper  turning  points  each  extending  for  about  20  km  within  the  upper  500  m.  The 
net  travel  time  variation  is  about 


T 


5c  x 2 X 20  km 


10-1  sec . 


Thus  we  require  a system  capable  of  0.1  sec  accuracy  in  travel  time  in  order  to  detect  the 
mesoscale  feature.  Multipaths  will  arrive  with  a separation  of  about  0.1  sec  (including 
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only  early  arrivals  and  excluding  the  closely  spaced  late  arrivals)  so  that  a time  ^so/ufion 
of  0.1  sec  is  required  to  separate  paths.  This  implies  a signal  bandwidth  of  10  Hz  or  more. 

Consider  now  the  errors  associated  with  the  required  travel  time  measurements.  At 
short  ranges,  or  low  frequencies,  micromultipathing  will  contribute  negligible  jitter  to  the 
received  signal.  Error  in  travel  time  measurement  arises  solely  through  additive  ambient 
noise.  The  rms  travel  time  error  due  to  noise  is  approximately 
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y/SNtl 


— • — — - 10' 


sec 


for  an  SNR  of  20  dB  and  a transmitted  band  of  10  Hz.  If  phase  coherent  processing  is 
used  this  error  can  be  reduced  to 


ot  - 


~ _L_  . - 1 — « 0.5  X 10~3  sec 

220 


y/SNR 


for  a 220  Hz  carrier  and  a SNR  of  20  dB.  In  either  case  it  is  clear  that  the  requisite 
accuracy  and  resolution  of  10' 1 sec  can  be  achieved  in  each  transmission. 

At  long  ranges  and  higher  frequencies  where  scattering  is  saturated,  micromultipaths 
may  cause  significant  travel  time  jitter.  Experimental  evidence  at  206  Hz  and  500  km 
(MIMI)  indicates  that  time  jitter  due  to  micromultipaths  will  be  small  (^10  sec)  imp  y g 
an  ocean  bandwidth  of  about  100  Hz.  Theoretical  formulations  (Flatte)  Predicta  so  e 
what  smaller  bandwidth,  -20  Hz,  but  still  sufficiently  wide  to  support  a 10  Hz  bandwidth 
transmission  and  a resolution  of  10' 1 sec.  If  we  take  the  lower  of  the  two  bandwidths 
resolution  will  be  of  order  5 X 10"2  sec.  This  can  be  improved  through  simple  averaging 
over  multiple  transmissions.  Since  decorrelation  times  are  of  the  order  of  200  sec,  both 
from  CW  measurements  (WHOI)  and  theory  (Flatte)  the  tune  bandwidth  product  of  the 
ocean  for  500  km  ranges  and  frequencies  of  several  hundred  Hz  is  about  2 X 10  Hz  e . 

A 10  day  experiment  with  a 10  Hz  transmitted  band  results  in 

10  days  X B6.4  X 104  sec/day  X 10  Hz  ^ 4Q() 

^ 2 X 104  Hz-sec 

independent  estimates  of  travel  time  so  that  the  achievable  accuracy  will  be 


Boo 


5 X 10" 3 sec 


These  simple  calculations  indicate  that  it  should  be  entirely  feasible  to  measure 
oceanic  mesoscale  processes  by  acoustic  methods. 


WHOI-IGPP  EXPERIMENTS 

A preliminary  test  of  these  lotions  is  to  be  conducted  during  the  spring  of  1978.  A 
moored  version  of  the  IGPP  computer  controlled  high  frequency  transceivers  will  be 
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deployed  at  25  km  spacing  together  with  WHOI  continuous  wave  and  10  Hz  bandwidth 
sources  and  receivers.  The  WHOI  10  Hz  band  source  has  a carrier  frequency  of  220  Hz 
modulated  by  a coded  sequence.  Reception  is  via  a microprocessor  controlled  moored 
receiving  capsule  at  ranges  of  25  and  125  km.  In  addition,  reception  is  scheduled  to 
occur  at  remote  hydrophones  some  1000  km  distant.  The  experiment  is  to  last  10  days. 
The  range  of  acoustic  frequencies  and  distances  spans  predicted  oceanic  scattering  regimes 
from  weak  to  full  saturation.  A rather  complete  environmental  measuring  program  in- 
cluding current  meter  and  CTD  stations  will  complement  the  acoustic  experiment. 

In  the  summer  of  1978  WHOI  will  conduct  another  experiment  in  conjunction  with 
an  intense  physical  oceanographic  experiment  known  as  the  Polymode  Local  Dynamics 
Array.  The  goal  of  the  oceanographic  experiment  is  to  monitor,  in  detail,  mesoscale 
processes  in  the  Atlantic.  The  WHOI  acoustic  experiment  consists  of  a continuous  tone 
source  placed  as  near  to  the  center  of  a mesoscale  eddy  as  possible,  with  three  receivers 
on  a 100  km  radius  from  the  source.  The  intent  is  to  measure  travel  time  anomalies, 
through  phase  variations,  of  a mesoscale  process. 


RANGE  INDEPENDENT  FLUCTUATIONS  AND  PATTERN  RECOGNITION 


OF  VERTICAL  ANGLE  OF  ARRIVAL  STRUCTURE 


Fred  H.  Fisher 
(See  Volume  2) 
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OMNIDIRECTIONAL  AMBIENT  NOISE  AS  A FUNCTION 
OF  DEPTH  AND  FREQUENCY  IN  THE  DEEP  NORTHEAST  PACIFIC 

J.  A.  Shooter 
(See  Volume  2) 


MEASUREMENT  OF  CHARACTERISTICS  OF  AN  ACOUSTIC  PROPAGATION 
CHANNEL  BY  INVERSE  FILTER 


Harry  A.  De Ferrari  and  Robert  F.  Tusting 
University  of  Miami 

Rosenthal  School  of  Marine  and  Atmospheric  Science 
Miami,  Florida  33149 

I 


Several  characteristics  of  a propagation  channel,  such  as  Doppler  spectrum,  scattering 
function,  or  pulse  response,  can  be  conveniently  and  rapidly  computed  in  real-time,  em- 
ploying a small  computer  and  the  Fast-Fourier  algorithm  by  the  method  of  inverse  filtering. 
We  have  used  the  method  with  good  success  in  several  experiments,  both  with  on-line  proc- 
essing and  post-processing  of  analog  tapes. 

We  report  results  and  examples  from  three  broad-band  experiments.  The  first  is  for 
transmission  over  60  km  between  a fixed  bottom-mounted  source  and  receiver  (750  Hz) 
carrier.  The  second  is  for  towed  source  transmission  at  a range  of  from  35  to  100  km  at 
speeds  of  from  3 to  8 knots.  And  the  third  is  for  transmission  from  a moored  source  to 
a moored  set  of  receivers  in  the  water  column  at  a range  of  90  km  and,  simultaneously, 

I to  a second  moored  hydrophone  at  a range  of  270  km. 

Inverse  filtering  consists  of  dividing  the  spectrum  of  the  received  signal  R(oj)  by  the 
spectrum  of  the  transmitted  signal  S(a>)  in  order  to  compute  the  system  response  of  H(gj) 
of  the  channel. 


H(co)  = R(cj)/S(gj)  , S(o>)  =£  0 . 

The  system  response  function  of  the  channel  can  then  be  transformed  over  the  band- 
width of  the  transmitted  signal  to  give  a pulse  response. 

+Aw 

p(t)  = / R(oj)/S(oj)  do)  . 

•'-Aw 

The  inverse  filter  is  known  to  perform  poorly  to  noise,  since  the  lower  values  of  the 
spectrum  are  contaminated.  However,  for  many  ocean  acoustic  measurements,  high  signal- 
to-noise  is  available  and  desirable  since  the  objective  is  to  measure  the  range  of  fluctuation. 

Augmented-linear-maximal-pseudo-random  sequences  (recommended  to  us  by  Metzger 
and  Birdsall,  University  of  Michigan)  have  the  necessary  nonzero  spectrum  for  the  trans- 
mitted signal.  The  sequences  have  length  N = 2n.  The  convenience  and  computational 
speed  of  the  method  is  achieved  by  storing  complex  demodulates  from  M successive 
sequences  to  form  a time  series  of  length  L. 


M = 2m 
L = 2m+n  . 

Processing  consists  of  first,  performing  a Fast-Fourier  transform  on  the  time  series  of 
length  L,  resulting  in  N X M spectral  line  (demodulates  are  also  stored  on  magnetic  tape 
for  post-processing).  Any  desired  characteristics  of  the  propagation  channel  can  then  be 
computed  from  spectral  lines. 

Denoting  the  index  of  each  spectral  line  by 

S; , i — 1,  2,  3,  L, 


we  note  that 

1.  i = 1 identifies  CW  carrier  line  assuming  0 Doppler. 

2.  The  broad-band  signal  strength  is  computed  by  averc  ring  the  spectral  lines 

1 N- 

N E S(pxm  + 1)  • 

P=0 

3.  The  pulse  response  at  zero  Doppler  is  obtained  by  performing  an  inverse  filter 
on  the  set  of  lines  given  by 

®pxm  + l » P - 1»  2,  ...,  N. 

4.  The  Doppler  shifted  pulse  response  may  reside  in  any  of  the  lines  separated  by 
an  interval  of  M lines. 

5.  The  scattering  function  is  computed  by  inverse  filtering  all  sets  of  lines  separated 
by  an  interval  of  M lines. 

Doppler  shift  is  readily  observable,  since  the  set  of  lines  containing  the  signal  energy 
can  be  observed  to  shift  to  another  set.  Doppler  compensation  can  be  accomplished  by 
adjusting  the  demodulating  clock  frequency.  In  experiments  where  the  source-motion  is 
measured  by  other  means,  such  as  bottom  navigation  systems,  the  effects  can  be  removed 
from  the  spectral  response  by  adjusting  each  of  the  spectral  lines  for  a time-varying  phase, 
computed  from  the  motion. 

Several  examples  of  time  history  of  pulse  responses,  scattering  functions,  Doppler 
compensation  and  the  removal  of  measured  source  and  receiver  motions  are  given  and 
discussed. 
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A REVIEW  OF  SIGNIFICANT  PAPERS  AT  94th  ASA  MEETING 

Sam  Hanish 

PREFACE 

From  the  very  beginning,  the  Invited  Papers  and  Evening  Workshop  on  Ocean  Fluc- 
tuations of  the  94th  Meeting  of  the  Acoustical  Society  held  14  December  1977  were  so 
j highly  specialized  that  any  brief  review  in  retrospect  of  the  major  papers  will  prove  of  little 

advantage,  except  to  the  deeply  entrenched  expert.  It  seems  more  appropriate  to  dwell  a bit 
on  the  key  issues  here,  in  order  to  bring  them  into  some  focus  and  thus  provide  a back- 
ground for  appreciation  of  Workshop  accomplishments.  However,  to  keep  the  main  ideas 
within  bounds  only  the  three  principal  invited  papers,  namely  by  Flatt^,  Dashen  and  Ellin- 
thorpe  will  receive  attention  here.  They  cover  a good  portion  of  fluctuation  modeling  and 
experiment  based  on  physics  now  current  in  this  country.  If  one  adds  to  them  the  work  of 
Beran  and  McCoy,  a rather  representative  total  of  the  fluctuation  modeling  in  the  ocean  be- 
comes available.  Although  McCoy  made  some  remarks  at  the  Workshop  the  effort  of  the 
Session  was  clearly  directed  toward  review  of  the  work  of  the  three  people  quoted.  We 
begin  then  with  some  general  remarks  on  the  subject. 

SCHOOLS  OF  RESEARCH  IN  FLUCTUATIONS  OF  ACOUSTIC 
SIGNALS  IN  THE  OCEAN 

Several  schools  of  research  are  currently  playing  a role  in  ocean  acoustic  fluctuations. 
Each  professes  to  have  a reasonably  complete  phenomenological  theory  of  the  physics. 
However,  upon  examination  it  is  found  that  the  completeness  of  the  theory  is  potential 
rather  than  proven.  In  actual  work  only  small  portions  of  the  theoretical  models  are  adapted 
to  numerical  calculation  and/or  experimental  verification.  The  distinction  between  schools 
therefore  rests  on  proven  parts  of  the  models  to  date  rather  than  on  their  potential  com- 
pleteness, which  is  about  the  same  for  all. 

The  first  school  is  that  of  the  JASON  Group  (Fatte,  Dashen,  Munk,  Dyson,  etc.). 
This  group  has  exercised  their  model  in  the  regime  of  point  temporal  statistics  of  fluctua- 
tions. This  means  they  have  used  their  model  to  predict  (or  explain)  the  temporal  character 
of  received  signals  at  a field  point,  or  temporal  coherence  between  field  points,  mostly  in 
the  region  of  saturation,  that  is,  at  ranges  where  the  fluctuating  part  of  the  received  CW 
signal  vector  executes  a random  walk  in  time.  They  consider  a sum  of  many  such  arrivals 
(=  multipath  arrivals)  at  a single  hydrophone,  and  couch  their  results  in  terms  of  the  natural 
log  of  the  amplitude  squared  and  in  terms  of  the  phase  of  the  multipath  sum.  In  experiment 
these  entities  are  time  series  and  their  power  spectra  can  be  modeled.  At  long  enough  ranges 
the  modeling  is  that  of  a random  walk.  In  a recent  paper  Dyson  et  al  JASA  59, 1121  (1976) 
have  compared  predicted  power  spectra  versus  frequency  with  project  MIMI  experimental 
data.  They  found  reasonable  agreement  at  frequencies  above  a parameter  7"1  = y/ o ^2, 
where  2 is  the  variance  of  the  fluctuating  time  derivative  of  phase,  but  at  very  low  fre- 
quencies the  power  spectra  actually  measured  differed  considerably  from  the  random  walk 
model. 


The  second  school  is  that  of  Beran  and  McCoy.  They  have  exercised  their  model  in 
connection  with  experiments  on  spatial  coherence  of  acoustic  wave-fronts  across  antenna 
apertures  of  many,  wavelengths  in  isotropic  and  anistropic  media.  This  means  they  com- 
pared the  statistics  of  phase  variation  across  a wavefront  as  a function  of  position,  aver- 
aging these  over  time.  They  found  good  agreement  with  theory  in  experiments  with  long 
horizontal  arrays  in  the  ocean  receiving  low  frequency  waves.  They  are  currently  working 
on  temporal  coherence  between  field  points  in  an  anistropic  medium. 

Other  schools  currently  working  in  ocean  fluctuations  in  the  USA  use  distinctly  dif- 
ferent approaches  than  the  two  noted  above.  These  are  (1)  those  that  use  the  Dyson  equa- 
tion and  the  Bethe-Salpeter  equation,  or  its  variants  in  the  form  of  perturbation  series  (2) 

* those  that  use  Monte  Carlo  methods  (3)  the  optical  school  who  use  the  concepts  of  optical 

propagation  through  a turbulent  medium  in  which  the  statistics  of  medium  randomness  are 
described  by  Kolmogoroff  spectra  of  velocity  turbulence  (4)  those  that  use  “phase  screen” 
methods  (or  path  integrals). 

As  noted  above  each  approach  can  be  formally  generalized  to  a single  all-embracing 
formalism.  A typical  example  is  the  formalism  of  the  path  integral  which  leads  to  almost 
all  known  results  by  symbolic  manipulation.  However,  the  all-inclusiveness  is  primarily 
potential,  actual  numerical  consequences  being  difficult  to  obtain  except  in  those  special 
cases  for  which  other  methods  have  supplied  the  needed  algorithms.  As  to  testing  models, 
agreement  with  experiment  over  all  ranges  of  parameters  is  not  achievable  with  known 
models  to  date.  What  is  at  fault  in  the  current  failure  to  make  a model  that  covers  every- 
thing is  the  complexity  of  the  ocean  environment.  This  requires  a special  note  which  is 
discussed  next. 

REGIMES  OF  FLUCTUATION  STATISTICS 

The  physical  parameters  of  ocean  fluctuation  are  so  diverse  and  cover  such  large 
range  of  values  that  theoretical  modeling  must  be  piecemeal  approximated  to  limited  por- 
tions of  the  total.  A first  problem  in  modeling  is  to  classify  the  regimes  in  which  the  sta- 
tistics of  models  are  thought  to  be  valid.  A convenient  classification  scheme  adopted  by 
the  JASON  Group  is  based  on  two  parameters,  4>  and  £2.  The  parameter  4>  measures  the 
severity  of  phase  randomization  of  a wavefront  caused  by  random. changes  of  sound  speed. 
The  second  parameter  ft  measures  interference  of  effects  of  a blob  of  inhomogenity  in  the 
path  of  a wavefront.  It  is  essentially  a statement  of  coherence  of  signals  radiated  from  the 
extreme  ends  of  a blob  and  crossing  each  other  at  some  distant  field  point.  A plot  of  ft  as 
ordinate  and  $ as  abscessa  provides  a 4-quadrant  classification  grid  in  <l>-ft  space  (Fig.  1). 

A typical  line  of  demarcation  is  4>-ft  which  divides  the  regime  of  saturation  (shown  shaded 
here)  from  all  other  regimes.  Indeed,  if  ft  /4>  « 1 the  coherence  (in  time)  is  neglible  rela- 
tive to  the  randomizing  effect  of  medium  inhomogeneities  taken  to  be  strong  in  this  case. 
The  remaining  quadrants  cover  regimes  of  partial  saturation  (Rytov  approximation)  and 
geometric  acoustics  (WKB,  Bom  etc.  approximations). 
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Figure  1 


Let  us  examine  the  regime  of  saturation  since  it  looms  so  large  in  the  following  dis- 
cussion of  the  ASA  Workshop.  The  term  saturation  is  best  understood  in  connection  with 
early  (circa  1965)  experiments  in  the  transmission  of  light  through  a random  medium.  A 
plausible  theory  suggests  that  the  appropriate  equation  in  the  scalar  electric  field  E is  V2E  + 
k2m2  E = 0.  A Russian  analyst  Rytov  (1937!)  suggested  solutions  in  the  form  E = exp(i//)  = 
AeiS,  = X + iS  where  x = In  A.  This  is  the  origin  of  the  famous  log  amplitude  fluctuation 
parameter,  and  its  physical  counterpart,  the  variance  o2  of  fluctuation  intensity  I, 


Later  Tatarski(1961, 1967)  used  a Kolmogoroff  spectrum  for  velocity  turbulence  to  develop 
an  explicit  formula, 

o2lni  = 1.23  C2n  k7'6  L11'6 


in  which  k = ^-  = -v1,  X is  the  wavelength  of  light,  L is  the  range  at  which  the  fluctuation 
of  light  intensity  I is  being  measured,  and  Cn  is  a constant  dependent  on  the  chosen  model. 
The  object  of  experiment  is  to  measure  o2^  by  measuring  <I2  > - <I>2  of  the  received 
light,  and  plotting  this  number  as  ordinate  against  the  theoretical  expected  value  as  abscissa. 
The  result  of  this  plotting  showed  that  for  small  ranges  L the  theoretical  formula  and  experi- 
ment were  linearly  related  (Fig.  2).  However,  when  the  range  became  great  enough  this 
linear  relation  failed.  In  fact  the  variance  of  log  amplitude  squared  no  longer  depended  on 
range.  It  was  then  realized  that  the  variance  had  become  saturated  (with  range),  and  the  sta- 
tistics of  the  field  resembled  that  of  a random  walk. 
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Experimental  a2jnI  versus  theoretical  ff2jni  of  Propagation 
of  Light  in  a Turbulent  Atmosphere 
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Saturation  of  the  variance  of  acoustic  intensity  is  an  experimentally  verified  fact,  as 
demonstrated  by  AFAR  (See  Ellinthorpe’s  conclusions  below)  and  MIMI  (Dyson’s  paper 
shows  this  to  be  true).  A good  portion  of  modeling-approaches  now  current  attempt  to 
explain  its  statistics. 

With  these  remarks  on  some  key  issues  we  could  go  directly  into  discussion  of  the 
ASA  papers.  But  it  will  be  still  valuable  to  try  to  outline  the  basic  mathematics  common  to 
all  the  models,  and  so  to  allow  an  appreciation  of  the  thread  that  ties  them  all  together.  We 
begin  by  noting  that  while  mathematical  models  of  the  fluctuation  of  acoustic  signals  in  the 
ocean  are  differentiated  one  from  the  other  by  what  restrictions  are  placed  on  the  source, 
on  the  medium  and  on  the  receiver,  a single  operator-type  equation  covers  them  all: 

Mp}  = Mp)  + A,  iQl  (1) 

in  which  p is  the  acoustic  pressure  at  the  receiver,  £b  is  the  background  operator  (that  is, 
it  describes  the  propagation  of  against  a non-fluctuating  background  medium),  £f  is  ficti- 
tious source  operator,  £q  a true  source  operator,  and  Q is  a true  source. 

The  meaning  of  this  equation  may  be  clarified  by  the  following.  In  most  models  the 
background  operator  £b  is  taken  to  be  a deterministic  function  of  space  and  time,  appear- 
ing as  a collection  of  differential  and  integral  forms,  plus  constants.  Spatial  variations  of 
the  medium  are  considered  small  over  a distance  of  a wavelength  of  the  acoustic  signal,  and 
temporal  variations  at  a particular  field  point  are  taken  to  be  slow  over  a period  at  the  fre- 
quency of  the  acoustic  signal.  A typical  example  is  the  Helmholtz  operator  describing  the 
propagation  of  sound  in  a homogenous  isotropic  medium  £b  = V2  + • Another  example 

is  the  parabolic  operator  which  appears  frequently  in  current  models  of  fluctuating  signals. 

In  contrast  to  £b  the  operator  £f  describing  fictitious  sources  is  taken  to  be  random 
(or  stochastic)  when  modeling  fluctuations.  Typical  fictitious  sources  are  signals  reflected 
from  rough  ocean  surfaces,  or  diffracted  by  inhomogeneous  blobs  in  the  ocean.  In  these 
cases  the  forms  of  £f  are  typically  algebraic,  with  space  and  time  dependency  that  are 
statistical. 

The  true  source  operator  £q  is  taken  to  be  deterministic  if  the  source  and  receiver 
are  fixed  in  space  and  communicate  with  each  other  with  deterministic  signals.  It  becomes 
a random  operator  if  the  source  and  receiver  are  in  relative  motion,  and  the  motion  is  ran- 
dom. Specific  forms  are  not  treated  in  the  following  discussions  since  the  ASA  meeting 
was  not  concerned  with  source/receiver  in  motion. 

The  operator  equation  is  father  to  a vast  literature.  The  number  and  variety  of  solu- 
tions testify  to  a multitude  of  applications.  A widely  used  solution  takes  the  form  of  the 
Green’s  function  G which  is  itself  a solution  of  the  differential  equation  in  the  operator  £b, 

£b  { G } = -6 

in  which  6 is  the  unit  amplitude  delta  function  souce.  By  inversion  of  £b  one  can  write 

G = £b'1{5} 
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A complete  knowledge  of  G allows  us  to  write  the  symbolic  solution  of  Eq.  (1) 


■BHH 


i 


p = GXf{p}  +G£q{Q}  (2) 

The  two  equations  (1),  (2)  are  the  customary  starting  points  for  modeling  acoustic  fluc- 
tuations in  the  ocean.  They  are  rewritten  in  a more  suggestive  form: 

(4  + £f)  {p}  = £f  {Q}  (3) 

(U  +GjCf){p}  = G£q  {Q}  (4) 

. 

in  which  i is  the  identity  operator.  Except  in  special  (trival)  cases  both  equations  must  be 
solved  by  approximation. 

It  is  to  be  appreciated  that  £f  and  £q  are  random  operators.  This  makes  the  solu- 
tions p to  be  random  variables  which  can  be  identified  only  through  their  moments  <p>, 
<p(l)p*(2)>,  <p(l)p*(2)p(3)p*(4)>  etc.  in  which  1,  2,  3,  4,  etc.  are  location  vectors  in 
the  field  of  p,  and  * signifies  conjugate. 

We  now  use  Eqs.  (1),  (2),  (3),  (4)  as  convenient  tools  to  review  the  principal  papers 
of  the  ASA  meeting. 

Paper  LI  Sound  transmission  through  a fluctuating  ocean.  S.M.  Flatt£  (UCLA). 


Flatt^  considers  a case  of  ocean  fluctuations  which  is  equivalent  to  Eq.  (3)  above.  He 
takes  Q = 0 (meaning  there  are  no  true  sources  in  the  volume  of  the  ocean  under  analysis) 
and  by  special  approximation  (=  parabolic  approximation)  derives  the  specific  forms, 

/ J_  _92_  ,_3_\ 

^ 2k0  3z2  * di  J 


£f  = k0  5c/c,  c = c(r,z,t),  k0  = — , Cj  = const. 

ci 

in  which  z is  the  depth  coordinate,  r is  the  range  coordinate,  5c  is  the  random  component 
of  the  sound  speed,  and  c is  the  sound  speed  (total)  taken  as  both  depth  and  range  depend- 
ent. The  field  'I'  (r,  z,  t)  is  thus  a solution  of 


1 32 

2k0  3z2 


* = ko 


3c 

c 


(5) 


The  fluctuations  in  sound  speed  (total)  c is  a slowly  varying  function  of  time.  Various 
time  scales  have  been  observed.  Of  great  interest  is  the  time  scale  1/3  h to  24  h.  The  pri- 
mary cause  of  random  fluctuations  in  c over  this  time  scale  is  considered  to  be  internal  waves. 
The  ratio  5 c/c  varies  with  space  and  time.  In  specific  models  of  fluctuations  the  calculation 
of  5 c/c  must  be  accurately  done  because  small  deviations  in  the  statistics  of  its  value  are 
important.  The  calculations  of  c is  taken  to  be  deterministic,  while  that  of  6 c by  definition, 
is  statistical. 
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In  the  absence  of  values  of  c experimentally  measured  everywhere  in  the  volume  of 
the  ocean  under  analysis  it  is  convenient  to  choose  a representative  background  sound  speed 
profile.  Flatt&  chooses  the  Munk  “canonical”  profile, 


C„  = c. 


{l  + e [e-’Ml  -r?)  ]}  , t?  = 2(z-za)/B 


Here  zA  is  the  depth  of  the  minimum  value  of  c,  labelled  CCp,  B is  an  estimate  of  the  depth 
thickness  over  which  Ccp  is  substantially  minimum,  cx  = const  = 1500  ms-1 , and  e measures 
the  order  of  magnitude  of  deviation  from  the  minimum  CCp.  A typical  choice  of  parameters 


is: 


zA  = 1000  m,  B = 1000  m,  e = 0.57  X 10"  2 


It  is  next  required  to  calculate  the  random  part  5c  caused  by  internal  waves.  In  this 
type  of  wave  motion  ocean  fluid  particles  move  in  elliptical  paths  in  the  vertical  plane  (as  do 
ocean-surface  waves)  with  velocity  magnitudes  that  are  functions  of  depth  in  the  ocean  and 
time.  Internal  waves  must  satisfy  boundary  conditions  at  the  surface  and  bottom  of  the 
ocean  (analogous  to  standing  waves  on  finite  strings).  A single  time  harmonic  component  of 
these  waves  has  the  components  of  displacement,  u,  v,  w.  The  vertical  component  is  assumed 
to  be  separable  in  the  coordinates, 

w(x,  y,  z,  t)  = W(z)  F(x,  y)  e'itJt 

The  satisfaction  of  the  boundary  conditions  restricts  W to  specific;  modes,  W(j,  k,  z)  in 
which  j is  the  mode  number,  k is  the  “horizontal  wave  number”  , and  Ah  is  the  horizon- 
tal wavelength  of  internal  waves.  The  dimensions  of  k are  cycles  per  kilometer.  An  exam- 
ple of  W(z)  is  W(3,  0.1,  z)e‘<kx  - w (3>  0 •1>t  where  k = 1/10.  This  is  shown  in  Fig.  3a. 

This  pattern  of  displacement  is  moving  to  the  right  (x  - direction),  Fig.  3b:  At  any  fixed 
distance  x = const.  The  pattern  of  Fig.  3 oscillates  with  time-regulated  phase  cut  so  that  for 
odd  multiples  of  it  the  displacement  arrows  shown  reverse  direction.  These  modes  are  key 
elements  in  the  analysis.  The  basic  premise  is  that  5c(z)  can  be  expanded  in  modes  W(j,  k,  z). 
The  expansion  coefficients  are  written  as  random  variables  A (j,  k), 

6c  internal  waves  = 6jur  = Ci  — ~ (Re  A + 3mA) 

•s/2 

A=  £ N2(z)  A(j,  k)  W(j,  k,  z)  exp  [i  (kx  - cu(j,  k)t)] 
j,  k 

Here  it  is  assumed  that  the  internal  wave  is  traveling  parallel  to  the  sound  wave  (x-direction). 

Attention  is  focused  on  amplitudes  A.  They  are  taken  to  be  complex  Gaussian  ran- 
dom variables.  For  any  mode  j the  power  spectrum  of  internal  waves  is  <|A(j,k)|2>,  and 
can  be  plotted  as  a function  of  horizontal  wavenumber  k of  the  internal  wave.  Models  of 
internal  wave  spectra  have  been  constructed  on  experimental  bases.  An  oft-quoted  model 
is  that  of  Garrett  and  Munk:  they  take. 


<A>  = 0,  <A(j,  k)  A*(j\  k')>  = 5jj.,  Skk',  02  H(j)  B(j,  k) 

H(j)  = (^j)2  ; B(j’  k)  = n (k2  +lj2)i/2  ;1i=T  n!-  J 

in  which  j = 0, 1,  2. . N0  = const.  = 3 cycles  per  hours,  CJj  is  the  inertial  frequency  of  the 
earth’s  rotation  (Coriolis  effect)  and  taken  to  be  one  cycle  per  day,  and  B = const  = 1000  m. 
The  spectrum  is  normalized  such  that 

E f H(j)  B(j,  k)  dk  = 1; 


y*00 

/ B(j,  k)  dk  = 1 

^OO 


In  numerical  work,  for  each  choice  j,  k the  complex  random  expansion  coefficient 
I A(j,  k)  I e‘*  is  generated  as  follows:  | A | is  a random  number  chosen  from  a population  in 
which  the  probability  of  occurrence  of  magnitude  lying  between  |A|  and  | A + dA|  is  Ray- 
leigh distributed,  and  <f>  is  a phase  selected  randomly  from  the  range  0 to  2ir  with  equal 
probability.  A computer  algorithm  is  available  for  calculating  complex  A according  to  this 
model.  Typical  curves  are  shown  in  Fig.  4.  Once  being  calculated  and  available,  the  numbers 
5cjW  are  added  to  CCp  to  form  the  total  c. 


C (r,  t)  - Cep  (z)  + 5c  iw 

5 c iw 

Sc/c  = e te-,)  - (1  - Vi)  1 + — “ 

Ci 

Thus  a randomized  sound  speed  profile  can  be  generated  at  any  plane  irl  = x = const,  at  any 
time  t,  as  a function  of  depth  z in  any  part  of  the  ocean,  recognizing  that  such  a construc- 
tion is  a model. 
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Figure  4 

Squared-Magnitude  Expansion  Coefficients 
| A|2  vs  k with  j = 1,  2 ...  as  Parameter 
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Assume  that  one  realization  of  a random  sound  speed  profile  caused  by  the  entire 
spectrum  of  wavenumbers  of  the  internal  waves  at  a particular  instant  of  time  is  made 
available  by  numerical  computation  as  described.  This  quantity  is  next  inserted  into  Eq. 
(5)  (i.e.,  the  parabolic  equation).  It  is  then  directly  found  that  the  solution  for  the  field  p 
which  can  generate  the  form  in  Eq.  (5),  has  itself  a form. 


P = V?  expi  (kor~  wt>- 


kD  = — , cx  = const. 


Thus  the  function  'k  (r,  z,  t)  also  satisfies  the  parabolic  equation  and  it  is  the  entity  to  be 
determined  from  the  differential  equation.  To  solve  for  'k  via  Eq.  (5)  Tappert  and  Hardin 
use  the  “split-step-Fourier”  method.  This  is  a numerical  “marching  algorithm”,  meaning 
that  given  'l'(r,  z)  one  is  required  to  find  >P(r  + dr,  z).  The  rationale  of  their  procedure  is 
this:  as  a result  of  increasing  range  r by  dr,  a phase  increment  exp  (-ikG— dr)  occurs  in  the 
wavefront.  The  product  of  the  field  and  phase  increment,  namely  ^(r,  z)C exp  (-ikc  ^ dr) 
is  then  Fourier  transformed  with  respect  to  z,  using  the  transform  pair  k2z,  where  kz  is  the 
z-component  of  the  field  wavenumber.  In  “kz  space”  the  term  32/dz2  adds  its  phase, 
namely  exp  i (kz2/2ko)dr,  which  multiplies  the  previously  obtainecf  transform . This  prod- 
uct is  finally  retransformed  back  into  real  space;  the  result  is, 


I (r  + dr,  z)  = 


kz2 

2k0  dr 


-i  kQ  dr  1 

|_e  °c  *(r,z)Jj 


In  actual  numerical  work  the  field  of  a single  point  source  at  Z = 1000  m was  used  to  give  a 
starting  value  4*  (0,  z).  512  increments  in  depth  z were  used  to  numerically  evaluate  the 
Fourier  transforms  (using  FFT).  The  increment  in  range  dr  was  1/2  km.  The  time  record 
allotted  for  calculation  was  128  h during  which  6 c/c  varied  randomly  according  to  the 
Garrett-Munk  model. 


Two  types  of  receivers  were  used  in  the  model  calculation.  In  the  first  a single 
hydrophone  received  several  acoustic  rays  simultaneously,  resulting  in  a model  of  multi- 
path  interference.  Fluctuations  of  a 100  Hz  signal  over  a period  of  128  hours  showed  5 to 
30  dB  changes  in  the  intensity  of  the  received  field.  In  the  second  type  of  receiver  a verti- 
cal array  of  hydrophones  700  meters  long,  centered  on  the  channel  axis,  was  used.  A 
beamformer  in  the  signal  processing  separated  out  single  rays  (or  arrivals).  The  result  was  a 
considerable  reduction  in  fluctuations  of  intensity  of  received  pressure. 

Paper  L2.  Path  Integrals  for  Waves  in  Random  Media,  R.  Dashen  (Inst,  for 
Advanced  Studies,  Princeton). 


Dashen  takes  Q - 0 in  Eq.  (1),  meaning  there  are  no  true  sources  in  the  volume  of 
ocean  under  consideration.  He  also  makes  the  parabolic  approximation, 


V2 

2k 


£{ 


kp(x,  t) 


a2 

2y2 


Here  z is  the  range  coordinate,  x is  the  depth  coordinate,  and  the  plane  xy  is  normal  to  z 
(direction  of  propagation).  In  the  plane  z = 0 (that  is,  at  zero  range)  there  is  a distribution 
of  field  p(x,  y,  z = 0)  due  to  some  true  point-source  external  to  the  volume  under  considera- 
tion. The  quantity  p is  the  perturbation  of  the  index  of  refraction  from  its  background 
value,  and  k = co/c  where  is  the  CW  frequency  of  the  acoustic  signal,  and  c is  the  local 
sound  speed. 
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The  random  variable  n has  a characteristic  length  L,  and  a characteristic  time  T over 
which  it  is  considered  effectively  fixed  in  value,  and  beyond  which  it  is  zero  in  magnitude 
and  purely  random  in  phase.  L is  analogous  to  the  size  of  an  inhomogeneous  blob  in  the 
ocean,  if  spherical.  If  the  blob  has  different  sizes  in  different  directions,  Lv,  Lh  are  used 
to  designate  characteristic  sizes  in  the  vertical  and  horizontal.  T is  analogous  to  the  dura- 
tion of  a “blob”  in  time. 

The  acoustic  size  kL  of  a blob  in  the  direction  of  propagation  is  taken  to  be  much 
greater  than  unity.  This  means  that  the  scattering  (or  diffraction)  of  a wave  X = ^ will  be 
sharply  in  the  forward  direction.  Similarly  the  time  required  for  a wave  of  frequency  oj 
and  speed  c to  pass  over  a distance  L is  taken  to  be  much  smaller  than  T.  This  means  the 
ocean  is  locally  “frozen”  during  passage  of  the  sound  wave.  Finally,  in  a range  r = R over 
which  the  blobs  follow  each  other  closely  there  will  be  R/L  of  them.  Assuming  the  angle 
of  deflection  of  the  normal  of  the  acoustic  wavefront  is  proportional  to  <p2  >1^2  it  is 
required  that  the  rms  multiple  scattering  angle  (=  <p2>1/2  R/L)  be  small.  This  is  a basic 
assumption  of  Dashen’s  treatment  of  fluctuations. 

Now  as  to  approach.  Dashen  concentrates  on  the  method  of  the  path  integral.  By 
all  accounts  this  is  an  intricate  tool.  As  a first  effort  to  elucidate  it  we  use  the  simpler  idea 
of  a 2-dimensional  “diffracting  screen”.  Consider  a 2-dimensional  antenna  lying  in  the 
plane  z = 0,  across  which  there  is  an  “aperture  function”  F (x0,  yo)eiwt  (optical  people  call 
this  “illumination”).  In  the  Fraunhofer  approximation  the  distance  from  an  elementary 
area  dx0,  dy0  in  the  aperture  to  the  far  field  point  P(xx , yx , zx ) is 


i*|:*z  + 2z  {(xo  -xi)2  +(yo  -yi)2} 


Thus  the  field  itself  is, 


P(xl,yi,zi)  = ^^)eik(z1-z o)j0T  dxQdyo  F(Xo>yo)e-2zi  K*0-*1)2  Mvo-yi)2] 

To  make  the  notation  symbolic  let  the  location  vector  in  the  xy  plane  be 
ri  = xi  K + Yi  ®y> 

ex,  ey  are  unit  vectors  in  which  i = 0,  1,  2.  . . .,  and  i designates  the  range  planes  z = 0,  z = 1. 
. . . .Then,  given  the  aperture  function  in  plane  z = 0,  one  can  write  symbolically^ 

ikeik(zl  “ zo)  f 

P(rx,Zil  r0,z0)  = -2n  (Zi  _Zq)  J F(r0)  G(rx  |r0)  dr0 


in  which  G (r  x |r  0 ) is  a symbolic  representation  of  the  phase  factor  in  the  Fraunhofer 
approximation.  In  a similar  way  the  field  at  plane  z2  can  be  determined  by  diffraction 
from  plane  z = zx . The  result  (again  symbolically  written)  is 


ik  eik((z2 -zi)  + (zi  _zo)) 

P (r2>  z2>  =7—  

27r  (zx -z0)  (z2  - zx) 


JJ drj>  d7x  F(r^)  X G (r^  | ^ ) G(?x  |^) 


Let  us  carry  this  out  to  N (equal)  increments  of  range:  then,  symbolically, 

p(rn,z„)  m J ••••drn-l  «£_i) 

In  the  limit  as  the  number  of  diffracting  apertures  (or  screens)  in  a given  range  increase  to 
infinity  one  arrives  at  the  infinite  dimensional  integral, 

p(r,z)  = J D [r(x,  y)]  J(r(x,y)) 

in  which  the  symbol  D specifies  an  infinite  dimensional  integration,  namely  a sum  over  all 
individual  screens  integrated  completely  over  x,  y in  each  screen.  Eq.  (6)  is  called  the  path 
integral  and  represents  a single  sweep  over  the  volume  between  planes  z = 0 and  z = Zn.  A 
sketch  of  source  (S)  and  receiver  (R)  situation  calculated  by  a path  integral  is  shown  in 


Path  Integral  Explained  as  a Sequence  of  Phase  Screens 


It  is  thus  seen  that  if  the  source  is  at  (r0,  z0)  and  the  receiver  at  (r3,  r3)  every  point  on 
every  conceivable  trajectory  between  planes  z0,  z3  is  in  the  volume  of  integration.  Thus 
Eq.  (6),  originally  a functional  integral  because  of  the  dependence  r(x,  y),  is  also  a path 
integral. 

The  particular  form  of  symbol  J depends  on  the  particular  form  of  £b,  £f  of  Eq.  (1). 

In  the  parabolic  approximation  treated  by  Dashen  the  diffraction  screen  is  taken  to  have  a 
complex  amplitude  distribution  as  an  aperture  function  (that  is,  F(x0,  y0)  is  complex), 
and  the  propagation  function  is  taken  to  be  the  phase  (and  amplitude)  associated  with  the 
Fraunhofer  approximation.  In  a homogeneous  medium  (meaning  the  index  of  refraction 
of  the  background  is  constant  and  the  fluctuations  are  dependent  only  on  distance  fx’1  — x0 1) 
the  phase  at  the  end  of  the  N increments  of  range  (or  N screens)  is 


- H (rj  + ez  Zj,  t) 


in  which  R is  the  total  range,  and  R/N  = Az  the  increment  of  range.  Note  that  this  is  a sum 
of  N increments  of  phase  caused  by  Fraunhofer  diffraction  of  a finite  aperture  with  aperture 
function  F(x0,  y0)  and  also  caused  by  the  random  index  of  refraction  p.  In  the  continuing 
limit  as  the  number  of  screens  between  source  and  receiver  becomes  infinite  the  phase  func- 
tion becomes 


L /*  i(*y 

I Jo  2 \dz/ 


p (r(z)  + gzz)  t 


Then  field  is  then. 


p(it  r*> , R,  t)  = ^ Jb  (r(z))  exp  j"-ik  J'  dz  - n (r(z),  t)j 


=J^0  Di  (r(z))  exp £ dz  (dl) 


where  by  definition. 


lim 

N-^-oo 


/ Di  (r*(z))  =y. dr^2 dr^ 


The  field  p is  stochastic.  The  amplitude  of  the  moments  of  this  random  variable  are 
defined  as 

<A2n>  = <(pp*)n> 

For  example,  the  mean  value  of  field  p in  a homogeneous  medium  is 

<p>  = <pp*>1/2 

- /d?„  F,r0)yD,  ,n.,)«P[A  £ (£)2  -*££ 


— 


p(|z-z’l)dzdz^D(r(z))  exp^^J' 


dz 


= &0  exp 


_ — <J)2 
. 2 


where 


/R  f 00 

p (Sz  z,  t)  dz)2  > =*  k2  R / p(|z|,0)dz 

J-co 


andg0  is  the  field  of  a homogeneous  medium  without  fluctuations  p.  Note  that  in  this 
approximation  any  vector  connecting  r(z)  and  r(z')  (that  is,  connecting  points  at  two 
range  planes)  is  considered  to  have  a magnitude  |r(z)  - r(z')  + (z  - z')2  I not  much  different 
from  | z - z'  |2.  This  once  more  is  the  parabolic  approximation  which  assumes  normals  to 
the  wavefronts  point  in  directions  always  close  to  the  z-axis  (axis  of  propagation). 


The  symbol  <£>  is  the  same  that  appeared  earlier  in  this  review.  It  is  a measure  of  the 
strength  of  fluctuations  caused  by  random  index  of  refraction  p.  By  itself  the  parameter  $ 
does  not  account  for  an  additional  deflection  of  the  normals  to  the  wavefront  caused  by  dif- 
fraction of  these  fronts  around  blobs  of  inhomogeneities  (in  p).  This  is  treated  by  a second 
parameter  SI,  or  the  squared  ratio  of  scale  size  L of  tho  inhomogeneous  blob  to  the  extent 
of  the  Fresnel  zone  ( ^XR)  as  measured  at  range  R, 

SI  = 6kL2  /R 


From  its  definition  it  is  seen  that  when  is  large  the  angle  of  diffraction  of  a wave  caused 
by  the  blob  is  very  small,  that  is,  geometrical  acoustics  applies.  Small  SI  therefore  corre- 
sponds to  large  diffraction  effects.  The  reciprocal  A = fi_1  is  sometimes  used  to  express 
the  same  effects,  namely  small  A means  small  diffraction  effects.  The  ratio  SI /$  (called  a) 
is  thus  a measure  of  both  fluctuation  strength  and  diffraction.  For  example  if  a«  1 
the  fluctuations  fo  the  field  p are  very  nearly  purely  random.  The  field  is  then  said  to 
saturate,  that  is,  the  field  is  determined  by  the  statistics  of  a random  walk  and  is  independent 
of  range.  However,  when  is  less  than  unity  and  is  greater  than  unity  the  vector 
representing  the  complex  field  retains  its  orientation  in  vector  space,  making  only  small 
random  deviations  about  its  orientation.  The  field  is  then  partially  saturated  or  unsaturated. 

The  parameter  a is  very  important  in  discussing  the  higher  moments  of  p.  To  see  its 
significance  we  again  examine  <p>.  In  the  phase  integral  method  there  is  only  one  volume 
of  integration  between  initial  and  final  planes  transverse  to  the  direction  of  propagation 
(Fig.  5).  The  actual  rays  from  r0  to  rlt  are  limited  by  the  small  scattering  angle  assumption 
to  the  angle  ~A/L  relative  to  the  geometric  acoustic  ray.  The  points  along  all  paths  that 
contribute  to  p will  extend  at  the  outermost  a transverse  distance  OR  = AR/L  from  the 
geometric  path.  But  AR  fL  = s/\R  ^/XR/L,  and  by  assumption  v'AR/L  «1.  Hence  only 
a very  small  volume  of  paths  (=  rays)  actually  contribute.  If  we  take  <1>  to  be  large  the  mean 
field  <p>  will  be  exponentially  small  as  the  range  increases.  Similarly  the  average  intensity, 
<pp*>  will  also  be  small  for  the  same  reasons. 


1 
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A more  important  second  moment  of  the  fluctuating  field  is  the  mutual  coherence 
function,  <p(r  lt  r 01,  tx)p*(r  2,  r 02,  t2)>.  The  path  integral  formulation  of  this  quantity 
is 

<p(l)p*(2)>  = ^ JJ  D(r1(z))D(72(z))exp  ^ (rx  '2(z)  - r2(z)2)dz  - V 

in  which  r ' = The  symbol  V expresses  a sum  of  three  quantities:  the  self  correlation 
of  field  along  path  1,  self-correlation  along  path  2,  and  cross-correlation  between  path  1 
and  path  2.  To  define  the  word  “path”  we  present  Fig.  6. 


PATH  1 


PATH  2 


Figure  6 


Sketch  of  Path  Integrals  in  Propagation 
of  Mutual  Coherence 


The  figure  shows  that  the  path  integral  formulation  of  <p(l)  p*(2)>  requires  a double 
volume  integration  over  the  same  range  between  the  initial  and  final  plane.  The  actual 
calculation  depends  on  the  type  of  field  equation  used.  Assuming  for  example  the  para- 
bolic approximation  one  arrives  at  the  result 


<p(l)  p*(2)>  = p*0(2)  exp 


'_D(1,  2)j 


in  which  Po(l),  p0(2)  are  solutions  for  the  field  in  a homogeneous  deterministic  medium, 
and  D (1,  2)  is  the  structure  function,  where  for  large  <f>, 


L,  2)  = 4>2  C 


Or,  -r2)2  +(r01  -r02)2  + (r,  -r2)(r01  -r02) 


t2  ! 


It  is  seen  that  only  when  r j~r  2,  roi~*02  the  mutual  coherence  has  significant  magni- 
tude. Thus  the  value  of  the  double  volume  integral  is  significant  only  when  the  paths  are 
close  together. 


The  fourth  moment  of  the  random  field  is  also  important, 


: p*(4)  p(3)  P*(2)  p(l)  > = (2ir)-4  J D(ra (z))  D(r2  (z))  D(r3  (z))  D(r4  (z)) 

i-1  -« 


in  which  M is  a matrix  of  all  self  and  cross  correlations  of  four  paths.  A sketch  of  the  four 
paths  is  shown  in  Fig.  7.  The  crucial  parameter  in  calculating  this  quadruple  volume  inte- 
gration is  the  distance  d between  pairs  of  paths  expressed  in  units  of  scale  length  L. 
Roughly 


d $R  1.  _ 

~L  6kL2  a ” « 

Thus  when  a is  small  (that  is,  when  the  strength  of  fluctuations  is  large)  the  separation 
between  pairs  of  paths  is  large  relative  to  L,  with  the  result  that  the  pairs  are  uncorrelated. 
Thus  the  quadruple  moment  becomes  the  sum  of  products  of  second  moments, 

< p*(4)  p(3)  p*(2)  p(l)  > a < p*(l)  p(3)  > < p*(2)  p(4)  > 

+ < p(2)  p*(3)  > < p(l)  p*(4)  > 

This  shows  that  for  4>  large  all  the  moments  belong  to  a Gaussian  distribution. 


Calculated  by  Path  Integrals 
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CLASSIFICATION  OF  REGIONS  IN  Ji<!>  SPACE 


By  definition  the  correlation  for  index  of  refraction  in  a homogeneous  medium  is 
<p(x,  t)  p(x,'t)  > = p(|x  - x'|,  t - t') 

Let  p (q)  be  the  dimensional  fourier  transform  of  p,  that  is,  p(q)  is  the  power  spectrum  of 
p in  the  wavenumber  q.  Assume  that  p(q)  goes  as  |q|-2  “P  where  for  large  q,  4 < p < 1.  To 
account  for  time  variations  suppose  that  the  medium  p is  time-dependent  so  that  plane 
waves  propagating  through  it  have  a dispersion  relation  cj  ~ k6 . Then  the  Fourier  transform 
of  the  second  derivative  of  is, 

p (|q|)  = const  X |q|  -(2  + P “ 

Four  parameters  are  now  available  for  classifying  signal  statistics,  4»,  S2,  p and  5.  They  will 
be  applied  to  collections  of  paths  that  connect  source  to  receiver.  Any  group  of  paths  that 
satisfy  a perturbed  ray  equation  of  the  form 


d2  r^z) 
dz2 


92  92  \ -*■ 

*3^  )<‘W*>  + Vl  ■ 0 


will  be  called  Fermat  paths,  or  nearly  geometric  paths.  If  in  a given  space  and  time  varying 
medium  there  are  meaningful  multiple  Fermat  paths  all  lying  within  a correlation  length  L 
of  each  other,  the  regime  of  fluctuations  of  the  field  p in  the  plane  will  be  called  partially 
saturated.  If  the  Fermat  paths  are  spaced  more  than  L units  apart  the  regime  will  be  called 
fully  saturated.  The  exact  location  of  regimes  in  the  4>£2  plane  depends  on  the  magnitude 
of  p.  The  following  table  lays  out  the  regime  for  the  case  when  temporal  changes  in  p can  be 
neglected. 


Case 
2 < p < 4 


1 <p<2 


Partially  Saturated 
*>1,  *4/p/n>i,  4>/fi<l 


* > i,  $4/ p in  > i,  <*>2/p  /n<i 


Fully  Saturated 

4>  > 1,  <t>/n>i 


$>1,  4>2'P/n>l 


Now  let  us  suppose  that  the  medium  has  significant  time-varying  p,  and  choose  p < 2 for  the 
partially  saturated  region.  Then  the  statistics  of  acoustic  pressure  pa  are  given  by  the  table: 


Statistics 


Rayleigh 


Variation  of  Signal  with 

Space 

Time 

IIQER5353I 

p - 6 < 0 

p - 6 > 0 

Gaussian 

Gaussian 

Phase 

Wrapping 

Phase 
times  a 
Guassian 

.51 


The  meaning  is:  (1)  as  long  as  p < 2 the  intensity  I of  field  pa  is  Rayleigh  distributed,  that 
is  the  probability  distribution  is 

1 I 

P(I)  exp  

<I>  F <I> 

<In>  = n!  <I>n 

(2)  if  the  Fourier  power  spectrum  of  p falls  of  as  |q  |“*  where  £ < 4,  then  the  signal  variation 
with  random  spatial  variations  of  p is  Gaussian.  (3)  if  p ~ |q  |-t , £ < 4,  then  if  p is  varying 
in  time  such  that  plane  waves  are  dispersed  in  frequency  as  w ~ k6 , and  if  p - 5 < 0,  the 
track  of  the  tip  of  the  vector  pa  representing  the  complex  signal  will  have  a Gaussian  distri- 
bution in  angle  (that  is,  all  angles  can  occur  with  equal  probability).  The  tip  then  executes  a 
path  that  looks  like  a random  walk.  (4)  if  condition  (3)  holds,  except  that  p - 5 > 0,  the 
track  tip  will  wrap  around  the  origin,  changing  only  slightly  in  amplitude.  Random  walk 
and  phase  wrapping  of  the  complex  vector  of  the  field  as  shown  in  Fig.  8. 


IMAG  p 


IMAG  p 


REAL  p 


"RANDOM  WALK" 


REAL  p 


PHASE  WRAPPING 


Figure  8 

Sketches  of  Random  Walk  and  Phase  Wrapping  of  Phase  Angle  of  a Propagating  Wave 

(5)  When  the  field  pa  varies  in  frequency  the  statistics  of  correspond  to 

Pa  (w)  = Po  (w)  e'"*  X(w) 

in  which  \p  is  a real  Gaussian  random  variable  with  <\p>  = 0 and  <i//2>  = Rp(0)  X C-2,  X is 
an  independent  complex  Gaussian  random  variable  with  zero  mean,  and 

< X*(w)  X(w')  > = A (w'1-  w) 

Thus  in  the  partially  saturated  regime, 

<p*(w)p(w')>  1 , „ „ „ 1 . 

— = exp  --  (w  - w )2  R p c 2 A co  - w) 

P0(w  P0(w  ) 2 


The  quantity  p will  be  explained  later.  This  explains  why,  in  the  table,  one  finds  the  entry 
that  tiie  statistics  of  pa  for  variable  frequency  correspond  to  that  of  a field  of  phase  (i.e. 
exp  ic o\l/ ) multiplied  by  a Gaussian  (namely  X(oj)). 


So  much  for  the  partially  saturated  regime.  In  the  fully  saturated  regime  (4>  > 1,  a -*■  0) 
the  exponential  factor  leads  to  factors  of  unity  in  higher  order  statistics  of  pa.  Thus  in 
the  second  order  moment  of  pa, 


p^g)p,(— ;»_  A(UI. 

* (w)  Po  (w  ) 


The  absence  of  the  phase  factor  (see  in  the  above  table  in  the  entry  under  variation  of  the 
signal  with  frequency)  constitutes  the  principal  difference  between  partially  saturated  and 
fully  saturated  regions. 


INHOMOGENEOUS  AND  ANISOTROPIC  MEDIA 


The  random  index  of  refraction  p in  the  case  of  inhomogeneous  media  is  considered 
to  have  a mean  <p>  which  is  a function  of  range  and  depth.  In  addition  the  covariance 
p(=  <p(xi,  t)  p(x2,  t)  > - t)  ><  p(x2,  t')  depends  on  x = 1/2  (?x  + x2)  as  well  as 

on  the  shorter  length  vector  - x2 . Finally,  in  anisotropic  media  p depends  on  the  orienta- 
tion of  xx  - x2. 

We  begin  with  an  inhomogeneous  medium.  To  obtain  tractable  path  integrals  one 
requires  (1)  variations  in  p over  lengths  of  order  L4>  /J 2 to  be  negligible  (2)  variations  in 
<p>  over  L4>/J2  to  be  also  negligible  (3)  the  parabolic  equation  is  valid.  With  these  assump- 
tions the  statistics  of  the  field  pa  in  an  inhomogeneous  medium  are, 

<p>  = p0  e"1'2*  2;  <p*(2)  p(l)  > = p0(2)  p0(l)  e-l/2D 


Evidently,  the  form  of  these  equations  is  the  same  as  for  a homogeneous  medium  but  the 
explicit  forms  of  <l>  and  D are  different, 


<h2  = k2 


rK 

J p((f,  0;z)dz;  D = 2 


d(w(z),t!  - t2;  z)  dz 


in  which 


d(w,  t;  z)  = k2  [p(0,  0;z)  -p(w,  t;  z)] 

All  formulas  depend  on  the  quantity  p.  For  a homogeneous  medium  p is  the  transverse 
correlation  function  in  the  Markov  approximation,  i.e. 


In  an  inhomogeneous  medium  p is  more  complicated  but  retains  the  same  basic  definition. 
Explicit  forms  for  special  choices  of  p (x,  t)  are  available. 

Next  one  considers  anisotropic  media.  As  a measure  of  anisotropy  one  adopts  the 
angle  0O  which  is  the  ratio  of  the  small  dimension  of  an  inhomogeneous  blob  to  the  large 
dimension.  The  most  important  case  is  propagation  in  a channel  wherein  the  unperturbed 
I ray  makes  loops  and  wherein  the  long  axis  of  the  inhomogeneities  is  parallel  to  the  channel 

axis.  This  is  shown  in  the  figure.  (Figure  9) 


A 


Figure  9 

Sketch  of  the  Propagation  of  a Ray  in  Anisotropic  Media 

The  scattering  of  the  wave  is  strongest  at  the  “turning  points”  of  the  ray  pa  h A,  and  weakest 
in  between,  B.  To  understand  the  cumulative  effects  of  scattering  we  consider  Fig.  10. 
in  which  we  show  a wavefront  WF  passing  through  a homogeneous  random  medium  from 
Zj  to  z2.  We  assume  p is  well  correlated  in  the  plane  z = const.  However,  in  the  direction  ^ 
of  propagation  the  correlation  is  a delta  function  5(z  - z'):  this  means  there  is  no  ‘memory” 
of  previous  states  of  p with  change  in  the  z-coordinate.  The  wavefront  itself  moves  nearly 
unaffected.  Hence  the  coupling  between  random  field  pa  and  random  index  of  refraction 
is  the  same  as  that  of  two  statistically  independent  entities  viz  <Pa/i>  = <Pa><#i>-  A 
medium  which  is  thus  delta  function  correlated  in  the  direction  of  propagation  generates  a 
field  p which  is  Markov  approximated.  In  isotropic  media  the  scattering  does  not  depend 
on  angle  of  the  normal  of  the  wave  front  relative  to  some  reference.  In  the  Markov  approxi- 
mation (still  considering  isotropic  media)  the  multiple  scattering  angle  (i.e.  its  total  accumu- 
lation) is  small.  The  picture  is  that  of  a wave  moving  nearly  in  the  same  direction,  whose 
normals  bob  back  and  forth  randomly,  making  small  angles  to  the  direction  of  propagation. 

In  anisotropic  media  scattering  of  the  field  pa  is  highly  dependent  on  angle  of  incidence  as 
soon  as  the  angle  is  of  order  0O.  Hence  the  scattering  is  reinforced  in  special  directions  and 
the  wavefront  normals  “remember”  previous  history.  Hence  the  Markov  approximation  fails 
at  a shorter  range  than  for  isotropic  media.  But  suppose  the  medium  is  anisotropic  and  at 
the  same  time  inhomogeneous  (as  pictured  in  Fig.  9).  The  scattering  effectively  occurs  only 
at  the  turning  points  (points  A).  At  each  turning  point  we  assume  the  Markov  approximation 
holds.  Generally,  the  turning  points  are  separated  by  more  than  a coherence  length.  Hence 
the  effect  of  previous  scatterings  will  be  equivalent  to  a random  modulation.  Thus  the 
Markov  approximation  will  be  valid  at  long  ranges. 
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To  this  point  it  has  been  shown  that  formulas  for  <pa>  in  inhomogeneous  aniso- 
tropic media  have  the  same  form  as  those  of  homogeneous  insotropic  media,  but  requiring 
only  change  in  the  definition  of  <l>  and  D.  In  the  same  manner  it  can  be  shown  that  by  re- 
! defining  ft  and  4>  the  table  of  statistical  regimes  of  the  field  pa  in  the  partially  saturated 

regime  have  the  same  form  but  that  the  statistics  of  pa  in  the  saturated  regime  where 
frequency  changes  are  in  question  is  Gaussian  with  the  added  phase  described  before.  The 
great  importance  of  the  parameter  ft  is  once  again  brought  out  by  this  redefinition.  It  is 
this:  ft  measures  the  rate  at  which  the  phase  of  the  oscillating  factor  in  the  path  integral 
formulation  of  the  field  pa  varies  as  the  rays  of  the  path  move  away  from  the  unperturbed 
ray.  In  different  terms,  ft  is  a measure  of  the  phase  change  required  to  move  a ray  a dis- 
tance L away  from  the  unperturbed  ray.  A specific  form  of  ft  for  the  case  of  anisotropic 
inhomogeneous  media  has  been  derived.  Together  with  the  similarly  modified  form  of  A 
(co  - to')  it  can  be  shown  that  the  statistics  of  the  fully  saturated  regime  are  completely 


Gaussian. 


Paper  L3  “AFAR  Results:  Comparison  of  Experiment  and  Theory.” 

A.  W.  Ellinthorpe  and  R.  Dashen 

(Reviewer’s  Note:  An  appreciation  of  this  paper  requires  a rather  complete 
knowledge  of  a set  of  environmental  and  acoustic  experiments  originally 
conducted  by  Ellinthorpe  in  1975.  A summary  of  these  results  presented 
below). 

The  Azores  Fixed  Acoustic  Range  (AFAR)  is  an  underseas  basin  some  8500  ft  deep 
located  between  three  submerged  sea  mounts  at  the  comer  of  a loose  triangle  in  the  geo- 
graphic area  of  the  Azores  some  1800  ft  to  2500  ft  below  sea  level.  The  underseas  facili- 
ties of  this  range  (called  Echo,  Sierra,  November  and  Bravo)  are  located  on  this  mounts 
(Fig.ll).  The  distances  between  them  are:  Echo  to  Sierra,  35000M;  Echo  to  Bravo,  2990 
meters;  Echo  to  November,  about  34000M.  Echo,  Sierra,  and  November  have  receiver 
masts  mounted  with  directional  receivers  which  receive  on  30°  beams  at  600  HZ  and  8° 
beams  at  3500  HZ.  Bravo  is  an  omnidirectional  hydrophone  about  1000  ft.  deep.  Of  all 
the  facilities  only  Echo  can  operate  in  the  transmit  mode.  A typical  receiver  mast  with 
three  directional  antenna  is  shown  in  Fig.12.  Typical  theoretical  ray  paths  based  on  histori- 
cal temperature  data  are  shown  for  Echo  to  Sierra,  Fig.13  and  Echo  to  Bravo,  Fig.  14. 

An  acoustic  experiment  was  conducted  as  follows:  A CW  signal  at  400  HZ  was 
radiated  by  Echo  and  received  by  the  (directional)  receiver  at  Sierra,  and  by  the  (omni- 
directional) receiver  at  Bravo.  A 5-HR  record  of  reception  is  shown  in  Fig.15.  The  received 
signal  at  Sierra  is  a recognizable  sinusoid,  that  is,  a rotating  vector  whose  relative  amplitude 
on  dB  shows  slow  variations  of  up  to  40  dB  over  the  5 hour  period.  The  phase  of  this  vec- 
tor as  measured  against  a local  clock  changed  by  some  0.009  - 0.001  = 0.008  sec  at  400  Hz 
(or  6.4  7T  radians).  The  same  CW  signal  at  Bravo  shows  very  rapid  amplitude  fluctuations 
of  some  1 dB  superimposed  on  a slow  fluctuation  of  some  5 dB.  The  phase  change  at 
Bravo  relative  to  a local  clock  was  about  0.0088  - 0.0048  = 0.004  sec.  or  3 ir  radians. 

A second  experiment  repeated  the  first  at  4400  Hz.  A 5 hour  record  is  shown  in 
Fig.  16.  This  time  the  signal  amplitude  received  at  Sierra  showed  very  rapid  fluctuations  of 
from  0 dB  to  40  dB  while  the  phase  varied  by  about  0.004  sec  or  35.2  n rad.)  relative  to  a 
local  clock.  The  received  amplitude  at  Bravo  showed  rapid  severe  fluctuations  of  some  5 
to  30  dB  over  the  period  and  a phase  change  of  0.004  sec  or  35.2  n radians  relative  to  a 
local  clock. 

In  order  to  explain  these  received  signals  a set  of  environmental  data  was  taken, 
which  was  then  reduced  to  statistics  of  sound  speed  profiles  as  functions  of  space  and  time 
ill  the  AFAR  basin,  as  well  as  data  on  buoyance  frequency.  Typical  calculations  in  the 
signal  processing  were:  (1)  the  variation  of  sound  speed  along  a horizontal  track  of  10 
miles  was  Fourier  transformed  to  a spatial  power  spectrum  versus  wavenumber  “k”  in  units 
of  cycles  per  meter.  A typical  plot  is  shown  in  Fig.17  in  which  the  noticeable  feature  is 
the  breakpoint  between  two  regimes.  (2)  the  spatial  power  spectrum  for  a 90  mi  track  is 
shown  in  Fig.18.  Note  the  three  breakpoints.  (3)  the  coherence  in  time  of  two  instruments 
measuring  sound  speed  in  the  horizontal  plane  separated  by  1100  meters  versus  wavenum- 
bers k was  measured  and  is  shown  in  Fig. 19.  Also  the  coherence  of  two  instruments  in  the 
vertical  plane,  Fig.  20.  (4)  Typical  sound  speed  profiles  are  shown  in  Fig.  21.  (5)  The 


CONTOUR  CHART  OF  BOTTOM  TOPOGRAPHY  OF  AFAR  Figure  11 


Fig.  12  —Physical  Size  and  Proportions  of  AFAR  MAST 


Fig.  14  —Direct  Path  and  One  Reflected  Path  of  AFAR  Echo-Source  and  Bravo  Hydrophones 
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LOG  CPM 

Spectrum  of  Sound  Speed  Fluctuations  Along  a 90  mile  Track 


Temporal  Coherence  of  Sound  Speed  Fluctuations  as  a Function  of  Spatial  Wavenumber 


10  LOG  (CPM) 

Coherence  of  Two  Sensors  in  Vertical  Plane  as  a Function  of  Spatial  Wavenumber 


power  spectrum  of  random  variations  in  vertical  sound  speed  profile  with  mean  removed  is 
shown  in  Fig.  22.  (6)  Power  spectra  of  water  currents  (similar  to  internal  wave  spectra) 
were  obtained.  These  are  shown  in  Figs.  23,  24.  (7)  A composite  spectrum  of  horizontal 
and  vertical  fluctuations  in  sound  speed  profile  is  shown  in  Fig.  25. 

ELLINTHORPE’S  CONCLUSIONS 

(1)  Data  from  Echo-Sierra  and  Echo-Bravo  were  processed  to  find  the  “log  variance” 

<I2>  - <I>2  /<I>2of  the  fluctuations  of  intensity  of  the  arrived  signal.  These  are  shown  in 
Fig.  26.  From  it  one  concludes  that  the  log  amplitude  variance  of  pulse  amplitude  exhibits 
classical  saturation  at  these  ranges. 

(2)  The  vertical  spatial  correlation  length  behaves  approximately  as  fc~a(a  — 1/2)  for 
the  saturation  regime.  Fig.  27.  ( fc  is  the  carrier  frequency).  The  horizontal  spatial  correla- 
tion behaves  approximately  as  fc~a  (a  ~ 3),  Fig.  28. 

(3)  The  correlation  time  (time  to  decorrelate  to  e1  of  original  value)  behaves  as  /c_1 
for  complex  amplitude.  It  behaves  as  /c-1  for  amplitudes  in  the  saturation  regime,  and  f~a 
(a  < 1)  in  the  unsaturated  regime,  Figs.  29,  30. 

1 

(4)  The  structure  function  for  phase  of  the  received  signal  namely 

< {^(t  + t)  - 0(t)  - <<t>(t  + r)  - 0(t)  >}  2 > 

for  both  Echo-Sierra  and  Echo-Bravo  obeys  the  slope  r5/3 , Fig.  31. 

(5)  Comparisons  of  log  variance  of  intensity  of  AFAR  experiments  with  the  equa- 
tions of  Tatarski  and  Zachariasen  are  shown  in  Fig.  32.  This  figure  will  be  discussed  in 
some  detail. 

Visual  inspection  of  the  results  in  the  unsaturated  regime  ( o2^!  < 1.0)  leaves  little 
to  choose  between  equations  (1)  and  (2).  A least  square  fit  of  the  expression  o2  = cf°  to 
these  data  points  yeilds  an  exponent  a = 1.19  which  lends  support  to  equation  (1).  The 
fluctuations  are  obviously  saturated  at  all  frequencies  on  the  35  km  path.  Applying  the 
Student’s  t test  at  a 5%  significance  level  to  the  mean  of  these  saturated  observations 
leads  to  rejection  of  the  hypothesis  that  the  mean  = 7 r2/6.  The  analysis  which  predicts 
7r2/6  (Dyson  and  Munk)  assumes  that  the  propagation  is  characterized  by  several  inde- 
pendent paths  (n  > 3)  with  no  single  dominant  path,  leading  to  a Rayleigh  distribution 
of  received  signal  amplitude.  The  result  here  indicates  that  these  assumptions  may  not 
be  quite  true  or  that  noise  other  than  the  simple  Gaussian  variety  may  be  contributing 
to  the  log  variance. 

(6)  Direct  path  convariance  of  fluctuations  in  frequency  between  Echo  and  Sierra, 
namely 

<p*(w  + Aw)  p(w)> 

c(Aw)  

Po*(w)  Po(w) 

has  been  measured  and  appears  in  Fig.  33. 
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Fig.  22  — Power  Spectrum  of  Random  Variations 
in  Vertical  Sound  Speed  Profile 


FREQUENCY  CYCLES 


Fig.  23  - Power  Spectra  of  Water  Currents 


Correlation  Time  versus  Frequency 


f - FREQUENCY,  L = RANGE 

Fig.  32  — Comparisons  of  Theoretical  Log-variances  with  Experimental  Data 


'variance  of  Fluctuations 


COMPARISON  OF  JASON  THEORY  WITH  AFAR  EXPERIMENTAL  DATA 


At  the  completion  of  the  expositions  of  Flatt^,  Ellinthorpe  and  Dashen,  the  final 
portion  of  the  invited  papers  was  reached.  In  it  Dashen  presented  a detailed  comparison 
between  predictions  of  the  Jason  Group  and  actual  measurements  at  the  AFAR  facility. 

A summary  of  these  comparisons  is  presented  below. 

A single,  fully-refracted  path  of  the  Echo  to  Sierra  stations  over  the  35  km  range  was 
used  to  check  theory.  The  first  task  was  to  obtain  a sound  speed  profile  which  was  the 
sum  of  a background  profile  plus  fluctuations.  Profiles  of  sound-speed  measurements  were 
made  between  depths  of  240  m to  840  m.  Upon  averaging  these  to  obtain  a background 
profile,  it  was  found  that  the  average  did  not  agree  with  the  canonical  model  of  Munk. 

The  experimental  profile  was  then  fitted  to  a sixth-order  polynomial,  and  the  result  was 
used  in  the  analysis  as  the  deterministic  background  sound  speed  profile. 


A second  task  was  to  obtain  the  fluctuating  part  of  the  sound  speed,  </i2(z)>.  It 
was  assumed,  to  start,  that  this  part  was  caused  by  internal  waves.  To  calculate  it  from  the 
internal  wave  model  required  measurements  of  the  buoyancy  frequency  m(z)  as  a function 
of  depth.  Actually,  temperature  and  salinity  measurements  were  first  made,  and  then  from 
these  a profile  of  buoyancy  frequency  versus  depth  was  constructed.  It  roughly  approxi- 
mates the  form  N0e  z/B,  Furthermore,  in  addition  to  m(z)  one  needs  the  gradient  of  the 
sound  speed,  9Cp/9z  = + CX7A,  where  C is  the  background  sound  speed  profile,  Cx  = 

1500ms-1  and  7A  is  a stochastic  variable  from  the  theory.  These  quantities  were  all  assem- 
bled from  the  constructed  background  profile  and  from  the  internal  wave  model,  and  the 


prediction  was  made  that 


< jLt2  (z)  > = const.  X n (z) 


2 


Measurement  of  sound  speed  fluctuations  were  compared  with  this  prediction  for  the  pur- 
pose of  testing  it  for  depth  variation.  It  was  found  that  the  internal  wave  model  holds  well 
in  the  depth  region  between  240m  and  790m. 


Thus  there  was  on  hand  an  experimentally  certified  model  of  the  random  sound 
speed  profile.  One  was  then  prepared  for  comparison  of  the  acoustic  experiments  with 
theory. 


First  the  CW  experiments.  Pulses  centered  at  several  frequencies  with  widths  of  130 
msec  were  used  to  simulate  CW  transmission.  The  received  signal  was  a CW  sinusoid  modu- 
lated by  the  random  properties  of  the  medium.  A first  question  was,  how  severely  did  the 
random  sound  speed  affect  the  phase  and  amplitude  of  the  received  signal,  visualized  as  a 
rotating  phasor,  Ae'0?  The  answer  required  calculation  of  the  two  parameters,  4>  and 
j2-l(=A),  as  explained  earlier  in  the  expository  discussions  written  above.  The  calculation 
of  4>  depends  on  the  numerical  evaluation  of  an  integral  of  </j2(s)>/1(s)  along  the  refrac- 
ted path,  where  /2  (s)  is  a theoretical  quantity  known  from  the  theory  of  Munk  and 
Zachariasen.  The  result  of  the  integration  was 


* 


in  which  a was  the  acoustic  frequency.  The  calculation  of  n1  = A was  made  from  the 
simple  formula  A = R/6kLv2.  in  which  R is  the  range,  K =gj/c,  and  Lv  is  the  vertical 
correlation  length  of  the  random  inhomogeneities  in  index  of  refraction  in  the  ocean.  The 
result  was 

25 

A = — ; HZ,  Lv  = 240  M 
a y 

Thus  the  all-important  product,  <f>A,  had  the  value  0.63,  and  the  product  A<J>2  was Hz. 
At  the  frequencies  in  question  (412  to  4672  HZ)  this  combination  in  the  A4>  plane  indi- 
cated that  the  regime  of  statistics  of  Ae'0  was  partially  saturated. 

Upon  certification  of  this  regime  it  was  then  predicted  that  the  amplitude  statistics 
of  Ae'0  would  be  approximated  Rayleigh  distributed  (rather  than  log  normal).  This  was 
experimentally  verified.  Fig.  34.  In  addition  it  was  predicted  that  the  variance  of  log-inten- 
sity distribution  would  follow  the  classical  plot,  that  is,  increase  with  A<J>  up  to  an  asymp- 
totic value,  calculated  in  this  model  to  be  7r2/6.  This  prediction  was  verified,  Fig.  35.  A 
separate  issue  was  how  well  the  higher  moments  of  received  intensity  agreed  with  the 
modified  Rayleigh  distribution  derived  by  Dashen,  namely 

<In>  = n!  <I>n  (1  + - n(n-l)7 

2 

This  prediction  was  verified  at  the  frequency  4671  Hz,  but  no  verification  was  possible 
with  data  at  the  lower  frequencies  because  of  insufficiently  long  experiment-time-records. 

Another  group  of  comparisons  between  theory  and  experiment  was  concerned  with 
temporal  fluctuations  of  the  two-point  correlations  of  the  field  i// . Dashen ’s  theory  pre- 
dicts that 

( a tl 

< (t  + At)  ^*(t)  > e 2 

and  that  for  small  At  time  separations  one  should  have  D— $2(At)2.  Thus  one  requires  an 
environmental  parameter  ^ to  be  constructed.  Actually  4>  was  obtained  from  the  inter- 
nal wave  model,  and  <t>  itself  was  calculated  in  the  manner  discussed  above.  A set  of  pre- 
dictions of  the  various  2-point  temporal  statistics  of  is  shown  in  Fig.  36.  The  experimen- 
tal points  verify  the  predictions  very  well,  incidentially  certifying  that  the  regime  is  par- 
tially saturated,  as  predicted  by  other  parameters.  In  particular  the  time  variation  of  the 
phase  of  \p  shows  phase  wrapping  as  predicted. 

A final  group  of  comparisons  was  concerned  with  the  temporal  statistics  of 
arriving  trains  of  short  pulses.  These  pulses  were  centered  at  3200Hz  and  were  1.28  ms 
wide.  As  they  traveled  through  the  medium  they  were  randomly  dispersed  in  frequency 
so  that  the  intensity  on  arrival  varied  randomly  with  time.  The  first  order  statistic  of 
intensity,  namely  the  mean  value  <I(r)>  was  measured,  and  then  compared  with  pre- 
dictions based  on  the  AFAR  model  of  the  sound  channel.  The  agreement  with  experi- 
ment is  good  over  the  time  interval  0 to  40  msec.  Fig.  37.  Two  theoretical  models  are 
used  to  compare  intensity  distribution  experimentally  obtained.  In  one,  the  quantity 
<I>  is  considered  proportioned  to  exp  (-a/4>)  Q(a),  in  which  a is  the  acoustic  fre- 
quency; a = At/t?0  is  a parameter  describing  the  frequency  spread  Ao  relative  to  a 
reference  bandwidth  rj0  (=  40HZ)  provided  by  the  model;  and  Q(a)  is  a factor  derived 
from  an  internal  wave  model.  In  the  other  model,  <I>  is  proportional  only  to  exp 
(-  ct/4>).  Fig.  37  shows  that  the  inclusion  of  Q(a)  brings  experiment  and  theory  much 
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closer  to  agreement  on  the  variation  of  <I>  with  time.  The  second  statistic  checked  is 
the  auto  correlation  (in  time)  of  the  received  intensity  < I(t  + r)I(t)>.  Theory  predicts 
that  the  Fourier  transform  of  this  quantity  should  be  proportional  to  the  Fourier  trans- 
form of  |Q(a)|2  obtainable  from  the  internal  wave  model.  Comparison  with  experiment 
shows  good  agreement.  Fig.  38. 

The  success  in  these  comparisons  is  commendable.  However,  it  is  recognized  that 
the  isolation  of  a single  path  in  these  AFAR  experiments  excludes  the  serious  difficulty 
of  multipath  interference.  This  problem  is  still  under  investigation. 
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